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AERODYNAMIC SYMBOLS 


1. GENERAL 

m Mass 

t Time 

V Resultant linear velocity 

Q Resultant angular velocity 

0 Density, o relative density 

y Kinematic coefficient of viscosity 

R ‘Reynolds number, R = /V/y (where / is a suitable linear dimension) 
Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. 
For air under these ( pe = 0-002378 slug/cu. ft. 

conditions | » =.1:59 « 10=4 sq. ft. /sec. 

The slug is taken to be 32-2 Ib.-mass. 

OL Angle of incidence 

é Angle of downwash 

5 Area 

b Span 

C Chord 

A Aspect ratio, A = 0?/S 

1 Lift, with coefficient C, = L/$oV?S 

D Drag, with coefficient Cp = D/teV°S 

y Gliding angle, tan y = D/L 

L Rolling moment, with coefficient G = L/4oV20S 

M Pitching moment, with coefficient C,, = M/teVcS 

N Yawing moment, with coefficient (, = N/}eV*0S 


2. AIRSCREWS 
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Kevolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient ky = T/en?D* 
Torque, with coefficient kg = Q/en?D° 
Hficiency, 7 = TV/P = J[kp/2aky 
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Summary.—Three methods of obtaining the frequencies of the overtone vibrations of elastic 
bodies are described. The first consists in the direct application of the Rayleigh or Galerkin 
methods. Experience shows that in order to obtain a good approximation to the frequency of 
the 7 overtone it is necessary to use (vy + 2) independent displacement functions or degrees of 
freedom, and the labour becomes heavy when 7 exceeds 1. This labour may be reduced by 
supplementing the Rayleigh or Galerkin methods by the matrix solution given in §3 of the text. 
However, a method depending on an application of the orthogonal properties of the modes is 
developed which allows the overtone frequencies to be obtained with satisfactory accuracy when 
only two degrees of freedom are used (§4). It is also shown that overtone frequencies can be 
obtained by a modification of the well known iterative process. The ordinary iterative method 
fails to yield the overtones on account of the phenomenon of regression to the fundamental 
(see §5), but this difficulty is avoided by making use of the orthogonal properties of the modes, 
and without any assumptions regarding the positions of the nodes. The various methods are 
illustrated by worked examples. 
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$1. Introduction.—There is no difficulty of principle in obtaining the frequencies of the overtone 
vibrations of elastic bodies by means of the well known Rayleigh-Ritz and Galerkin methods. 
However, it is necessary according to experience to use not less than (vy + 2) independent 
displacement functions when it is desired to obtain the frequency of the 7 overtone with 
satisfactory accuracy. The labour entailed in solving (v7 + 2) simultaneous equations containing 
a frequency parameter (characteristic number) by the ordinary methods when 7 exceeds | is 
excessive, and it is therefore desirable to find some method or device by which the labour may be 
reduced. One such device consists in solving the set of equations yielded by the Rayleigh or 
Galerkin methods by the application of the iterative matrix solution already given elsewhere * °, 
and an example of this is given in §3. The advantage of the matrix process becomes more 
pronounced as the number of equations in the set increases, but the method cannot be used 
until the rules for the basic arithmetical operations with matrices have been mastered ; these 
are in fact simple and can easily be learned by those having no knowledge of the mathematical 
theory of matrices. However, it was thought desirable to formulate an easy process for obtaining 
overtone frequencies without recourse to matrices, and this has been done by making use of the 
orthogonal properties of the modes. 


A displacement of an elastic body which is consistent with the conditions of support, but 
otherwise quite arbitrary, can be analysed into a set of displacements each of which is a natural 
mode of vibration of the body with an amplitude appropriate to the particular assumed 
displacement. But the amplitude of a selected mode of vibration will be zero when the assumed 
displacement renders a certain integral containing the displacement function of the selected mode 
zero;.and when this is so the assumed displacement is said to be orthogonal to the selected mode. 
Suppose now that it is desired to find the first overtone frequency, the fundamental having 
already been found with sufficient accuracy (and it is important to observe that this can usually 
be done without excessive labour). Select a pair of independent displacement functions which 
are both orthogonal to the known fundamental mode, and employ these in the ordinary Rayleigh 
or Galerkin processes. The result will be a quadratic equation for the frequency parameter, and 
the lower of its roots will be an approximation to the value of the parameter for the first overtone, 
for no linear combination of the assumed displacement functions can contain the fundamental 
mode. The foregoing is a general description of a method which is explained in detail in §4, 
where the extension of the process to obtain the higher overtones 1s also given. It should be added 
that the restriction to the use of two orthogonal displacement functions is non-essential, but all 
the examples hitherto tested show that two independent functions are sufficient for obtaining 
the first and second overtones. 


A well known method for obtaining the fundamental frequency is that in which an arbitrary 
displacement is adopted, the inertia forces deduced from this and then taken as a static load 
system to yield a second or derived displacement for the body. The process is repeated with the 
derived displacement substituted for that first used, and continued until the type of displacement 
is sensibly constant ; this is the fundamental mode of displacement and the frequency can be 
calculated. This iterative solution fails for the overtones because of the phenomenon of regression 
to the fundamental (see §5), but a successful modification of it which makes use of the orthogonal 
properties of the modes has been devised, and is illustrated in §5. 


The choice of the basic set of displacement functions is a matter of great importance, and the 
form of the functions depends on the nature of the body at its boundaries. For example, in the 
case of a cantilever beam the section at the tip may be finite, or the beam may taper down to 
zero at the tip in a variety of different ways (for instance the tip may be rounded or pointed). 
The geometrical conditions to be satisfied by the displacement functions will in general be different 
in these cases, and this question is investigated in the Appendix, where suitable functions for 
various types of cantilever are given for the cases of flexural and torsional motion. 


It must be pointed out that it is useless to employ any refined process for calculating frequencies 
unless an adequate standard of numerical accuracy is maintained in the computations. The 


3 
numerical accuracy in the illustrative examples included in this paper is for the most part higher 


than would usually be considered necessary, but it was desired to obtain really precise comparisons 
of the several methods used, and reliable checks against known solutions. 


All the examples given in the paper refer to the oscillations of cantilever beams. The methods 
described are, however, applicable to elastic bodies in general. 


§2. Direct Use of the Rayleigh or Galerkin Methods.—2.1. General.—Galerkin’s method has been 
described and illustrated in R. & M. 17981, and it has been shown in R. & M. 1848" that this 
method necessarily yields results which are identical with those given by the application of 
Rayleigh’s principle in all problems concerning the free oscillations of elastic bodies. It is 
therefore unnecessary here to give more than a very brief summary of these two standard methods, 
and to state the equations to which they lead for the cases of the flexural and torsional oscillations 
of beams. 


The first step in the solution of a problem, whichever method is used, is to write down a sequence 
of functions, say Y,, Y,, Ys, etc., which satisfy all the end conditions or boundary conditions 
of the problem. The actual displacement is then assumed to be 


Veer gy 4\-(or Ga Yoh t= etc) 8 en . eh Shel sh) 


where the quantities g are constants whose ratios are to is HeRenied: The functions Y are 
conveniently polynomials of rising degree, or trigonometric functions of rising periodicity. 
Functions suitable for the treatment of a number of oscillation problems are given in the Appendix. 


In the Galerkin process the differential equation (or equations) of the particular problem is 
multiplied in succession by a set of functions which are usually the functions Y,, Yo, etc., 
themselves,* and then integrated over the entire oscillating body. This yields a set of linear 
homogeneous equations in the coefficients g, and their eliminant is an algebraic equation for the 
frequency of oscillation. 


In the application of the Rayleigh principle the equation of energy is used to express the 
square of the frequency in any pure mode as a function of the coefficients g. In accordance 
with the principle the frequency is stationary for small variations of the coefficients g from their 
true values. Hence the equation of energy can be differentiated with respect to each gq in turn, 
while the frequency is treated as constant. In this way a set of linear homogeneous equations 
in the coefficients g is obtained, and it can be demonstrated (as in R. & M. 1848) that these are 
actually identical with the equations yielded by the Galerkin method, provided that the 
displacement functions Y are the same for both methods. 


Suppose that a beam has flexural rigidity EI and mass per unit length m, where in general 
these will both be functions of the co-ordinate « measured along the axis of the beam. Then 
Rayleigh’s principle leads to a set of equations typified. by 


dy d°Yx Balla fagee 
[EI oat ax — | myY dx = 0 
eet Vell), ots De eel 
where w is 2x times the frequency of flexural oscillation, and the integrals extend over the whole 
beam. Since y is given by (2.1.1), the last equation is linear and homogeneous in the coefficients q, 
which may properly be regarded as the Lagrangian generalised co-ordinates of the beam. The 
corresponding typical equation given by the Galerkin method is 
(iv 12 12 ) dx — w { MYNX vax QO, ee gs = (Dele) 
but this can be shown to be identical with (2.1.2) by integration by parts and use of the end 
conditions. + 


* Rules for the choice of the multipliers are given in R. & M. 1848. 
+ It is assumed that the ends of the beam are free, simply supported or rigidly built in. 
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Next, for the case of torsional oscillation of a beam, let the angular displacement be 8, and 
assume that 


§ = 4,0, + 9.0, + etc., -.. e a tld) 


where ©,, @,, etc., satisfy the end conditions. Then if C is the torsional stiffness of unit length 
of the beam, and Ji is the polar moment of inertia of unit length about an axis through the centroid 
of the section, the typical Rayleigh equation is 


ao dQr 4 aid ak 
| CoAT dx — ot | J0Odx = 0 
(7a el ee eons Y 37 Zale 
The corresponding Galerkin equation 1s 
d ao ‘ a 
[ O £(C& Jax — o8[ J0Odx=0, .  .. .. @18) 


and the identity of (2.1.5) and (2.1.6) can be established by integration by parts. 


2.2. Examples of the Direct Use of the Rayleigh and Galerkin Methods.—A few examples will 
be given in order to show the degree of accuracy of the fundamental and lower overtone frequencies 
obtained by the use of varying numbers of independent functions (degrees of freedom) in the 
Rayleigh and Galerkin methods. 


Example 1. Flexural Oscillations of a Uniform Cantilever. 
use of the functions 


: Ye=% +2) 43) 4d — Fooly 48) BH2 4 erect 1) 8, 2. (28) 


The problem will be treated by 


where €= x/l, and / is the overhang. This set of functions is suitable for the treatment of the 
oscillations of any cantilever beam which does not taper to zero thickness at the tip (see the 
Appendix). When the flexural displacement is assumed to be given by 


y= 19, Y5 5-05 eis Lan ee ry. Ae = 43 9 (2, 2a) 
the equations (2.1.2) or (2.1.3) reduce to (compare §10 of R. & M. 1798) 
16 104 8 1304 16 2159 
als — met) + %(3~ ger %) + 20 ao — topes On 2 ae 
8 1304 80 1304 S022 
a Ede BY les e Sgt Se epee y) 
ai( 5670 “/ + 22\ 9) ~ 6237 “7 da( Be MESS aa 
16259 SOM 27 40. 227 
a 7 ~jog03 *) + (F — Figg) + CG jaye) =O 225) 
2/4 
where a _ ; Ke i. a be ie S: me nae (2228 


m == mass per unit length, 
«) = 27 times frequency, 
EI = flexural rigidity. 
The eliminant of (2.2.3) — (2.2.5) is exactly 


85103 — 4,159,51802 + 1,865,239,200« — 22,424.688,000 = 0, me Bae 
of which the roots are 
n= a2 362307 
C5 == eet OO: 


a = 4389-88. 
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In the following table these results, together with those obtained by the use of Y, and Y, 
and of Y, only, are compared with exact values of «, which are the fourth powers of the roots of 


1 + cosh mcos m = 0 ae ee 2 ae e ee a) 
TABLE 1 
Values of Non-Dimensional Frequency Parameter for the Flexural Oscillations of a Uniform 
Cantilever 
Root («) 
Fundamental First Overtone | Second Overtone 
Value obtained by Rayleigh or age . 12-46 = _- 7 
Galerkin method with displace- Y,and-Y, 12-3625 515-86 | ~- 
ment functions as stated. Nay io Alicy Y g 12- 362367 485-56 | 4389-88 
Exact.value 2. As ; 2a a 12-3623* 485 -522 | 3806-55 


* Last digit doubtful. The root of (2.2.8) is only known to six figures. 


It will be seen that one function yields a fairly correct value for the fundamental, two functions 
give a very accurate fundamental and a rather inaccurate first overtone, while three functions 
lead to very accurate values for the fundamental and first overtone, but a very inaccurate second 
overtone. 

Example 2. Torsional Oscillations of a Parabolic Cantilever.—The sections of the beam are all 
geometrically similar, and their linear dimensions are proportional to the square root of the 
distance from the tip. Hence both the torsional stiffness of unit length C and polar moment 
of inertia J are proportional to the square of the distance from the tip, so that 

Ce Cate ac) ae ts Oa: ¥ 2 " A Gare) 
and Jit =a heh seat ee Be Pee 22a) 
where. C, and J, are the values of C and J aagesaensk at oe athe at e is x/2 as usual. 
Accordingly the differential equation for the torsional displacement @ in any pure mode of 


oscillation is 
d dO 
J+ a 5a gee Se = 99 
“ake é) 7+ o(1 £)?0 = 0, See vipa ae a?" (2.2 “Tas 


2 2 
where. ee eo es -. et P. fog A acon. Le) 
Co 


The exact solution has been given by P. B. Walker® in terms of Bessel functions of order +3, 
but this can be more conveniently expressed as 


p= sen va PE te ee elt 2 13) 
The values of the non-dimensional frequency parameter « are 7, 47?, On”, . 2 2. n?r? . 


In the approximate treatment of the problem by the Rayleigh or Galerkin methods it will 
be assumed that 


0 = q,9, +- Jo, + 930s oe o. oe oe e Vic (2.2.14) 
wheref @O,= (r+ Ie — re. he se os i oe (eel 
Equations (2.1.5) or (2.1.6) become when ues 

4 8 I ne in 
(5 108") + (3 ie x) aye Cae is Cee edie 


G - 0 x) + qe = ~ 0) + 95 eee ate = Oiachn, Wag 21a) 


a = 29 
n (5 - 38 x) ae Es 530 *) + 3 es 3465 %) = Ot Meda 2. on kS) 
+ See the Appendix. ; 


op) 


The eliminant of these equations is exactly 
203 — 3270" + 12,0600 — 89,100 = 0; ie Ms = om Veen 
of which the roots are 9-870019, 39-5742 and 114-056. A comparison of these with the true 


values of « for the fundamental and first two overtones is given in Table 2, which also contains 
the results obtained by the use of ©, only, and of ©, and ©,. 


TABLE®2 
Values of Non-Dimensional Frequency Parameter for the Torsional Oscillations of a Parabolic 
Cantilever 
Root («) 
Fundamental | First Overtone | Second Overtone 
Value obtained by Rayleigh or 6, 10-5 Be = 
Galerkin method with displace- @, and ®, 9-87076 46-069 a 
ment functions as stated. @,, ®, and ©, 9-870019 39° 5742 114-056 
Exact value .. ..  .. .. .. | 9869604 | = 30-4784 || 88-8264 


Again the rule holds that an accurate value for the m“ root is obtained by the use of (7 + 1) 
functions. For instance, the error in the frequency* of the first overtone (second root) obtained 
with two functions or degrees of freedom is 8 per cent., but the error falls to 0-12 per cent. when 
three functions are used. 


2.3. Conclusions regarding the Direct Use of the Rayleigh and Galerkin Methods.—The examples 
given above and many others show that when it is desired to obtain a closely correct approxima- 
tion to the frequency in the m" mode} it is necessary and according to experience also sufficientt 
to use (7 + 1) distinct functions or degrees of freedom. Now the labour entailed in the 
computations increases very rapidly as the number of degrees of freedom rises, and it becomes 
of importance to devise means for reducing this labour. Two entirely different methods are 
described in §3 and §4. The second of these appears to be particularly effective and has the 
advantage that it involves no generally unfamiliar notation or process of computation. 


§3. Matrix Solution of the Equations obtained by the Rayleigh or Galerkin Methods. || 
—3.1. General.—The Rayleigh and Galerkin methods lead to sets of linear homogeneous equations 
containing a characteristic number (non-dimensional frequency parameter), and a large part of 
the labour entailed in any particular problem is involved in the solution of these equations. 
A process for solving the equations without the laborious expansion of the Lagrangian determinant, 
and making use of the notation for matrices, has already been published*®, and its application 
to the present problem will be briefly described. 


The equations yielded by the direct application of the Galerkin and Rayleigh methods may be 
written 


(5  00y4)9, 4 (Cis = OAqo)Gs 4 ELCs ==10 | 


(2o1 — %Ag1)41 + (C22 — OAo9)qo + etc. = i 
ele; 


(3.1.1) 


This set of scalar equations can be represented by the single matrix equation 
Cf > Cag V); “it AY Sea ae 


* Note that the frequency is proportional to the square root of «. 

+ The fundamental is the first mode, the first overtone is the second mode, and so on. 

{ Hitherto the methods have not been applied to modes higher than the third. 

| Readers who are unacquainted with matrices may omit this section without detriment to their understanding of 
the rest of the paper. ; 


yi 


where c and a are the square matrices formed from the elements cj and ajj respectively, and ¢ is 
the column matrix or modal column formed with the elements g. When (3.1.2) is pre-multiplied 
by clit becomes 
g — ac ‘aq = 0, 

or Dg og; = ss a ms 7 *: i uc, repales) 
where Dec *a ‘Br BE ng a oe a mz meme Cea 9 
is the dynamical matrix. Equation (3.1.3) shows that «71 is a latent root of the dynamical matrix, 
and q is the associated modal column. 


The matrix method of solution of (3.1.3) consists in taking an arbitrary column 2, and repeatedly 
multiplying this by D to obtain in succession columns 2), 2, ... %. After a time the columns 
so obtained become proportional to one another and to the modal column for the fundamental 
root, and the fundamental root is then given by the ratio of an element of z, to the corresponding 
element of 2,,,;. The higher or overtone roots are obtained in succession by using modified 
dynamical matrices possessing all the latent roots of D except the fundamental and other roots 
already found. This process will be illustrated in an example. It is fully explained in §§4.18 
and 10.2 of Reference 5, and in Reference 4. 


3.2. Example of the Matrix Solution.—The matrix method will be applied to the problem of 
the torsional oscillation of a parabolic cantilever, of which the direct solution is given in Example 2 
of §2.2, The equations to be solved are (2.2.16)—(2.2.18), and these when multiplied 
throughout* by 13,860 to clear of fractions become 


g, (11088 — 1056) + 9, (5544 — 627) + gy (3168 — 396x) = 0, 
gy (5544 — 6272) + 9g, (4752 — 418a) + 9g, (3564 — 286x) = 0, 
g, (3168 — 396z) + g, (3564 — 286%) + gy (3168 — 208) = 0. 


Hence == le LOSGS O27) “S96 
627 418 286 / 
396 286 208 


and c=/[ 11088 5544 3168 
9944 4752 3564 
3168 3564 3168 


The reciprocal of the lastft is found to be 


tis ats eae 
2640 990 1320. 
EY he a ey 
990 495 1980 
ice yay 
1320 1980 1980 
my pace de ees: ath 
Hence D=cla= 15 790 1980 
1 9 2 
15 45 99 
a9 Prt ae 4 
60 120 165 | 


* The same multiplier must be used for all the equations. It can be shown that when different multipliers are 
used the dynamical matrix is unaltered, but errors are liable to arise in the calculation of the overtones since the 
matrices a and c are unsymmetrical and are not the true inertia and stiffness matrices (see Chap. 10 of Reference 5). 

+ For methods of calculating reciprocal matrices see Chap. 4 of Reference 5. 
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The iterative solution was begun with the quite arbitrarily chosen column having elements 
0-9, 0-95, 1:0. After sixteen multiplications* by D a column is obtained with elements 
proportional to ; 

1-0 
1- 139151138 
—0- 093075413. 


When this is multiplied again by D the result is 


1-0 
0- 101316929 1-139151138 | 
—0-093075413 
Hence for the fundamental 


I 
% = 9.401316929 = 9°870019 


and the modal column is 
fe. () 
1- 139151138 
—(- 093075413. 
The root agrees exactly with that obtained by direct solution (see § 2.2, Example 2). 
When the fundamental is absent the numbers q satisfy the equation 


¥1191 + %oiJ2 + %3i973 = 9, -- vs - ». (3.2.1) 
where the coefficients x can be obtained as a row by repeated postmmtuplica eae of an arbitrarily 
chosen row by D until proportionate rows are obtained, or with less labour as a column from the 


product ag, where g is the modal column for the fundamental already obtained. In either of 
these ways it is found that 


Ee ster aii] — | 50 0-621064 0- 405239 |, 


where for convenience the first element has been made equal to unity by multiplying throughout 


by the appropriate factor. The modified dynamical matrix from which the fundamental has 
been eliminated is 


TD), = gee oe | 1-0 0-621064 0405239 | 


0 -+ 0-003040195 — 0-006813927 
0 + 0-018684395 -+ 0-030996411 


Repeated multiplication of an arbitrary column by this matrix yields the first overtone root 


28 ° — (:009459805 — n-ooes1927 | 


1 
in agreement with the direct mine and the corresponding modal column is 
it) 
1- 426740 
— 4-654280. 


* The modal column repeats to six figures after only eleven multiplications. 


ba] 
If in addition to the fundamental the first overtone is also absent, the numbers g satisfy an 
ge relation 
%1291 + Xo092 + Xgo3 = 9. -- +e. -» (3.2.2) 
The octicients x can be obtained by repeated saciorabivau tienen nee an chi trary, row by D,, 
or as the row ag, where g is now the already obtained modal column for the first overtone. It i is 
found that 


EE ich se | x [ 1-0 Se 00g561) = 1-526384 |, 


When (3.2.1) and (3.2.2) are both satisfied it follows that 

91° Go: 93 = 1: — 3:565825 : 2-997253. 
These then are the ratios appropriate to the second overtone. Let D operate once on the column 
formed from these numbers. The ratio of corresponding elements in the original and derived 
columns is 114-056, which is therefore the second overtone root in accordance with the direct 
solution given in §2.2. 


$4. Method Depending on the Orthogonal Properties of the Modes.—4.1. Nature of the Ortho- 
gonal Properties.—The displacements of any conservative dynamical system performing small 
oscillations in any pure mode possess certain orthogonal properties® ®, which are of the following 
general nature. The sum (or integral) taken over the whole system of the continued products of 
the corresponding displacements in two distinct pure modes and the appropriate inertia coefficient 
is zero, and there is a somewhat similar orthogonal property in which the inertia coefficient is 
replaced by a stiffness coefficient. In the present paper the examples are concerned with the 
flexural and torsional oscillations of beams or blades; accordingly the orthogonal properties will 
be explicitly stated for these cases.* First consider flexural motion, and suppose that yr and ys 
are the amplitudes of the displacement for any two distinct pure modes of oscillation, while 7 is 
the mass per unit length of the beam (in general a function of x, the longitudinal co-ordinate). 
Then 


| myrysax =O, .. = s; 2: a a ‘. oan, (4.858) 
if a? 7 an Ss € 
and | El ax = 0, it ES i, te = ON S= 3 a (4 BIO) 


where the integrals extend over the whole beam. For the case of torsional motion let the 
amplitudes in the distinct modes be 6; and 0s, and let the polar moment of inertia of unit length 
of the beam and the local torsional stiffness of unit length be J and C respectively. Then ihe 
two orthogonal properties are represented by 


[ JO:Odx = 0, prin SS arte Saad (55 BRN Saree tale ender eB Bs) 


(ice ¢ Uo: OA 


and Bees he a =: re ms a ec ae 45.4) 


For convenience the seem property depending on the inertias (e.g. (4.1.1) and (4.1.3)) will 
be called the first orthogonal property, and that depending on the stiffnesses the second 
orthogonal property. 


4.2. Outline of the Method.—The method is applicable to any body or system whose equations 
of motion are linear with constant coefficients, but for definiteness it will be explained in relation 
to the flexural oscillations of a beam. As in §2.1 let Y,, Y,, Ys, etc., be a sequence of functions 
satisfying the end conditions. The first step consists in finding the fundamental mode by the 
ordinary Rayleigh or Galerkin method (see §§2 and 3), on the assumption*t that this mode is-a 
linear combination of Y, and Y,. Let the fundamental mode determined in this way be 


Sieve BPR a Ate on 7 = AS 


* The ends of the beam are re supposed to be free, simply supported or rigidly built-in, 

+ Any number of degrees of freedom can be used when determining the fundamental and overtones. However, 
experience indicates that two degrees of freedom are sufficient unless very high accuracy is required. Throughout 
the explanation two degrees of freedom will be assumed, but the modifications in the method required for a greater 
number are obvious. 
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Then any possible mode of seeieeren y from which the fundamental is entirely absent can 
be expressed linearly in terms of the first, second and higher overtone modes, and Saat by 
the orthogonal properties of the modes— 


{ mpydi 00° SO Se a a ee 
and BS: =O. ee Ok ee 
Equations (4.2.2.) and (4.2.3) express the first and second orthogonal properties respectively. 
Next, take the functions 

Tea Vat e Nee tae os A “3 us .. (4.2.4) 
and Mo SY pt Neier ee 5 nt: o. (4.2.5) 


and choose the constants g and / so that y, and xy, are both arhoeerd to the fundamental. 
Then the fundamental will be absent from any linear combination of y, and yp», so that if 4, and 
fy are used as displacement functions in the ordinary Rayleigh or Galerkin process, the lower 
root of the resulting frequency equation must correspond to the jirst overtone of the beam. The 
coefficients g and / are very easily found. It is entirely a matter of convenience which orthogonal 
property is used, but if it is, say, the first, then 


[ mr, dx (0); 


and MV ax = 0. 
Hence | mY ,Vy,dx% + g| TIN a VGN A) eee Me abies ma in es a Als) 
and = [ mY,y,dx +h { mY sydk = 0, ta 


from which g and / can be determined, since y,, Y,;, Y. and Ys; are all known. 


The foregoing process yields a good approximation to the first overtone frequency and mode. 
In order to obtain the second overtone frequency and mode the procedure is as follows. 


Take the functions 
tee Ve ae ee x - ne .. (4.2.8) 
and Has Ng OY, aoe - >) 5 JA. 2o} 
and choose the constants 7, k, / and m so that y, and y, are orthogonal to both the fundamental 
and first overtone modes as already found. These conditions lead to a pair of simultaneous 
linear equations for j and k, and to a similar set in / and m, which are readily solved. Finally x, 
and y, are used as displacement functions in the ordinary Rayleigh or Galerkin method, and 
the second overtone frequency is obtained as the lower root of the resulting frequency equation. 
Evidently the method indicated above can be extended to find any desired overtone. At 
each stage only two displacement functions are used, and both of these are made orthogonal 
to all the lower modes already found. However, it must be remembered that the fundamental 
and lower overtone modes are only determined approximately, and there must be a tendency 
for the errors to accumulate as the order of the overtone rises. Hence it may be anticipated 
that it would be necessary to determine the fundamental and lower overtones with special accuracy 
if it were desired to investigate an overtone of high order. This might require the use of more 
than two degrees of freedom at each stage of the process, but even so the labour would be much 
less than that entailed in a direct attack on the problem. 
‘ A number of illustrative examples are given in §4.3. 


4.3. Examples of the Method Depending on the Orthogonal Properties of the Modes.—Example 1 
Flexural Oscillations of a Uniform Cantilever.—The first step is to obtain an approximation to 
the fundamental mode, which is assumed to ie a pian combination of 


Y, =22 Fe 45%, 


and ania so 
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The Rayleigh and Galerkin methods lead to the equations (compare Example 1 of §2.2) 
16 «104 8 1304 
Ce cn anit 


405 5670 
8 1304 20 1304 
at is aaa Jt 4s gi — aoag® ) Os 


The result of eliminating g, and gq, is a Bete equation in « of which the smaller root is 
12-3625. When this value of «is substituted in either of the foregoing equations it is found that 


2 — 0144100, 
11 
so that the fundamental mode is approximately 


y, = Y, +0-144100Y,,. if ) yards (43.0) 


The next step is to find the coefficients g and which appear in the equations (4.2.4) and (4.2.5) 
defining the functions y, and y, which are orthogonal to the fundamental mode. If the /irst 
orthogonal property is used, then, since m is here a constant, (4.2.6) yields 


[Yd 
Soe = — |-114644, 
[Yond 
while from (4.2.7) 
| Yooax 
SS ee 105769. 
[Yandx 
Hence 
= Y, — 1:114644Y, 
= 26% — §-048813&3 + 4-04881324 — 1-114644&5 bo sae (45382) 
and 
- Ne = 3°333333E3 — 8-86217854 + 7-634614—° — 2-211538e wae (4 Seal 


When 7, and xy, are adopted as displacement functions the Rayleigh and Galerkin methods lead 
to the following quadratic equation for the frequency parameter 


a” — 4861-1840 + 2,124,577 = 0. 

The roots of this are 485-55 and 4,375-6, and the mode corresponding to the smaller root is 
| Yo = m + 0-661018y,. eg * a a Be (4.3.4) 
The approximation 485-55 to the first overtone root nla be compared with the value 485-56 
Ae ie by the use of three degrees of freedom (see §2.2, Example 1), and with the true root 
The results obtained when the second orthogonal property is used instead of the first are :— 

g = — 1-114644*, 
tee a LO ONY; 
a = 485-62 (or 4340-1), 
Vo = M1 + 0° 659083 y,. 

Since the differences from the former results are very small it is evidently immaterial which of the 


orthogonal properties is used. In the present case the labour is reduced by using the second 
property, and this will be adopted in working out the second overtone. 


*It will be noted that this is identical with the result obtained from the first “orthogonal property. The two 
properties will always yield identical results when the function made orthogonal to a mode is a linear combination of 
exactly the same displacement functions as those of which the mode itself is compounded ; for instance, here both 
the fundamental mode and », are linear combinations of Y, and Y, only. Another illustration occurs in Example 4 
below. 
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In order that 7; shall be orthogonal to the fundamental mode and Ist overtone it is necessary 
that the following equations be satisfied 
7 + 0-902870k + 1-114643 = 0, 
7+ 1-746696k — 0- 144100 = 0. 


Hence q = — 2:461462, 
and k = .1-491709, 
Similarly 1 = — 2-440966, 
and m = 1:466118, 


when 7, is orthogonal to the fundamental mode and Ist overtone. When yz and yy are used 
as displacement functions, the Rayleigh and Galerkin methods lead to a quadratic equation for 
x whose roots are 3,802: 15 and 82,393. The true value of « for the second overtone is 3,808-55, 
to which the lower root of the quadratic is a close approximation. 


Example 2. Flexural Oscillations of a Parabolic Cantilever.—The sections of the beam are 
geometrically similar and the linear dimensions are proportional to the square root of the distance 
measured from the tip. Let 

x — distance measured from root, 
/ = overhang, 
2 = x/L, 
m = mass per unit length, 
EI = flexural rigidity of the section, 
and denote values measured at the root by addition of a zero suffix. Then clearly 


m= mM, (1 — &), 
and 18 Fao hea 
The differential equation for the displacement in any pure mode is 
a? a? 
calla Ms pt oa | = (111F\y SO), at. vee lo a oe 
__ moorl4 
where C= a - 5 me ee 3 # oe .. (4.3.6) 
The on conditions are 
dy 
Ae dz = 0 when & = 0, ee eR 
dy é vt 
a =5e( E158 when = 1 4.. 6 = | er 


No exact solution of this problem has been obtained, but it follows from results given on p. 108 
of Watson’s “‘ Bessel Functions’’’ that the solution is associated with a series of a generalised 
hypergeometric type ; in Watson’s notation the series 1s 


2 eee 
Fs (5, La apie ce) ne Uk ae ee aa eam er 
where ee he ts  Eeeak oY Ae = eps i a = .. (4.3.10) 


The series solution has been discussed by Morris®, and his results when converted to the present 
notation give the following values of « for the fundamental and first overtone :— 


39-674, 514-8. 
The parabolic beam tapers to zero at the tip, and it follows from the discussion given in the 
Appendix that the conditions’ (4.3.8) reduce to d*y/d@? = Owhen €=1. Accordingly the 
displacement functions 


‘ver a 1 (r a8 2) Ertl __ dy ert? 1; Sh. ys des sen ee (4.3.11) 
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are suitable for use in the Rayleigh and Galerkin methods. If the displacement is assumed to 
be given by 


y= Y,' + bY’, dey i Re. re iS % ¢ (4.3.12) 
these methods lead to the equations 
9 7 6 7, Q 1 
ee ER oe, SER CC) 
os 


The quadratic eliminant reduces to 
170? — 10,962 + 408,240 = 0, 


of which the roots are 39-684 and 605-14. The smaller root agrees well with Morris’s result for 
the fundamental. On substitution of the smaller root in (4.3.13) or (4.3.14) it is found that 
b = 0-461. 


The displacement functions 


Neat a ols. 
and he NG. Se 
are orthogonal to the fundamental mode already obtained (if the second orthogonal property is 
used) when g = — 1-284 and # = — 1-264. With the orthogonal displacement functions in use 


the Rayleigh and Galerkin methods lead to a quadratic equation for « whose roots are 513-0 and 
1564-3. The smaller of these agrees closely with the value 514-8 obtained for the first overtone 
by Morris. 


Example 3.—Flexural Oscillations of a Conical Cantilever.—The beam is in the form of a cone 
or pyramid tapering to a point at the tip. The sections are geometrically similar, and 


m = my (1 — &?, 
eels (lee )4, 


Hence the differential equation for the displacement is 


za [( 8) ‘oe | — al — By = aes Manat Acres 15) 


where « is defined in equation (4.3.6). 


This problem was first discussed by Kirchhoff and later by several other writers, including 
Wrinch.? These authors succeeded in obtaining a solution in Bessel functions which it is not 
necessary to quote here.* The numerical results which will be compared with the approximations 
obtained by the present method are taken from the paper by Wrinch. 


It is shown in the Appendix that for a beam of this type neither d?y/d& nor d%y/d&3 vanish at 
the tip. Hence the displacement functions 


Vea are es: a in ue We e et a2 (4:0; 16) 
may be adopted. Accordingly, assume that 
ik a hae ee Bh ie! aie oe :. 8 ee (4:33.17) 


The Rayleigh and Galerkin methods lead to a quadratic equation for « (as defined in equation 
(4.3.6)) whose roots are 76-304 and 634-1, whereas according to Wrinch the value for the 
fundamental should be 76-042. Also the value of 6 corresponding to the smaller root is 1-3526,. 


* On this point the reader may also consult Watson’s ‘“‘ Bessel Functions ”’, p. 408, or Maclachlan’s “ Bessel Functions 
for Engineers’, p. 114. 
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It is found that the functions & — 1-5525842° and @ — 1-498569&4 are orthogonal (second 
property) to the fundamental mode already obtained. When these are used as displacement 
functions a quadratic equation for « results, of which the roots are 460-4 and 2656-6. The 
correct value of « for the first overtone is 447-16. Since « is proportional to the square of the 
frequency, the error in the first overtone frequency is about 1-4 per cent. The comparatively 
low accuracy in this example may be attributed to the fact that the displacement functions 
Y,” do not have to satisfy any condition at the tip, and are therefore less tightly conditioned 
than the displacement functions for Examples 1 and 2. 


Example 4.—Torsional Oscillations of a Uniform Cantilever.—The differential equation for the 
displacement in a pure mode is 


oa t 00 =0, 1 ERE a SRC CS 


where « is defined in equation (2.2.12). The end conditions are 
6 =O when = 0; 


oe = 0 when oe ie 
Hence (see the Appendix) the displacement functions 
, O, = (ry + IE: — ert Pe iG gotte, ae Aa 


are suitable for the treatment of the problem. When the functions @, and @, are used, the 
Rayleigh and Galerkin methods yield a quadratic equation for « whose roots are 2-4680 and 
23-5625, whereas the correct values for the fundamental and first overtone are 


a Eetoe Or* LAO <9 

q 2:4674 and m 22° 207. 
The corresponding approximation to the fundamental mode is 9 = 0, + 0-2452790,. Displace- 
ment functions orthogonal to this are 


w= ©, — 1219632 O31" or 2°? property), 
and . No = YO, — 1-189888 0, (1* property), 
or Ne ae 0, Sa ie 1827590, (2™' property). 


When 7, and 7, corresponding to the first orthogonal property are used as displacement fuitiétions 
the results obtained are 


a = 22-402 or 69-661 
Mode, 0 =3 41210-92557: 


If y,; and 7, correspond to the second orthogonal property the results are 


a = 22-397 or 69-158 
Mode 0 = 7, + 0:9328y7p9. 


The error in the first overtone frequency is about 0-4 per cent. 
In order to obtain the value of « for the second overtone, displacement functions 


ns = O, +70, + FO,, 
and N= ©, + 10, + mO,. 


both made orthogonal to the fundamental and to the first overtone mode (as obtained by use 
of the second orthogonal property) will be adopted. The equation for « obtained by use of 
these functions has the roots 63-317 and 136-185, whereas the correct value for the second over- 
tone is 257?/4 = 61-685. Hence the error in the second overtone frequency is about 1-2 per cent. 
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Example 5. Torsional Oscillations of a Parabolic Cantilever.—The problem has.already been 
treated by the direct method in Example 2 of §2.2, and the application of the matrix method 
to the direct solution has been given in §3.2. The differential equation is (2.2.11), and as shown 
in the Appendix the functions ©, as defined in equation (4.3.19) are applicable. When the 
displacement is assumed to be 


6) = @, + b@,, 
the results given by the direct application of the Rayleigh and Galerkin methods are 
a = 9-87076 or 46-069 


Mode 0 = ©, + 1:030266,. 


The correct fundamental root is m? = 9-869604, so that the error in the fundamental root is 
only 0-005 per cent. To obtain the first overtone functions y, and y, constructed as in 
' Example 4 and made orthogonal to the fundamental mode according to the second property 
are adopted. The results obtained are 


a = 39-568 or 114-196 
Mode 0 = », + 3:032527y, 
= @, + 1-4233230, — 4-628607 ©, | 
= 2& + 3-269969&2 — 21-361074E3 + 13-885821 24. 


Since the correct value of « for the first overtone is 47? — 39-478, the error in the first overtone 
frequency is about 0-12 per cent. Lastly, when the second overtone is obtained by use of the 
displacement functions y,; and y, as in Example 4 the results are 


a = 88-880 or 238-958 
Mode 0 = 22 + 5-5882? — 71-68828 + 128-0772 — 62-684é5. 


The true value of « for the second overtone is 9x? = 88-826, so that the error in the second over- 
tone frequency is only about 0-02 per cent. 


§5. Direct Iterative Solution.—5.1. General.—One of the oldest processes for obtaining the 
natural frequencies of vibration of elastic bodies may be called the direct iterative method.* 
According to this system a mode of oscillation, consistent with the conditions of support, is 
assumed, and from this the distribution of the inertia forces in a presumed sinusoidal motion 
having the same frequency at all points is calculated... The inertia forces are taken as a load 
system for the body, and the corresponding static deflexions are calculated from the known elastic 
properties of the body. These static deflexions are then taken as the assumed mode, and the 
process is repeated. In this way a sequence of modes is obtained, and these converge to a 
definite mode,? which is in fact the fundamental mode for the body. The fundamental frequency 
can then be calculated from the mode in the ordinary way, but it can also be found from the 
displacements in a pair of consecutive derived modes as explained below. 


The process of iteration easily leads to the fundamental mode and frequency, but a serious 
difficulty arises when an attempt is made to derive any overtone mode and frequency in the same 
manner. If a mode closely approximating to a certain overtone mode be assumed, it is found 
that the successive derived modes diverge from the true overtone mode, and if the process be 
sufficiently continued they again converge upon the fundamental. This phenomenon will be 
described as regression to the fundamental, and an illustration of it is given below. It may be 
remarked that systems possessing only a finite number of degrees of freedom behave in exactly 
the same way, as can be proved by the theory of matrices *°. This proof can be extended to a 
continuous elastic body if it be regarded as a dynamical system possessing a very large number of 
degrees of freedom. 


* The earliest reference to this method known to the writers is Vianello, Z.V.D.I., 1898, p. 1436. 
+ The mere amplitude does not converge, but the ratios of displacements converge. 
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Attempts have been made to evade the difficulty in obtaining the overtones by invoking the 
principle, inherent in Rayleigh’s principle itself, that the frequency is stationary for small 
departures of the nodes from their true positions. _Hohenemser! makes a preliminary estimate 
of the nodal positions. by an approximate method, assumes the nodes fixed, and finds the funda- 
mental frequency of the body so constrained. Hansen and Mesmer!! !? regard the nodal positions 
as known, and use a modification of the iterative method. For example, their procedure in 
finding the first overtone is as follows. Take an assumed mode having its node at the chosen 
position, apply the iterative process, and so obtain a derived mode which will in general not have 
the assumed node. This derived mode is combined linearly with the fundamental in such a pro- 
portion that the node is brought to the chosen position, and the frequency is found for this 
modified mode. The second overtone can similarly be found by combining the mode derived by 
iteration with proportions of the fundamental and first overtone so chosen as to bring the two 
nodes to the selected positions. 


A modification of the iterative solution is now proposed which dispenses with all preliminary 
knowledge of the nodal positions. This modification makes use of the orthogonal properties of 
the modes (see §4.1), and is as follows. To find the first overtone assume an initial mode, which, 
in order to shorten the work, should be a guessed approximation to the first overtone mode.* 
The usual process of derivation is applied once or oftener to the assumed mode. The derived 
mode is then combined linearly with such a proportion of the fundamental that the result is 
orthogonal to the fundamental, and the frequency is calculated for the orthogonal mode. The 
whole process may be repeated until the calculated frequency becomes constant. The procedure 
in determining the second overtone is similar, but the derived mode is combined with proportions 
of the fundamental and first overtone, and made orthogonal to both of these. The extension to 
higher overtones is obvious. 


5.2. The Process of Dervivation.—First consider the torsional oscillations of a beam or blade 
of overhang /, and measure the co-ordinate x along the blade with the origin at the root. Let 


Pee Ce) | 
cle) Ria: (5.2.1) 
a ‘ 
(eo) = J | (5. Ze 
Jo 
w?J [2 
os cs : (5.2.3) 
where & = x/l, C(&) is the local torsional stiffness of unit length of the beam, J(&) is the polar 
moment of inertia of unit length, and the zero suffix indicates a value measured at the root. 
Let the assumed initial mode be 0, and the successive derived modes 6,,. . . .0,-;, 9,. Define 
the function f ,_,(&) by the equation 
1 . 
Rate) =| (8) 0.48. m4 a e te + Sp ip (5.2.45 
“4! Pas chy) Jans (2) 9 
Then dz °(E) (5.2.5) 
é d0 
fe a nh 9 
and ae dé dé (5.2.6) 


The last three equations define the process by which each mode is obtained from its predecessor. 
It should be noted that 


$5 (Dis OF oe Bane er ele Pa 


and the ceed DSM cc) 2 Se Ss a a 


* A study of known solutions will enable the calculator to make a good choice of mode. 


ie . 

The approximation to « corresponding to the mode 9n as given by the usual energy quotation is 

1 aQ,, 2 

J c(@(S jae 
0 dé 
Cais aE oy De. ae aa - a ot 7 te A 82) 

[7 ©) 90 dé 

This can be put in a more convenient form. For by (5.2.5) 


fe(@(Se) a =| Se fe (O ae 


=[0f8 | — fo. Se az 
~ But 6, vanishes when & = 0 and f,_, (1) is zero by (5.2.7). Hence 
+1 a0, 2 7h ie x 
[ce 7) a= =i fol) i 
a [70-1 Fe IRE i RL Pe ee ee ap et RES BT) 


by (5.2.8). Therefore (5.2.9) becomes 


1 
| 7 (2) On1 Ond& 
a a ee se a oe ne < eos. Lb) 
[ 7 (2 9.2ae 
In the limit when 0, is proportional to 0,_, this shows that 
oe om a se ae og Ey af a eced (One re 
The truth of (5.2.11) also follows from a general argument which shows that a similar formula, 
with @ representing the displacement and 7(&) the appropriate inertia coefficient, must hold in 
other cases. (For instance, in a flexural oscillation 0 would be identified with the normal 
displacement y and 7(&) with the mass per unit length divided by the value of this quantity at 
the root.) For the numerator of (5.2.9) is twice the elastic energy,* while the numerator of 
(5.2.11) is twice the work done by the load system whose intensity is 7(€) 9,1, and these are 
necessarily equal.t 
The details of the process of derivation in the important case of the flexural oscillations of a 
beam deserve some comment. In the usual notation the inertia force per unit span is myw*, but in 
the derivation of the modes the factor o? is omitted (compare the torsional case already discussed). 
From this loading the bending moment M(x) is calculated in the ordinary way. The deflexion y 
is then always expressible by single integrals. For a cantilever or other beam built in at the end 
x% = 0 the formula is 


x Reel (cee ; = 
y= |e — 2") sey de’ Ree WE ace ey, ND 1) 
while for a beam of span / simply supported at the end x = 0. 
si yk eis NaN Se a eg aauilieag 2 
y= | (x ON ieee | =) iene Sleiibast9\05.20e 


9.3. Examples of the Iterative Solution.—Example 1.—Fundamental Torsional Mode of a Para- 
bolic Cantilever.—The symbols have the same meanings as in Example 2 of §2.2.. Take ©, (see 


* The constant factor C, is extracted. 


_} This argument is indicated by Hohenemser!®, 
(53486) B 
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equation (2.2.15)) as the first approximation 0, to the mode. Since 7(€) = (1 — &)? equation 
(5.2.4) yields 


f(@ =J,0 — &)° @, dé, 
(1 — QE- Bd 
(2 — 15& + 2588 — 15&4 + 3&), 


and by (5.2.5) 


AD, iG) ee 
A c(é) Lea) 
a at 
Hence 05 po z 4é, 


artede eee 
and by (5.2.11) 
1 
| (1 — £)76,8, dé 
%= = Saisie AP 
[ (1 — 862 dé 
The value of the frequency parameter corresponding to the mode 9, is 10-5 (see Table 2 of §2.2), 


and the true value is 9-8696. The process has been carried one stage further, and all of the 
results are summarised as follows :— 


Error im frequency 


Mode Value of « per cent. 
0, 10-5 3 
0, QAeD 0-15 
0, 9-870 0-005 


Example 2. Regression to the Fundamental.—The node for the first torsional overtone of the 
parabolic cantilever is at € = 0:5: 
The function 
0, = 6, — 1-50, 
vanishes when ¢€= 0-5 and yields «, = 44-8 as against the correct value 39-478 for the first 
overtone. Hence 6, is a fair approximation to the first overtone. The values of « for 0, and 
the successively derived modes are :— 


a, = 44-8 
ee ON 4 
en ee 


It will be seen that «is already approaching the value 9-87 appropriate to the fundamental mode. 
Example 3. Furst Torsional Overtone Mode for the Parabolic Cantilever.—In order to illustrate 
the convergence of « to the true value notwithstanding a bad choice of the initial mode this 
will be taken as 
Oe fier Oa Fic ee 


This has its node at the tip and yields «, = 21:6. The mode derived from 0, by use of the 
formulae (5.2.4) to (5.2.6) . 


a= gy (2E + 282 EP 4 Tet — 285) 
The fundamental mode as derived in Example 1 (there called 95) is 
Oat a (B2E 4 30e8 — 9458 Bro4 re aS eomaeey 


2520 ( 
Take 6,’ = (0, + k,9)) and make this orthogonal to Op. 
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Then 


f (1 <7 ©)" Qo (0, + RO)dé = 0, 
[a — 2)°,0a8 


1 
[ 1 — 82,2 dé 
= — 3-743824. 
The value of « obtained from the mode 0,’ by means of (5.2.9) is 40-362 (note that (5.2.11) 
cannot be used since 0,’ is not directly derived from 9,). The error in the frequency at this 
stage 1s 1-1 per cent. 
The process of derivation defined by (5.2.4) to (5.2.6) is applied to 0,’ and y ields a mode Q.. 
ec ben. 


so that i 


Galva Oe ela Ug 
is made orthogonal to 0), and the numerical value of k, is found to be 0:000344. The smallness 
of this shows that 0, is almost free from admixture with the fundamental. Finally the mode 
9,’ yields on application of (5.2.9) the value 39-645 for «. The corresponding error in the 
frequency is 0-2 per cent. 


In this example it probably would have been advantageous to apply the process of derivation 
again to 0, before making the mode orthogonal to the fundamental. It may also be noted that 
it was found convenient to use Simpson’s rule in working out the integrals concerned. 


APPENDIX 
The Choice of Displacement Functions for Cantilever Beams 
Flexural Oscillations.—The differential equation to be satisfied is 


Giver GV ae 1 a 
qatl as = mory ae Bc At ut ae as ae we (1) 
For convenience measure z from the tip, and suppose that 
che 3 
— q 
in c(z)z 2) 
an = A(e)z" 


where 0(0), c(0) and d(0) are neither zero nor infinite. Then the right hand side of (1) is 
proportional to z1 when z is small since y does not vanish at the tip, whereas the left-hand side 
is proportional to z?*™~2._ Hence 


De age a = G 
or Pi Pf Pe cr ee a ie (3) 


Flexure. Case 1. Beam not tapering to Zero at the Tip. Brice b=q=0. Hence (3) gives 
¥ = 2, s0 that near the tip 


ay 2 
qe Ue, 
he ; 
and 3 = Qzd(z) + 22d'(z). 
Therefore, 2 and 2 Z are both zero at the tip. Attherooty = & =(. The following functions 
satisfy the conditions 
Y= 4 +2) r+ ae — Lye payne F Yes i\ay ae ee ee 


where € = x/l and / is the overhang, while x is measured from the root as usual (z = / — x). 
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Flexure. Case 2. Beam with Rounded Tip. Ratio of Depth to Breadth of Section neither 
Zero nor Infinite near IT1p.—Here p = 2 and g=1. Hence by (3) 7 = 1, and 


d® 
5 = Ale), 4 = ala) + d'(@)2. 
It follows that d*y/dx* vanishes at the tip, whereas d*y/dx° is not zero there. With the same 
notation as before a suitable displacement function is 


Ye = dre Ze) — ne (5) 
Flexure. Case 3. Beam with Pointed Tip of F oie Maen Tet: p= 4 Ail ae = 2, eid by 
(3) 9 0) Hence nertherig/azamor oe ee vanish at the tip. oma et 


YN Reais te wv iz oF (6) 
is a suitable displacement Finer | 


Torsional Oscillations.—The differential equation 1s 


ad (a9 
(co) = — Jor®. ee o 8 i oe o 5S 54 a8 ote (7) 
Measure z from the tip and let 
Ce clz)z" 
ae (8) 
oe HZ\z 


where c(0), 7(0) and ¢(0) are neither zero nor infinite. Then the right-hand side of (7) contains 2% 
while the left-hand side contains z?**~'. Hence 
aes page | 

or y=q-—-ftl. .. : if (9) 
Now provided that the ratio of depth to Breedth of the eee remains finite near ihe tip Cand J 
will vary with z in the same manner, so that p= q. Hence by (9) 7 = 1, and d0/dx vanishes 
at the tip, but d?0/dx? does not vanish. It follows that the functions 

O=(7+ He —re ee * As De e : Ga 


are suitable to represent the torsional Fgilhresntaits 


The method used here to investigate the geometrical conditions at the tip of a cantilever was 
given in Appendices II and III of Reference 4. 


DIST (OF VREPERENGISS 


No. Author. Title. 
1 W.*]2 Duncan 22% oi + .. Galerkin’s Method in Mechanics and Differential Equations. 
R. & M. 1798, 1938. 
2 W. J. Duncan .. - - .. The Principles of the Galerkin Method. R. & M. 1848, 1939. 
3 Pais; Walkergerc. sc ye .. Simple Formulae for the Fundamental Natural Frequencies of 
Cantilevers. R. & M. 1831, 1938. 
4 W. J. Duncan and A. R. Collar .. A Method for the Solution of Oscillation Problems by Matrices. 


Phil. Mag., May 1934, p. 865. 
5 R. A. Frazer, W. J. Duncan and Elementary Matrices. Cambridge, 1938. 


A. R. Collar. 
6 HeLambs, e bes sus .. Higher Mechanics. Cambridge. 
7 G. N. Watson .. a bo .. Treatise on Bessel Functions. Cambridge. 
8 J. Morris ul. aa: ™ .. Airscrew Blade Vibration. R. & M. 1835, 1939. 
9 D. Wrinch ae a ip .. On the Lateral Vibrations of Bars of Conical Type. Proc. Roy. Soc., 
A., Vol. 101 (1922), p. 493. 
10 K. Hohenemser a 7 .. Beitrag zur Dynamik des elastischen Stabes mit Anwendung auf 
den Propeller. Z.F.M., 1932, p. 37. 
11 M. Hansen and G. Mesmer .. .. Luftschraubenschwingungen. Z.F.M., 1933, p. 298. 
12 M. Hansen and G. Mesmer .. .. Airscrew Oscillations. Aircraft Engineering, March 1935, p. 65. 


(53486) Wt.4 500 1/41 Hw. G.377/1 


SYSTEM OF AXES 


BB eat | 
Symbol . | | 
Axes pecceueet longitudinal lateral | normal 
poe lige forward starboard downward 
Irection 
Force Symbol Xx x Z 
A ee Symbol L M | N 
eck Designation rolling pitching yawing 
Angle of | 
Rotation yea! ? | y is 
i ET CURT Ba Nps Ma ASA ese ERS 
Vaiout Linear u v w 
y Angular p q r 
Moment of Inertia 


A | B | G 
POMS, UREN a) O67 SIGUE AES alae held Rhee EES le 

Components of linear velocity and force are positive in the positive cirection 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of 2. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 


The symbols for the control angles are :— 
€ aileron angle 
yn elevator angle 
ny tail setting angle 
¢ rudder angle 
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Hi psity 
7 ey (where / is a suitable linear dimension) 
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wa er 


b- .-mass. 


ETE | 96 Pent C= eo V aS 
ML Te "Pitching memene with coefficient Gn = M/deV%cs 
P Nye paw 8 moment, with coefficient C, = N/ZeV*0S 


ee 


Al TRSCREWS 
Wn Rev olutignis per second 
Oe, Diameter | 
Sed VjnD a: 

‘ P Power 


ee ieee ais -T Thrust with coefficient nya = - TJon2D* 
Pa ee aL ee O Torque, with coefficient kg = Q/on?D® 


EEE Wig ae steciency, oy ae tw Pies Teta tee 


24313 LIBRARY 


UNIVERSITY OF ILLINOIS 
. URBANA 
The Design of Wind ‘Tunnel Fans 
By 


AD PRE Collar bee Db. OCs 
of the Aerodynamics Department, N.P.L. 


Reports and Memoranda No. 1889 
10th August, 1940 


Summary.—The present report is one of a series dealing with N.P.L. methods for the design 
of return flow wind tunnels. Reports already issued deal with the design of the Compressed Air 
Tunnel! and Open Jet Tunnels,? with the design of corner cascades? and fan straighteners,* 
and with the improvement of velocity distribution by means of windmills’? and gauzes.® 


An explanation is now given of the method usually adopted in the design of:fans for return 
flow tunnels ; as an example of the method, the design of the fan for the N.P.L. non-turbulent 
tunnel is considered in detail. It is shown that for a wind tunnel a fan with blades of constant 
chord is advisable ; and a formula is given by means of which a rapid estimate can be made of 
the variation in rotational speed of the fan corresponding to a variation in power factor at 
constant power. 


1. Power Factor.—The power factor 1 of a wind tunnel is defined by the equation 

CUA. iad Bere. ft os soe J? i ey a. (1) 
where P is the power input (ft.lb./sec.), 9 the air density (slug/ft.*), A the area of the working 
section (ft.2) and V the wind speed in the working section (ft./sec.). The power factor is sometimes 
defined as the reciprocal of the quantity given by (1) ; but the form adopted here is in the nature 
of a drag coefficient, so that a low power factor corresponds to a low tunnel drag and, therefore, 
implies an efficient tunnel. It is, in fact, possible to obtain an estimate of the power factor and 
of the variation of power factor with wind speed by direct addition of the drags of the various 
internal surfaces of a tunnel at different Reynolds numbers; this has been done, for instance, 
by Wattendorf.” 8 


For a new wind tunnel of conventional type, the power factor can usually be estimated with 
sufficient accuracy by comparison of the design with those of other similar tunnels for which 
the power factors are known. A return flow tunnel with a contraction ratio of about 4 and a 
closed working section for which A and V are respectively of the order of 100 ft.? and 200 ft./sec. 
will normally have a power factor in the neighbourhood of 0-25; power factors as low as 0-2 
have been achieved in some tunnels which were designed to have the minimum possible internal 
drag. The provision of such features as a honeycomb and a settling length may raise the power 
factor to 0-3. If the tunnel has an open working section its power factor will usually lie between 
0-35 and 0-45. 


When the power factor of a new design has been assessed, equation (1) is used to obtain the 
wind speed resulting from a given power input, or the power input required to produce a given 
wind speed. It will be assumed in what follows that the power input P is fixed. 


If the area of the fan disc is Ay and the mean axial speed through the fan is w, then by continuity 

A,u = AV, x es ah 2 uf (2) 

whence “ may be found. It should be remarked that in the evaluation of Ay, account should 

be taken of the area of the fan boss; for reasons given in § 4, the boss diameters of the fans of 
N.P.L. tunnels are usually 0-3 to 0:4 of the fan diameters. 


Since uniformity of velocity in the working section is an essential in a wind tunnel, it is 
obviously advantageous to keep the axial speed across all sections as uniform as possible. In 
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AERODYNAMIC SYMBOLS 


1. GENERAL 
m Mass 
t Time 
V Resultant linear velocity 
Q Resultant angular velocity 
e Density, o relative density 
v Kinematic coefficient of viscosity 
R Reynolds number, R = /V/» (where / is a suitable linear dimension) 


R 


ify SS 


x 


Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. ; 
For air under these ( ep = 0:002378 slug/cu. ft. 
conditions | v = 1:59 x 107‘ sq. ft./sec. 

The slug is taken to be 32-2 Ib.-mass. 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio, A = 0?/S 

Lift, with coefficient C, = L/deV°S 

Drag, with coefficient Cy == D/4eV?2S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/4eV0S 

Pitching moment, with coefficient C,, = M/deV*cS 

Yawing moment, with coefficient C, = N/deV20S 


2. AIRSCREWS 


Revolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient ky = T/9n?D? 
Torque, with coefficient kg = Q/en?D® 
Hainciency, 7°s= PV 1 Ps teehee 
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V4 UNIVERSITY OF ILLINOIS 
{fs . URBANA 
0 IS$° ; 
The Design of Wind ‘Tunnel Fans 
By 


AIR COLLAR bee OGe 
of the Aerodynamics Department, N.P.L. 


Reports and Memoranda No. 1889 
10th August, 1940 


Summary.—The present report is one of a series dealing with N.P.L. methods for the design 
of return flow wind tunnels. Reports already issued deal with the design of the Compressed Air 
Tunnel! and Open Jet Tunnels,? with the design of corner cascades* and fan straighteners,‘ 
and with the improvement of velocity distribution by means of windmills’ and gauzes.® 


An explanation is now given of the method usually adopted in the design of:fans for return 
flow tunnels ; as an example of the method, the design of the fan for the N.P.L. non-turbulent 
tunnel is considered in detail. It is shown that for a wind tunnel a fan with blades of constant 
chord is advisable ; and a formula is given by means of which a rapid estimate can be made of 
the variation in rotational speed of the fan corresponding to a variation in power factor at 
constant power. 


1. Power Factor.—The power factor of a wind tunnel is defined by the equation 

TOA > toa Eeey nf * ot “t ~ a 7. (1) 
where P is the power input (ft.lb./sec.), 9 the air density (slug/ft.3), A the area of the working 
section (ft.2) and V the wind speed in the working section (ft./sec.). The power factor is sometimes 
defined as the reciprocal of the quantity given by (1) ; but the form adopted here is in the nature 
of a drag coefficient, so that a low power factor corresponds to a low tunnel drag and, therefore, 
implies an efficient tunnel. It is, in fact, possible to obtain an estimate of the power factor and 
of the variation of power factor with wind speed by direct addition of the drags of the various 
internal surfaces of a tunnel at different Reynolds numbers; this has been done, for instance, 
by Wattendorf.”® ) 


For a new wind tunnel of conventional type, the power factor can usually be estimated with 
sufficient accuracy by comparison of the design with those of other similar tunnels for which 
the power factors are known. A return flow tunnel with a contraction ratio of about 4 and a 
closed working section for which A and V are respectively of the order of 100 ft.? and 200 ft./sec. 
will normally have a power factor in the neighbourhood of 0:25 ; power factors as low as 0:2 
have been achieved in some tunnels which were designed to have the minimum possible internal 
drag. The provision of such features as a honeycomb and a settling length may raise the power 
factor to 0-3. If the tunnel has an open working section its power factor will usually lie between 
0-35 and 0-45. 


When the power factor of a new design has been assessed, equation (1) is used to obtain the 
wind speed resulting from a given power input, or the power input required to produce a given 
wind speed. It will be assumed in what follows that the power input P is fixed. 


If the area of the fan disc is Ay and the mean axial speed through the fan is w, then by continuity 

Aji AN, a a ae be. u ori (2) 

whence “ may be found. It should be remarked that in the evaluation of Ay, account should 

be taken of the area of the fan boss; for reasons given in § 4, the boss diameters of the fans of 
N.P.L. tunnels are usually 0:3 to 0:4 of the fan diameters. 


Since uniformity of velocity in the working section is an essential in a wind tunnel, it is 


obviously advantageous to keep the axial speed across all sections as uniform as possible. In 
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practice, the axial speed through the fan usually varies from point to point ; but the differences 
are not large, and in the absence of information on the distribution (as obtained for instance 
from experiments on a model of the turinel) it is assumed to be constant in the design of the fan. 


2. Strip Theory Formulae.—Corresponding values of w and P for the tunnel are found as in 
§1. To proceed to the fan design, the appropriate strip theory formulae are required. These 
are given in R. & M. 1293,9 and may be written in the form 


c= 4 oNcw*Gcoto, .. oy) bs oe $A (3) 
oe = } p NewHr cot 6, = & Ee =< At bf. (4) 
yy oat Cy 5 
i Sit) cos b ©) 
PSAP 28% 6 
ins sin © COS @ ’ 6) 
Uu 
tan o = , (7) 
ermeNe 5 75) 
1 Ne 
as Bats 8 
laa Payot 8) 
where T = thrust, 


O'= "torque; 
7¥ = local radius. 
N = number of blades, 
¢==aocal chard: 
C, = local lift coefficient, 
>» = local drag coefficient, 
(2 = angular velocity of fan, 
ad, = rotational inflow factor, 
6 = local blade angle, 
oa = local incidence. 
The exact application of these formulae will be explained later. To make a preliminary fan 
design, the equations are simplified by omission of the drag coefficient (which is usually quite 


small) and of the rotational inflow factor. Moreover, since wind tunnel fans are seldom of high 
pitch, # is small in comparison with 7Q, and it is, therefore, possible to write 


; Uu . 

tang =sind=o=—5, .. he “ - ‘© ne . =» ICL 

with cos 6 = 1. Accordingly (3) and (4) reduce to 

dv moat 2 2 

Wo 2 PNer? OPC, ae “Ne as “ aye pos be te (11) 
a = hp Noru OG, i rf = a ra Le >: ee le 
It will be seen that, to the degree of approximation adopted, the torque power input 2 Q and 
the thrust power output wT are equal. ~- The formula (11) can evidently be written down at 


sight if it is assumed that the velocity of the blade element relative to the wind is the rotational 
velocity yQ. Actually, the rotational inflow factor decreases the velocity below 7 Q, while 
the axial velocity component increases it; the corrections are usually of the same order of 


magnitude, so that (11) gives a sufficiently accurate representation of the thrust grading for 
preliminary design purposes. 
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3. Preliminary Design Considerations. —In order to maintain across the sections of the tunnel 
duct the constancy of total head which corresponds to uniformity of velocity in rectilinear flow, 
it is desirable that the pressure increment through the fan disc shall be constant, i.e. , independent 
of radius. In practice it is sometimes found advantageous to increase the pressure slightly 
towards the outer radu, more particularly for open jet tunnels, since a considerable part of the 
energy loss in the tunnel takes place at the boundaries of the stream ; however, the departure 
from constancy is not great. If the pressure increment is #, then by (11) 


fear skew 
ae Qnr dy 4x © Non Cr 
or 


e O? NevC, = 4rb. $i Bs: a i ay. ane 7b): 


If f is to be independent of radius, it follows from 13) t that crC, must be constant along the blades : 
a condition well-known in tunnel fan design. Also, to the present degree of approximation, 


Resp ys, . a e. 7 x Be “fe .. (14) 
so that (13) becomes 


pA gu 0? NevC, = 4rP:; _ cs +? + < = es ie lS) 


It remains to choose Q and the local values of Ne and C, to satisfy (15), in which the other 
quantities are known. 


Now, for several reasons, it is usual to make c constant along the blades. Thus, if the root 
chord is much smaller than the tip the stresses in the blade root may be too high ; moreover, 
a small root chord implies a lift coefficient which may be excessively high. On the other hand, 
if the root chord is large, blade root interference may appear. There is, however, a more important 
reason. Suppose the power factor achieved in a wind tunnel differs from the figure used in the 
fan design (the power factor will in any case vary slightly with wind speed) ; or alternatively, 
suppose that, by the addition of screens or a honeycomb, the power factor is materially altered 
from that applicable to the empty tunnel. This will be the case, for instance, in the non-turbulent 
tunnel at the N.P.L., in which the introduction of certain artificial turbulence screens may increase 
the power factor from 0-3 to 1-0 or even higher. Then the fan may, and in some cases certainly 
will, have to work under conditions different from those assumed in the design ; and the local 
lift coefficients will be different. If, for two conditions of working, the lift coefficients at radius 7 
are G, and C, + 9d¢€,, and ifin both cases the pressure increment through the fan is constant, it 
follows from (13) that both crC, and crSC, must be constant along the blades. Now for the 
linear portion of the lift curve it is possible to write 


(u + a), vf of 7 nf e mee {1G} 


where — yu is the no lift angle, and the slope of the lift curve is for all practical purposes the 
same at all radii. In view of (9) and (10), (16) becomes 


tt 
C= F(u+0 Pp Atlee Meh want. orsal Gab, (17) 


and hence, for the second working condition, 


ae ; aD SG) 


Ce eae a a ee 


Or 


Since crdC; is to be independent of radius, it follows from (18) that c must be constant along 
the blades. 


(58487) A2 


4 


Equation (15) may now be used to fix the values of Q and Nc. In the first place it is usual 
to choose C, = 0:6 at the tip radius 7 =7,. A lower value of C, would lead to a fan of less 
efficiency ; and it is an empirical fact (noted for instance in R. & M.1569") that a value in excess 
of 0-7 is not usually attainable at the tip. In any case, the value chosen leads to sufficiently 
high values of C,; at the root. Accordingly (13) gives 


ON ee Pe nth 


0-6 7,0A.u" 

Either Q or Nc may now be arbitrarily fixed, subject to the following limitations. If © is chosen 
too low, Ne will become large and blade root interference may result ; moreover, the pitch of 
the fan will become high, while the motor design may present difficulties. On the other hand, 
if Q is chosen to be large, the fan blades will be narrow even when N = 2, and the blade strength 
may be insufficient ; in addition, the tip speed may be so high that noise and compressibility 
effects become serious. On the last account, it is usual to limit the tip speed to be not greater 
than half the speed of sound; and as a result of the other considerations, the permissible range 
of Q is usually not great. When Q has been fixed, Ne is given by (19). The number of blades N 
is to some extent arbitrary ; it should not be too large, or the blades become unduly weak, while 
construction of the fan is complicated. For N.P.L. tunnels, N is usually chosen to be 2 or 4. 


Since c and crC; are constant along the blades, vC; is constant and equal to the value at the 
tip. Again, since 7C; is constant, it follows from (17) that r(u-+ 9) is also constant along the 
blades. 


4. The Preliminary Design.—sSo far, it has been shown how to determine preliminary values 
of Q,N,andcfor the fan. It has also been shown that rC, is constant along the blades ; hence, 
unless the boss diameter is fairly large, the blade roots will stall. For N.P.L. tunnels the fan 
boss diameter is usually chosen to lie between 0-3 and 0-4 of the fan diameter ; efficient working 
of the root sections is thus obtained. Moreover, a boss of this size is generally large enough 
to allow the fan motor to be placed inside the boss fairing ; the power unit is then very compact. 
The support of the boss fairing and motor on the fan straighteners is described in R. & M.1885.4 


The blade sections used for N.P.L. fans are usually sections C, D, E and F of the family of 
airscrews, details of which are given in R. & M.829.1° To suit the constant chord required for 
the present purpose, sections C, D and E have been slightly altered, while section F has relatively 
a much thinner trailing edge. The sections have flat undersurfaces, and are thus easy to construct. 
Details of the sections are given in Table 6 and Fig. 1. The lft and drag coefficients for infinite 
aspect ratio are given in Fig. 2 over a sufficiently wide range of incidence. 


To complete the preliminary design, therefore, it is usual to take four equally spaced radii, 
the innermost just outside the boss and the outermost just inside the tip. The sections C, D, 
I, and F are assigned to these radii in the order of increasing radius. From the local value of C,, 
the oes a of the section is found from Fig. 2; and the blade angle is then given by (9) 
and (10) as 


B= a+ anf See eM 


5. Strip Theory Calculations.—A more accurate estimate of the performance of the preliminary 
design of fan found by the methods of § 4 is now obtained by the use of equations (3) to (9). 
To apply these equations a successive approximation process is used, which in practice converges 
very rapidly. 

The values of wu, Q, N, c and 6 obtained in the preliminary design are treated as fixed. The 
local Reynolds number R at which the blade sections operate is given approximately by 
vQc/v, where vy is the kinematic viscosity ; it is usually sufficiently accurate to take a mean 
value of R for the whole blade. 


As a first approximation, put Cp = 0 and cos ¢ = 1 in equation (6), and use the local value 
of C; found in the preliminary design. The resulting first approximation H = C, is then used 
in (8) to find at each of the reference radii a first approximation to 1/(1 — a,), which in turn is 
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used in (7) to give a first approximation to 6. Equation (9) is next used to find «, and thence, 
from the curves of Fig. 2, values of C; and Cp at the appropriate Reynolds number are obtained. 
The values of 6, C, and Cp are then used in (6) to obtain a second approximation to H, which 
is inserted in (8) to obtain a second approximation to 1/(1 — a,) ; and thence from (7) a second 
approximation to d is found. The process of approximation is continued until $ repeats. In 
practice @ is computed to the nearest tenth of a degree, and usually it repeats to this order of 
accuracy after only two steps in the process of approximation. 


When final values of 6, C, and Cp have been found, equations (5) and (6) are used to obtain 
values of H and G, which are inserted in equations (3) and (4) to give curves of thrust and torque 
grading. These curves are integrated graphically or by Simpson’s rule to obtain the thrust 
and torque* ; and thence the torque power input, thrust power output, and fan efficiency are 
obtained. 


There is generally a difference of anything up to about 6 per cent. between the power input P, 
obtained in this manner and the input P assumed in making the preliminary design ; there is 
therefore (since the axial speed is common to both methods) the same difference in power factor ; 
an error of this magnitude is usually within the limits of possible error in the estimation of the 
power factor. To reduce P, to P, both w and Q must be multiplied by (P/P,)!/3, which usually 
represents a change of only about 2 per cent. and is thus not serious. If, however, the error in 
the power factor is considered to be serious, the blade angles or chord of the fan may be altered 
(the approximate extent required is obvious from the calculations) and the strip theory calcula- 
tions repeated. 


Again, from the thrust grading, the distribution of the pressure increment # with radius is 
found. If this is considered to be unsatisfactory, local alterations in blade angle can be made to 
improve it. The fan design is thus completed. 


As a numerical example, the design of the fan for the N.P.L. non-turbulent tunnel is discussed 
in detail in § 8. 

6. General Remarks.—In equation (20), two approximations are involved: the omission of 
the rotational inflow factor and the approximation tan 9 = 9. It so happens that for a blade of 
normal design the errors thus introduced, which are always of opposite sign, are of nearly equal 
magnitude, so that the value of 6 given by u/v Q is very nearly the same, for the whole length 
of the blade, as that obtained in the strip theory calculations. 


It may be argued that the sections of the family of airscrews are not so efficient as other sections 
developed in recent years. However, the sections are quite good for the ranges of C, employed, 
and give tunnel fan efficiencies of the order of 90 per cent.; since, in addition, the sections 
have such simple shapes, refinement seems hardly worth while. 

A short explanation of the derivation of the curves of Fig. 2 may be given here. The lift and 
drag coefficients of the sections of the family of airscrews, corresponding to infinite aspect ratio, 
are tabulated in R. & M. 8921! and R. & M. 17711"; the first report relates to tests in an atmos- 
pheric tunnel, and the second to C.A.T. tests. The range of Reynolds number covered is from 
about 0-1 10® to 6x 10°. The tests show that the scale effect on the lift coefficients in the 
working range (which for the most part is well below the stall) is not great. Accordingly, for 
design purposes it is sufficient to use the curves of Fig. 2 for values of R from about 0:2 10° 
to 610%. These curves are mean curves obtained from all the data of the two above reports ; 
* but the incidences for section F have been changed by the addition of a constant factor to give 
the theoretical no-lift angle for the altered section used in the present report. This no-lift angle 
was found by means of Glauert’s formula,!? which gives sufficiently good agreement with the 
experimental no-lift angles of sections C, D and E. 

The drag coefficient curves of Fig. 2 are also considered to be of ample accuracy for the purposes 
of tunnel fan design, though no great absolute accuracy can be attached to them. Examination 
of the experimental data for sections C, D and E showed that the drag coefficients corresponding 
to infinite aspect ratio decreased with increasing R in a manner which could be approximately 
represented by the equation 

Cp + 0-004 log, R = constant. 


* For usual spacings of the sections, the formula given in Appendix II may be used to obtain the integrated values 
of thrust and torque. 
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Mean curves of variation of Cp + 0-004 log,) (R x 10-®) against «, for the whole range of R 
covered by the experimental data, were therefore drawn; these are the curves C, D and E of 
Fig. 2. The values of minimum drag coefficient thus obtained are only a little higher than those 
calculated by Squire and Young in R. & M. 183814 for modern aerofoil shapes with the transition 
point at the leading edge; the rate of variation with log,) R in the neighbourhood of R = 10° 
is also in quite good agreement with the results of these writers. The curve shown in Fig. 2 for 
the modified section F used in the present report was estimated from the curves for the other 
sections, the difference being based on the results of R. & M. 1838.14 


The stresses and deflections in the blades due to centrifugal force, thrust, and twisting 
moment are briefly discussed in Appendix I. A representative blade is examined : and since 
the ratio of boss radius to tip radius, and the disposition of the sections along the blades, do not 
vary much in N.P.L. fans, the results for this representative blade can be used to obtain a rough 
estimate of the stresses and deflections to be expected in any other blade. If y is the density of 
the material of the blade and E is Young’s modulus, it is shown that the maximum tensile stress, 
which occurs near the root, is given approximately by 


6, i op = 0:25 yr777O? He 12 e5 or 27 eee - - cyt ea 


and that the extension due to centrifugal force and deflection due to thrust, measured at the tip, 
are respectively 


0, UAL ay Oh ee i: ie ae Sy ie a 5) 
CG, = 26 oO S47 ic) aie a a be 4 iy! i« \A2B) 
A rough estimate of the critical flutter speed of the blades of a new fan can be made on the 


basis of formulae given in R. & M. 1518,!° which discusses the flutter of solid cylindrical blades 
in a uniform airstream. 


7. Variation of Rotational Speed with Power Factor.—lf the power input P is fixed, equations 
(1) and (2) show that AV? and Aw? are constant. The variations in axial speed corresponding to 
small variations in A are thus not great. The speed obtained for a given power in a new wind 
tunnel is for this reason usually quite close to the anticipated figure. The rotational speed of the 


fan is usually also very close to the design figure, though the reason for this is not so immediately 
obvious. 


The calculation of the variation of Q with ’ by the exact strip theory equations is best effected 
indirectly. Any reasonable relative values of w and Q are assumed, and by the methods of § 5, 
the corresponding power input P, is calculated. The values of w and Q are then multiplied by 
(P/P,)1/3 to bring the power input to the fixed value P; corresponding values of Q and 4d are 
thus obtained. The calculations are repeated for other values of w and Q, and a curve of variation 
of Q with 4 can then be drawn. 


For most practical purposes, however, it is sufficiently accurate, besides being simpler and more 


instructive, to use a formula derived from the approximate equations of § 3. Substitution from 
(17) in (15) for C, gives 


> Agu O2 Ne a {r crocs ro ee 


and substitution from (1) and (2) for win terms of P then yields a quadratic for Q of which the 
solution is 


On a ree {U4 Vat Kp Kay H, + M: a : (24) 
where 
1/3 
HS =3(a ae Agr(u + 98) ’ 


ora HAS 
Ky Re a Cat) "et Op. 
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The quantities K, and K, are constants since, as shown above, c and 7(u + 9) are constant along 
the blades. As A increases, the first factor in the expression (24) for Q2 decreases while the second 
increases ; for the practical range of 4 the resulting variation in Q is usually quite small. 


By differentiation of (24), a very simple expression results which can be used directly to obtain 
the small variation $Q corresponding to a small variation 3A. The equation obtained, after a 
little reduction, is 

dQ Q F Or(u + 6) — a, 
dx 3A L 207(u + 6)—u 
Numerical examples of equations (24) and (25) are given in § 8. 


8. The Design of the Fan for the N.P.L. Non-Turbulent Tunnel.—As an illustration of the 
methods described above, details will now be given of the fan design for the N.P.L. non-turbulent 
tunnel. The tunnel has a working section in the shape of a regular 16-sided polygon, the diameter 
of the inscribed circle being 7 ft. Accordingly, 


AeecO-O it 2.23 ©. | a ile Bie Rl 28) 


The tunnel is unusually long ; the working section has a length of 45 ft. There follows a short 
expansion leading to the fan section ; there is then a long expanding cone, the diameter reached 
being 20 ft. A fine mesh honeycomb and a settling chamber about 30 ft. long are provided. 
The contraction ratio is about 8. Since the tunnel has so many unusual features, the estimate 
of power factor is comparatively rough. For the empty tunnel it may be as low as 0:3. When 
artificial turbulence screens are introduced, however, it is anticipated that the power factor may 
rise to 1-0 or even higher ; the fan will therefore have to work under widely varying conditions. 
It must evidently be designed to be unstalled at the highest power factor, and will then be 
further from the stall for the less resistant conditions. 


The fan diameter is 8-5 ft., and the boss diameter 3:0 ft.; a motor of 100 H.P. is installed 
in the boss fairing. Accordingly 


eect mee mee Oe me Ee ah (27) 
P=5-5 x 10 ft. Ib./sec. (28) 


For a power factor of A= 1, therefore, equations (1) and (2) give, with po = 23-8 x 10-4 
slug/ft.?, 


(25) 


Ver ioeGitycec, Ut 2 Me in -. - % Tee AS) 
Ve S51 SEC. a re Be: * 3 , cn eo 

Heiee trom (19); since'y = 4-25 tt., 
NP aeee 10192104 Thier? 4 = ne a oa ne emiod!|) 


It was decided that a rotational speed of about 1,000 r.p.m., giving a value of Q just in excess 
of 100 rad./sec. and a fan tip speed of about 450 ft./sec., would be suitable. Hence the following 
figure was chosen 


Nee out, | es ve Hf: x ee ae ie m (32) 
and therefore 


Sere. LOD* I tad seen aa 7 (33) 


The reference radii were chosen to be 1-75, 2:5, 3-25, and 4-0 ft., and the sections C, D, E 
and F were assigned to these radii, respectively. The lift coefficient at radius 7 1s given by 


ee 55) ee Sess. ee 8 (84) 


From the lift coefficient the local incidence « is found from the curves of Fig. 2; the blade angle 
is then given by (20), in which w and Q have the values given by (30) and (33). The results are 
summarised in Table 1. 
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TABLE 1 
u u 
r (ft.) Section Ge a (deg.) oa (rad.) are) (deg.) 6 (deg.) 
1-75 C 1-458 10-0 0-452 25-9 35-9 
2°5 D 1-020 | 4-8 0-316 18-1 22°9 
3:25 E 0-785 3°2 0-243 13-9 17*1 
4-0 F 0-638 2°6 0-198 i Weed 13-9 


The above is an illustration of the usual method for the preliminary design of a fan. In the 
case of the non-turbulent tunnel, however, the blade angles given by Table 1 were not finally 
adopted, since in view of the very long and narrow working section it was thought that the 
total head losses near the walls would be considerable. Accordingly, instead of keeping the 
pressure increment # constant, it was decided to increase it linearly with radius. For constant 
chord, (13) shows that 7C; must then increase linearly with radius also. The relation adopted 


was 
7¥Gy = 1264-690: 28 7 ae Be ra re is 3, Hees) 


so that at section F, p is roughly 30 per cent. greater than at section C. With the expression 
(35) for vC; and the values of wu, Q, and Ne given by (30), (33), and (32), the integral of the thrust 
power output found from (11) is 98-6 H.P., which is sufficiently accurate for the preliminary 
design. 


The blade angles were, therefore, determined as in Table 2. 


TABLES 
y(t.) | Section | & | a (deg.) o(deg.) | 6 (deg.) 
1:75 Coe 1-215 6-2 | 95-9 32-1 
2-5 D 0-935 3-9 18-1 29-0 
3.95 E 0-785 3.9 13-9 17-1 
4-0 F 0-690 3-1 11:3 14-4 


The preliminary design was then chosen to be a 2-bladed fan, of constant chord 15 in., and 
having the sections and blade angles given by Table 2. 


A computation of the performance of this fan by the exact strip theory equations will now be 
carried out exactly as in §3. A mean Reynolds number for the blades is R = 3 x 10® approxi- 
mately ; accordingly the curves of Fig. 2 represent Cp + 0-0019. Details of the process of 
successive approximation to d are given in Table 3. It will be seen that, except for the innermost 
radius, the second approximation to @ repeats the first ; and only one step further is required 
for the radius y = 1-75 ft. It will be noticed that the final values of o obtained in Table 3 are 
practically identical with those given in Table 2, which are found from the approximate equation 
(10) (see remarks in § 6). 


The final values of o, C, and Cp are used in Table 4 to obtain values of thrust and torque 
grading from equations (3) to Sy The pressure increment through the fan is 


p= - NG een ie Ae ee 3 «gf o8) 


and from the last column of Table 4 it will be seen that p increases steadily towards the outer 
radii. Integration of (3) and (4) using the values of Hy cot 6 and G cot contained in Table 4 
gives the torque power input and thrust power output respectively as 


Q QO = 56,600 ft. Ib./sec. 4 ne os a a aos) 
eal Zen) Elet | 
ul = 51,700 ft. lb./sec. oF ‘ 7 %, A £38) 
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The ratio of (38) to (37) gives the fan efficiency as 


7 = 0-914. = a (39) 
, TABLE 3 
: i Ne a u vib 1 : 
ry Section 60 | a | 5 c.. cn tan d 

1-75 e 32-10-0569 | 0-452 1-215 1-069 0-483 25-8 
2°5 D 22-0 00398 0-316 0-935 1-037 0-328 18-2 
3-25 E 17-1 — 0-0306 0-243 0-785 1-024 0-249 14-0 
4:0 F 14-4 0-0249 0-198 0-690 1017. | 0-201 11-4 

sin ¢ cos ¢ o Cr Gs H i = tan ¢ d 

ee: 

0-435 0-900 6-3 1-218 0-0196 1-398 1-079 0-487 26-0 
0-312 0-950 3:8 0-925 0-0130 1-015 1-040 0-328 18-2 
0-242 0-970 3-1 0-771 0-0108 0-840 1-026 0-249 14-0 
0-198 0-980 3-0 0-682 0-0098 0-746 1-019 0-202 11-4 
1 

sin ¢ cos d 0 C. G H | | tan d d 

l-a, 

0-438 0-899 6-1 1-202 | 0-0192 1-381 | 1-079 0-487 26-0 
ES me a! eae | us ESS per, == <= 
eo : ee ae ee 

TABLE 4 
y | H | G Hrcot? | Goot¢ | has 
| | | 
1758 | 1-381 2-725 ey 5-59 | 3-20 
Se = iieet 12015 2-950 7°73 8:99 | 3-60 
3-25 0-840 3:175 | 10-97 joey co Noe 3-94 
4-() 0-746 | 3-435 14-78 | 17-02 4:26 


The value of V given by (29) was calculated on the assumption that 4 = I and P = 100 H.P. 
Since the power is 102-9 H.P. it follows that A= 1-029. To reduce the power to 100 H.P., 


the speeds must all be divided by V1-029. The resulting value of Q is 
Oe a4" 1 racdeseemty ay. as =e + x See Gh" i 
the rotational speed being 995 r.p.m. 


Calculations have also been made for this fan for conditions giving power factors of 4 = 0:70 
and A—0:33. The fan was considered to give a satisfactory performance throughout the 
probable range of operation ; accordingly, no alteration was made in the preliminary design. 


The values of 2 corresponding to an input of 100 H.P. at the three power factors 1-03, 0-70 
and 0:33 are plotted in Fig. 3 as isolated points. It will be seen that the variation in Q is quite 
small. In Fig. 3 the approximate curve of © against A given by equation (24) is also drawn. 
To determine this curve, values of the slope of the lift curve and of 7 (u-+ 0) were required. In 
view of the altered distribution of #, 7 (u + 9) is not now independent of radius ; the variation 
is, however, slight, and an arithmetical mean 

foe Oe 1213 me x ‘gn ‘. ne i (41) 
was obtained from Table 5, which also gives values of the slope of the lift curve, estimated from 
the incidence range 0 — 5°. The mean value was taken to be 


le. gee 
eee) gens wack Me temas et ee, «(4 
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Equation (24) now becomes 
0 = 34-30" {14+ vipassana, sy Be oe .. (43) 


which is the curve of Fig. 3. The agreement between this curve and the three isolated points 
found from the exact strip theory equations is quite good. 


TABU 
r Section | u (deg.) 6 (deg.) u+6 (rad.) r(u + 6) | dC,/da 
1-75 $ 6:7 32:1 0-677 1-184 5-16 
2°5 D 5:5 22-0 0-480 1-200 5-56 
o:25 | E 4:5 17-1 0:377 1-225 5:73 
4:0 F 3°7 14-4 0-316 1-265 5:96 


For comparison, Fig. 3 gives also the curve of variation of w with A. It is 
Wi 83) lO Es rs a ae 7 Ae .. (44) 
and it will be seen that the variation in uw is much greater than the variation in Q. 


Through each of the three isolated points in Fig. 3 is drawn a short straight line having a slope 
given by (25). It is apparent that (25) gives a very good estimate of the rate of variation of 
Q with A. Fig. 3 shows that, for the extreme conditions of working of the non-turbulent tunnel, 
the variation in rotational speed of the fan for a power input of 100 H.P. is only from about 
975. 1.p.m. to 15129 r.p-m. 


The stresses in the fan may be estimated on the basis of the representative fan discussed in 
Appendix I, for which the stresses and deflections are given by equations (21), (22) and (23). 
For the fan discussed in the present paragraph 


7,/¢ = 3°40, a Be . a 4: me ue »»- (49) 
while the maximum value of Q will be about 120 rad./sec. (see Fig. 3) so that 
Oy =p LO Atesec, e : 7 ¢ - .. (46) 


It was decided that the fan should be made of spruce, for which the aera y and Young’s 
modulus E are approximately 


Sg MAU Suily pee of = es a ae oa .. agg 
Eis 240 Oath. * an * a al .. (48) 
Accordingly, equation (21) gives the maximum tensile stress as 


o, + o, = 6:5 x 104+ 8:9 x 104 
= [54a Os lion te bp ij - a: i .. (49) 


which is about one-fifth of the elastic limit of spruce. 
The extension due to the centrifugal force becomes 
Co ee 
= 0-006 in., .. wif - $F. ot si, 
approximately. Similarly, the tip deflection due to thrust is roughly 


é@= 112 x 10 * it 
=A Fale Seeley, - be od ae - 6 ” dee yt de Ne 
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APPENDIX I 


9. Stresses and Deflections in the Blades of a Representative Fan.—The stresses and strains 
occurring in an airscrew are discussed in some detail in R. & M.1274.!° The main stresses may 
be classified under the following heads :— 


(a) Tensile stress due to centrifugal force. 

(6) Tensile and compressive stress resulting from the bending moment due to thrust. 
(c) Shear stress due to thrust. 

(d) Shear stress resulting from the untwisting effect of the centrifugal pull. 

(e) Shear stress due to the twisting moment of the centrifugal forces. 

(f) Shear stress due to aerodynamical twisting moment. 


If the blades have any “ hang-forward”’ or “ hang-back ”’ there will be a contribution to the 
bending moment from the centrifugal forces ; in the present discussion it will be assumed that 
the blade sections are so disposed that there is no such contribution. 


The forces occurring in the blades of a fan obviously increase toward the root ; the thickness 
of the fan blades discussed in the present paper has been chosen to increase toward the root to 
minimise the resultant stresses. Consider a representative fan blade and regard it as being 
encastré at the root radius 


ee ae | ide ae yi = a ng = pee APA 
where 7, is the tip radius; and write 
et eee ' 


: we Me te, es a, (53) 
where ¢ is the current radius. If in the representative fan, sections C, D, E and F are assigned 
to the points x = 0-475, 0-625, 0-775 and 0-925 respectively, then to a close degree of approxi- 
mation 

es ONIN Meee Ry ee Fe eS ns 


C x 


where ¢ is the maximum thickness of the section. Moreover, the sections are derivable from 
each other by orthogonal projection ; and the area of cross-section and moment of inertia about 
the neutral axis parallel to the chord are respectively 


a = 0:745ci, 
be. 0. 05cr? 
which in view of (54) become 
(ee aye AS) aad AV oes ret fe ve re sp Ky SSO 
Deen” 6.5 1078 clxns*, ei: oe “ f: Be IDS 655) 


The material of the blade will be assumed to be isotropic and of density y. A caution should 
be added here that all densities must be expressed in the proper units; if the unit of force is 
one pound weight and the unit of length one foot, then the density must be expressed in 
slugs/ft.? ; if, however, the inch is chosen as the unit of length, the unit of mass becomes 12 slugs, 
and the density must be expressed as the number of these units per cubic foot. 


The centrifugal force on a slice of the blade of thickness 7,dx is yar,?Q?xdx ; hence, in view 
of (55), the total tension at any radius x is 


1 
nde aati OO 6 C10 Saye 2 Tex) ee an Re +s a ey) 
and the stress, 1.e., the tension per unit area, is therefore 
Gare V7, Sle) ee 2 ~ # Ce we DS) 
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The maximum value of this stress occurs when x = 0:5 and is of magnitude 
Gee) Osa see BL “y i = 7 : ee ee 


The variation in the stress over the inner third of the blade is very slight; for the range 
x = 0-4 to x = 0:6 it can be assumed that the stress is given by (59). 


If e, is the extension of the blade due to the centrifugal forces, then 


LAE ead 
ERS mare 


where E is Young’s modulus for the material of the blade, and o, is given by (58). The integral 
of this equation over the whole of the blade gives approximately 


ep Ol y7 2? Ca i sf 4 ss a. a) 


In the determination of the bending stresses due to the thrust the assumption will be made 
that the blade angles are small and that the thrust is normal to the chord. The stresses are 
then probably not smaller than those occurring in an actual twisted blade. 


Equation (11) gives the approximate value of the element of thrust at x = & on each blade as 


dT = for SOP, de 


and the total bending moment at the point x, due to the elements of thrust external to 4, is 
therefore 


M = $por,80? | EXE — ade + 
x 


Now since 7C; is constant along the blade, and C, = 0:6 at the tip, Cj = 0-6; hence (62) 
reduces to 


M0700 0 cr, tot LS eee ee we 5 x One 


Moreover, the neutral axis of the section is very approximately at 0-57 from the chord; hence 
the maximum stress at the section is 


epee Daya 

which on substitution from (54), (56) and (63) gives 
ie 

oy, = 78-1 pr 202(4) a2 (1 — x2 (2 + 9. sat, fected Ee 
This is itself a maximum for * = 0-525; as for the case of the centrifugal stress, however, its 
value is practically constant over the inner third of the blade, and is of magnitude, approximately 

Op 12 Oren rice -s ie: RC = A et (sis: 

If e, is the total deflection of the blade tip due to bending, then 

1 de, _ Mr, (1 — x) 
75 ax EI : 


Z. 


which on substitution from (56) and (63) gives approximately 


ay BY 
@ = 6 OT (Oe ee ee cS 
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The stress given by (65) is tensile on the flat face of the blade and compressive on the curved 
face. Actually, the maximum compressive stress is probably a little greater than the maximum 
tensile stress, since the neutral axis is, in fact, rather nearer the flat face than is assumed above. 
However, the centrifugal tensile stress has to be added to the bending stress ; since the maximum 
stresses due to both causes occur over the inner third of the blade, it is sufficient in practice to 
say that the maximum stress is tensile and of magnitude 


Ceri OE Org eO® hail SH eR (rule Meets epics axl. acs oerey ar (G7) 
In a similar manner it is readily shown that the shear stress due to the thrust is 
ope a Ne Woy el OL rcs 4 Ay uy % te ae COR) 


and the deflection due to shear is 
3Q2 
ey = 0-45 PTE (12) CP TRE PES Ee eee CD 


where wis the modulus of rigidity. As is usual for a loaded cantilever of slender form, however, 
these shearing actions are small in comparison with the bending actions. 


As regards the torsional actions classified above under the headings (d) and (e), these involve 
the blade angle and rate of twist, and cannot profitably be discussed for the representative 
blade. They can be evaluated by the methods described in R. & M. 1274.16 It may be remarked, 
however, that the twisting moments due to these two causes are of opposite sense: the direct 
centrifugal pull tends to untwist the blade, i.e., to increase the blade angles, while for normal 
blade angles the moment of the centrifugal forces tends to throw the blade sections into the 
plane of rotation, and thus to decrease the blade angles. For fans with relatively low tip speeds, 
such as are discussed in the present paper, the resultant moment is probably not great. 


Finally, the aerodynamical twisting moment is in practice negligibly small. Values of the 
aerodynamical pitching moment coefficient for the sections of the family of airscrews, measured 
about the quarter-chord position, are given in R. & M. 1771. For each section, the moment is 
almost independent of the value of C,; it follows that the lift force acts at the quarter-chord 
point. These pitching moment coefficients are, as might be expected, very nearly proportional 
to the thickness-chord ratio of the section, and for the representative blade, they can in view of 
(54) be taken to be given by — 0-05x-1. Now the centre of twist of the sections, estimated from 
the formulae of R. & M. 1444,1" is at 0-34 c from the leading edge. Hence the total pitching 
moment coefficient, measured about this axis of twist, is for the section at radius 7.x 


C,, = C, (0:34 — 0-25) — 0-05271. 
But xC, has the constant value 0-6 along the blade ; hence 


C., = 0:054x-1 — 0-05x71 : 
= 0:004x—1. ie oy oy apie ie: at * = {10} 


No accuracy can be attached to the numerical coefficient in (70) ; it can only be said that it is 
of this order. It is, however, so small that the resultant torsional moments are negligibly small. 
On the basis of equation (70), a blade of chord c = 1 ft., tip radius 7, = 5 ft., and having a tip 
speed of 500 ft./sec., will be subjected to a total aerodynamical twisting moment, measured at 
_ the root, of only about 2:5 Ib. ft. 


To sum up, it may be said that in practice equations (66) and (67) provide useful criteria for 
the strength of a fan blade of the type discussed in the present paper. If the deflection given 
by (66) is unimportant and the stress given by (67) is well within the elastic limit of the material, 
then in practice the blade is of adequate strength. 
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APPENDIX II 


10. A Formula for the Integration of the Thrust and Torque Grading Curves.—The thrust and 
torque grading are each defined by four points, one at each reference section. In drawing a 
smooth curve through four typical points, and integrating graphically or by Simpson’s rule from 
root to tip, an error of the order of + 0-5 per cent. may be introduced; thus the error in 
efficiency may be + 1 per cent. This error may be reduced by determining, for both thrust and 
torque grading, the cubic curve which passes through the four given points, and by analytical 
integration of this curve. When, as is nearly always the case, the distance from the root to 
section C is the same as the distance from section F to the tip, and when there is equal spacing s 
between the four reference sections, the integral of the cubic curve is given by equation (71), 
which is an extension of Simpson’s three-eighths rule. If the given ordinates are y,, ys, v3 and 
v4, and if the blade length 7, — 7, is /, then 


14 
l P P? ; 
[ var = ge (4-38) 0, +9) +( 27-3) ta +98) f nie - 2 eT) 
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TABET RG 
Details of Atrscrew Blade Sections 
n = fraction of chord from leading edge. 


R, = leading edge radius. 
R, = trailing edge radius. 


Ordinate, as fraction of chord 

: g 

C D E F 
0 0-0212 0-0163 0-0130 0-0108 
0:05 0- 1000 0-0768 0-0613 0-0512 
0-1 0- 1326 0-1018 0-0813 0-0679 
0:2 0-1611 0: 1237 0-0988 0-0825 
0:3 0- 1680 0- 1290 0-1030 0 - 0860 
0-4 0- 1660 0-1275 0-1018 0- 0850 
0:5 0-1591 0-1222 0-0975 0-0814 
0:6 0: 1460 0-1121 0-0895 0) -0747 
0:7 0-1253 0-0962 0-0768 0-0642 
0:8 0-0968 0-0743 0-0593 0-0495 
0-9 0: 0627 0-048] 0-0384 0-0321 
1:0 0:0124 0-0095 0-0076 0-0064 
Rj 0-0212 0:0163 0-0130 0-0108 
R, 0-0124 0-0095 0-0076 0-0064 


No. 
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Symbol : y i 
Axes _ Rete nee longitudinal lateral | normal 
| pe ioe forward starboard | downward 
, rection | 
Force Symbol. Nice Y Z 
M t Symbol L M oN 
rae Designation rolling pitching yawing 
Angle of : 
Rotation Symbol ® 4 Lf 
Vislocit Linear U v ww 
y Angular a gq a 
Moment of Inertia A B ¢ 


or ne eee ee reece en ne 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 


The symbols for the control angles are :— 
7 € aileron angle 
n elevator angle 
ny tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


1. GENERAL 


ZietHr pr >e eM R 


Mass 

Time 

Resultant linear velocity 

Resultant angular velocity 

Density, o relative density 

Kinematic coefficient of viscosity 

Reynolds number, R = /V/» (where /is a suitable linear dimension) 
Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. 
For air under these ¢ e = 0:002378 slug/cu. ft. 

conditions | » = 1:59 x 10-4 sq. ft./sec. 

The slug is taken to be 32-2 lb.-mass. 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio, A = b?/S 

Lift, with coefficient C, = L/4eV?S 

Drag, with coefficient Cp = D/teV2S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/4eV2bS 

Pitching moment, with coefficient C,, = M/eV2cS 

Yawing moment, with coefficient C, = N/ieV70S 


2. AIRSCREWS 
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n 


= to eS ae 


Revolutions per second 

Diameter 

V/nD- 

Power 

Thrust, with coefficient ky = T/on?D4 
Torque, with coefficient kg = Q/on?D5 
Efficiency, 7 = TV/P = Jky/22kq 
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Flow of an Ideal Fluid past a Cascade of Blades 
By 
W. Mercuant, M.A., S.M., 
| of Metropolitan-Vickers Electrical Co., Ltd. 
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87h July, 1940 


SUMMARY .—This is a mathematical investigation into the flow of the ideal fluid of hydrodynamics past a cascade of 
blades. A complete solution is given for a cascade of straight blades set at any angle to the cascade axis and for any 
pitch chord ratio. Previous work on the subject is reviewed and compared with the results of this investigation. 


Part I of this report contains an investigation of the general problem and a general relationship 
which must hold for all types of bladesis found. Part II gives a complete solution for any cascade 
of straight blades at any angle to the cascade axis and any pitch chord ratio. There are also 
three appendices. The first contains an alternative proof of the general relationship due to 
Dr. D. M. Smith. This is based on somewhat simpler principles than that given in the main 
body of the report and is only omitted there because it is not in the main line of the reasoning. 
The second deals with the relation between the lift coefficient for an aerofoil in a cascade and that 
of an individual aerofoil. The last reviews previous work on the subject. For the symbols 
used in relation to a cascade, see Fig. 4. 

Part I. General Theory.—The flow of an ideal fluid is governed by Laplace’s equation 

oN. atv. 
Bt ee ce Cee 
Neri aye 0. 
Both the velocity potential ® and the stream function ‘V obey this equation. 


Consider any function of the complex variable 
w = fle) = E+ in, 
where z = x + ty. 
Then 
VE= V7 = 0. 
Hence we can identify & and y with ® and V’, and from any function of the complex variable 
can obtain a solution of some flow problem. 


Let 
w= O+7P. 
Then 
we u —- 10, 

where 
pe tae at 
Tae ay? 

and 
ae oV 
NF ay ~ Ox’ 


and hence wu and v give the components of the velocity of the motion parallel to the axes. 
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It is also useful to remember that the deflections of a uniformly stretched thin elastic sheet, 
under no loading but having applied boundary conditions, obey Laplace’s equation. We thus 
can picture every flow problem as an elastic sheet problem and can equate contours with either 
flow lines or velocity potential lines as required. 


Consider any flow problem of which there is a known solution and imagine that we have the ®, 
’ lines drawn. These intersect orthogonally everywhere and if drawn for equal increments of 
® and then if one element forms a curvilinear square all form curvilinear squares. In such a 
case we may specify the position of a point either by x, y co-ordinates or by ®, Y’ co-ordinates. 
We can also consider a plane in which ® and are ordinary Cartesian co-ordinates. 


Consider any other flow in the x, y plane with a potential function 
w’ = MO’ + iF". 
Then 
W720 = VV 0: 
Now to every small square in the x, y plane corresponds a small square in the ®, Y plane. 
Hence 


02@' 02Q’ 
Qe | 
om? | ayy 
0 2" Q 2’ 
(6) 6) 
0) = A hh —. PTT STC 
aD oY? 


Therefore to every flow in the x, y plane there is a corresponding flow in the ®, V plane. We can 
thus obtain corresponding flows by plotting the flow lines in the one plane on the other using the 
correspondence between the x, vy and the ®, ‘’ co-ordinates as a mapping function. This process 
is known as conformal transformation. | 


The particular method used in this report is as follows. For any cascade of blades there is a 
direction of flow for which there is no deviation of the fluid stream. Let the angle of inclination 
of this direction with the normal to the cascade axis be 9; this direction is called the direction of 
neutral flow, or more shortly the neutral direction. Suppose the flow net is known for this 
direction. It can be replotted in the ®, ‘V plane as a simple square mesh. In this case the blades 
will transform into finite lengths of the VY lines. 


Now consider a series of laminas uniformly distributed along the x axis, as in Fig. 1. 


ig Le ia 


BiG. 


Further consider the flow past them, with no circulation round the laminas, whose direction 
at infinity makes the same angle @ with the normal to the cascade axis. Since there is no 
circulation round the blades there is no change in the direction of the velocity at infinity on both 
sides of the cascade; there is only a local distortion round the blades. This flow can also be 
replotted as a simple square mesh in the ®, VY plane. The blades again transform into finite 
lengths of the V’ lines. ' If these lengths are the same in the two cases and we know other flows 
round the blades on the real axis, we can use the ®, ’ plane to transfer these flows to flows round 
the actual series of blades in the cascade. 


o 


3 


As a partial problem we will therefore consider the oe flow Dee poe blades uniformly 
distributed along the cascade axis. i dpi JOO 


Consider the flow represented by the eNacon 
iy ly 
ee ee THEE V V2sin hz fe V 2cos } 


dz Vy --cosz Vr + camer" 


The first term represents a uniform flow parallel to the axis of x and of velocity U.. It obviously 
fits in with the boundary conditions for flow past blades on the cascade axis. 


(ple? iP). 


The second term is alternately real and imaginary for periods along the cascade axis as shown 
below in Fig. 2. 


\ Real Ima Real imag. Roel 
J SDs og fg 
. Fic. 2 
\ 


wh Now a real term represents flow parallel to the x axis and an imaginary term flow parallel to 
Aa the y axis. Therefore the flow is consistent with flow round blades along the real axis where the 
‘half length is given by + cos x = 0. 


w\ - Through the centre of a gap, e.g. at x = x, we have for the flow due to the second term 
x Uu— w= Ma v2 Sugden O) 
V7 + cos ie + wy) 
Vv 2coshdy 
Vy — cosh 


Now cosh y is always positive and greater than 1. Therefore the expression above is always 
purely imaginary ; and hence there is no flow across the line « = 7. 
Further 
Lt 
gis 
Therefore the velocity a large distance from the cascade is of constant value and parallel to the 


y axis. Also there are no singularities apart from the real axis. Therefore the term represents 
the flow of a uniform stream normal to the cascade and of velocity V at infinity. 


V 
(u — iv) = 7: 


In similar manner the third term can be shown to represent the flow due to equal circulations 
round the blades corresponding to velocities at infinity of --W parallel to the x axis. 


We thus have the solution for the general flow past the cascade of blades. 
For the flow with no circulation, W = 0; if this flow is at an angle © to the normal to the 
cascade we have tan 9 = U/V, and the stagnation points are given by 
V /-V2sin 2 sin sate 


Atos (ieee (p) aa See y, 
as /7 + cos Gy = oe 


whence 


1 — cos x 
tan? 9 = — eh 
rT - COS 24 


or 
COS %, == cos*o — 7 sin*. 
We thus know the point on the back face of the lamina which corresponds to the trailing edge 
of the blades. 


Now Joukowski’s hypothesis requires that for flow in any direction there shall be associated 
that amount of circulation which will keep the velocity at the trailing edge finite. Changing the 
problem to that of blades along the real axis this requires that the stagnation point should remain 
in the same place for other incident directions of flow than 9. 


Consider a direction of flow, free from the circulation imposed by the blades, of 0. Ifa 
circulation is added to this flow, the stagnation points are given by 


VV 2sin $x, Bs WV 2.cos $24 


Vie We = ad 6) amas eee Bee ane bat ae 
: Vr —eCOS 44 a Vr ST eCOS We 
with 
tan 9 = 5 . 


If these are the same points as for ¢ with no circulation, 
COs ¥%, = cos*o — rsin’g, 
or 
+ 
V2 SIN 9% 


tano = — 
‘ Vr + cos x, 


Hence 


P/Q uf 
0 =tan 0— tang + q ee 
r 


+ cos x 
and it therefore follows that 
V ie a: (J ee r) 7 Sy a 
W (tan 0 — tan 9) = J 1 + i) sec’o. 
Now 
tan == oa : 
tan 68 = i y ne 
tan 0 = e 
Therefore 


tanattan@—2tano /, , (1-7)... 
tana — tan B- = f1+t4 vik 


Put 
y= NEE +=) seco — I 


_ | 1+ cos x, 
1a COSEy 


Then 
y tan a + 2 tang 
y+2 
which is the general relationship mentioned in the summary to this report. 


tan 8 = 


We thus see that the general problem of flow past a cascade of blades resolves itself into two 
partial problems. The first is to determine the neutral direction for the cascade. The second to 
determine the length of the transformed laminas along the cascade axis. 


Part II. Complete Solution for Cascade of Straight Blades at any Angle to the Axis and any 
Pitch Chord Ratio.—Here is known being the angle the blades make with the normal to the 
cascade axis. 


We can determine 7 as follows—When we replot the flow in the 0, Y plane the length 
of the laminas is equal to the change in ® 
$ between the stagnation points. The width of 
the gap is the change in Y’ between each lamina. 
From Fig. 3 it is now evident that the change 
in Mis 
cos ‘{—7) — cos— eu 


udx + 


Gs 


Ag 
where x, is a stagnation at Also, 


4/2, sin 4% 


Uy eV tan e+ ere, 


The integrals can be readily evaluated by 
means of the substitution 


COS 3% = Ni 1 = 1 sec Ly 
Fic. 3 9 
anaes 


or, in terms of the parameter y of Part I, 


1 2 
Gea tan o tan? G aa tne) 4 r+? |. 


21 
$s COS © = [ V dx = nV. 


0 


the result being 


Also 


Therefore 


5 SC 9 


tan © tan— 1 (fate ice ( 2? y 


which completes the solution for the straight blades. A graph of y against s/c for 9 = 0, m/4 is 
appended (Fig. 5). 


In conclusion the writer would like to express his thanks to Messrs. Metropolitan Vickers 
Electrical Co., Ltd., for permission to publish this paper. 
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APPENDIX I 


Alternative Proof of General Relationship 
By Dr. D. M. Smith 


Consider flow of a perfect fluid, of uniform velocity at y = + o and at y = — o, pasta 
cascade of uniform blades arranged along the x axis. The most general possible type of uniform 
flow is compounded of the following velocities at y = + o: 


(1) — V in y direction. 

(2) U in x direction. | 

(3) + W at y = + o, and — W at y = — o, in x direction, due to a certain circulation 
round each blade. 

U+W. 


The resultant approaching velocity makes angle « with y direction, where tan « = V 


acca 
San 


Now consider the local velocity at any point P on the boundary of one blade. This local 
velocity may be written 


similarly for leaving angle 8, tan 8 = 


c= Uu+ Vv+ Wz, 


where uw, v, w are constant for any given point on a given blade at a given setting ; and the 
condition that P should be a stagnation point is that the values of U, V and W should be such 
that 
Uu + Vv + Ww = 0. 
The direction of the flow without circulation which makes P a stagnation point is given by 


tan 9 = —<. 


The condition Vv + Uu + Ww = 0 is equivalent to 


Dee Wig CUR 
Us aN Vu 


ACh ee ee 


—tang+5(1+%)tana+;(1—2)tan p=0. 


== (0) 


Or 


or 


These relations still hold good if P is a sharp trailing edge, because although wu, v, w become 
infinite their ratios remain finite ; hence the last relation is true if P is the sharp trailing edge of a 
cascade of blades the circulation round which is in accordance with Joukowski’s hypothesis. 

The last relation is equivalent to 
y tan «+ 2 tan 9, 

y+2 
wo 
where —y= = 
ye 


tan B = 


/\ 


APPENDIX II 
Effect on Lift Coefficient of Interference due to being in a Cascade 


Specify everything in terms of V, and 8; we have 


La V5. Sit. 0, 
Verge Vi COS.U, 
_ W+V, sind 
Le VE COSeUa ans 
Wee V og Sim 8 
san i Vee COs.Uenia 
ytan « + 2 tan Q 
tan 6 = - y 2 
Hence 
Ww — V. cos 9 (tan 0 — tan 9). 
(yf 1) 
Let the circulation per blade be [. Then 
Lee Ws 
_ 2s V. sin (9 — 9). 
cose (y+ 1)’ 


therefore I‘ varies as the velocity normal to the blades. The lift per blade, per unit length, is given 


by 
li2 — 4Cic elie: 


Lap 
20 4s_ sin (0 — 9) 


So Va. 46 cosoly+ 1)’ 
where for straight blades 
* sec 
wf tan © y + 2V 
21 — 
tan ¢ tan G = + x log a 
To find the value for an individual blade we have to take the limit as——> ©,i.e.asy—>. Now 
Age eres 
Lt y+2 Moe oes 
yo F) )log (5 =) = a —] Bie 
fate) 
Also | 
1 — tan 9 
tale eee 
Ce eae 
Lt re tol On by Me ak 
yon (vy ~|- 1) tan aoe, lat i 1 — tano. 
bs 


Therefore for a single straight aerofoil 
om _4sin (9 — 9) $mseco 
“cos o (tan? o + 1) 
which agrees with previous published results. 


= 2rsin (0 — 9), 
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Write C,. for the lift coefficient of the aerofoil in the cascade. Then 
Cic 2s 1 . 
C, me (y+ 1) cosa, 


where 
we 
ras 9 Sec 2 
aS __ tang 1 Veo can 
tan » tan 16 Th) + 510g ( y ) 
Graphs are given of Coc TT == 41) A (Fig. 6). 
L = 


APPENDIX III 


Brief Review of Previous Theory and Bibliography 
The problem has been under consideration for about thirty years. 


Kutta in 1911 obtained a solution for a set of straight blades normal to the cascade axis. 


Grammel in 1917 obtained a solution for a set of approximately straight blades at any angle 
to the axis. For o = 0, 7/2 the blades are straight and for o = 0 agree with Kutta. 


Glauert and McKinnon Wood, 1918-20, realised the interference due to the blades being in a 
series should be applied to airscrew theory and did apparently the first tests on blades in a cascade 
to obtain modified lft coefficients for use in airscrew theory. 


Konig, 1922, solved the problem in mathematical terms in a manner which would do for any 
shape of blade. The mathematics are very intractable and no subsequent use seems to have 
been made of his method. He gives figures for straight blades at 7/4 to the axis which agree 
with those of this paper. 


Numacht, 1929, calculated out the Grammel Formula as being the best available means of 
predicting the interference of blades in a series. His values for g = 0, z/2 agree with those of 
this paper. They do not agree for» = 7/4, but are of the same order. Thisis due to the Grammel 
blades not being straight for 9 = 7/4. 


Betz, 1931, gave a method of successive approximation to pass from an isolated aerofoil with 
all its properties known to an equivalent aerofoil in a cascade. The method is intended to be 
used for thin slightly curved aerofoils and needs considerable labour. 


Pistolest, 1937, gave a valid method for blades which are a small variation of straight lines 
normal to the axis. The method is an extension of that of Birnbaum-Glauert for single aerofoils. 
He applied the same method to small variations from oblique lines, but the mathematics become 
very complicated. He obtained Grammel’s formula, but for thick blades of an undetermined 
thickness. This made him think that the curved blades Grammel obtains should be reflected 
about their axes. 


Braun, 1932, gave a formula for the effect of varying s/c and « and backs the formula from 
Christiani’s experimental work. He stated that the theory would be published elsewhere. The 
formula does not agree in form with the general relationship given in this report. 
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SYSTEM OF AXES 


| Symbol : ea ee AES, | . 
Axes | yieacbtien | icneiidinal | Mae al normal 
| di erie | forward starboard downward 
| irection 
Force | Symbol | eae ae Z 
Merenee Symbol N 
Designation | rolling Hane yawing 
Rotation ae $ : ’ 
Velecit Linear u v w 
y | Angular b q ¥ 
Moment of Inertia | A B C 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 


The symbols for the control angles are :— 
€ aileron angle © 
y elevator angle 
ny tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


1. GENERAL 


jA~ eps s 


brreo FMR 


tne 


Mass 

Time 

Resultant linear velocity 

Resultant angular velocity 

Density, o relative density 

Kinematic coefficient of viscosity 

Reynolds number, R = /V/» (where / is a suitable linear dimension) 
Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. 
For air under these ( 9 = 0-002378 slug/cu. ft. 

ae: == ROORXLG7 oa bP SEG, 

The slug is taken to be 32-2 lb.-mass. 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio, A = 0?/S 

Lift, with coefficient C, = L/}eV?S 

Drag, with coefficient Cp = D/4eV7S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/4eV2bS 

Pitching moment, with coefficient C,, = M/toV?cS 

Yawing moment, with coefficient C, = N/4eV2bS 


2. AIRSCREWS 


SOs woos 


Revolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient 2, = T/en?D4 
Torque, with coefficient kg = Q/on®?D> 
Efficiency, 7 = TV/P = Jky/2zkq 
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_ Physical Properties of the Standard Atmosphere 
By 
Raw PANKHURST, Achk.G.0., B.SC. 
and 
mete CONN, .UoGr VLA, 
of the Aerodynamics Department, N.P.L. 


Reports and Memoranda No. 1891 
gth “Fanuary, 1941 


Summary.—The accompanying tables and figures give those properties of the Standard 
Atmosphere which are required for the prediction and analysis of aircraft performance. 


Definition of Standard Atmosphere.—In order to compare observations obtained in flight tests 
it is necessary to convert them to standard atmospheric conditions. In this panneenon an 
International Standard Atmosphere has been defined in the Official Bulletin of the International 
Commission on Air Navigation! and reproduced in Air Publication 1275?. The relevant portions 
of the definition may be summarised as follows :— 


a) Temperature.—The temperature is taken as 15°C. at sea level, with a constant rate 
of fall with height of 6-5° per km. up to 11 km. (the lower limit of the stratosphere) beyond 
which it remains constant at — 56-5°C. 


(b) Pressure.—The pressure at sea level is that corresponding to a barometric height of 
760 mm. of mercury at 0° C. with g = 980-62 cm. per sec. per sec. The variation of pressure 
with height is given by the following formulae, which may be deduced from the gas law 


b = eRT together with the relation dp/dz = — gg. 
be ( 288 — 6-5 2 \5250 
foal -) te Foe a0 heme nip 


z—11 
and logy ot = 6 


fOtus ero Lisle. be > on 4) (2) 


z being the height in km. 


(c) Density.—The density is 1-2257 kg. per cu. m. at sea level and varies according to the 
formulae 


988 — 4.256 
aa ae ae 7 ie ae ee aa 


and logy = 


Loree” LP et. He ts Ne a (4) 


Conversion to British s (column 3 of Table 1) 
are very close to round numbers, it was agreed between the Royal Aircraft Establishment and the 
National Physical Laboratory (May, 1925) that the rounded values of column 4 should be used in 
evaluating the temperature and pressure as functions of height; the corresponding values 
agreed upon for pressures in inches of mercury were used as the basis of the calculations of the 
data of Table 2, 


&) 
ye 


A reat 


Rate of fall of temperature - 6-5° per km. 1-9812° per 1,000 ft. | 1-98° per 1,000 ft. 


;| —56-5°C. at 11 km. | —56-458°C. at 36,089-3 ft. | —56-5°C. at 36,090 ft. 


Stratosphere (Eqns. (2) and (4)) .. | 14,6 nn 47,900°3 ft. 47,900 ft. 


For the present work the further conversions involved were 
980-62 cm. per Séc:- persseGyy == 2 ino 20 ee eee eseue ccs 
760.mm. of mercury; (=')5))=)2116<2 sib per saat 
]- 2257 kg..per cu. mi. { ==/p5) ==.0-0023/8earslus pemcusat, 


These were calculated using the following relations* between the fundamental units : 
Il kg. = 2: 204622 Ib. 
Wiaihy —= pee ePAU Chri. 


The density of mercury was taken as 13-5951 gm. per cu. cm. at 0°C. 


The velocity of sound at sea level was taken as the theoretical value 
J rile — 1,117 ft. per sec., 


where y, the ratio of the specific heats, was taken as 1-403 (International Critical Tables)?. 


Viscosity.—The coefficient of viscosity (u) or kinematic viscosity (v) is required for calculating 
Reynolds number. For this purpose the variation of » with temperature (T degrees Absolute) was 
assumed to follow Sutherland’s formula 


T3 
pe ary 
the values used for A and B being those given by Goldstein’, namely 


A = 3-059 x 10-% slug per ft. sec. 
B = 114 degrees, 


The experimental data used in deriving these values were obtained from the International 
Critical Tables® ; although the calculated values are quoted to four significant figures there may 
be some inaccuracy in the absolute value of the last figure, but it is considered that this 
uncertainty may be no more than the probable error of the available experimental data and in 
any case will have little practical effect on Reynolds number. 


Acknowledgement.—The authors thank Mr. F. A. Gould of the Metrology Department, N.P.L. 
for helpful discussion of various points involved in the work. 


REFERENCES 
1. International Commission on Air Navigation. Official Bulletin No. 7, p. 34, Dec. 1924. 
2. Instrument Manual, Section IV, Chapter 1, Appendix II. Air Publication 1275, 1, 1940. 
3. International Critical Tables 5, 81, 1929. 
4. S. Goldstein. Private Communication. 
5. International Critical Tables 5, 2, 1929. 


* The above relations received legal sanction by Order in Council of May, 1898 (Statutory Rules and Orders, 1898, 
No. 411). 


3 


TABLE 2 
Properties of the Standard Atmosphere 
List of symbols : h = height ¢ = temperature a = velocity of sound 
. p = pressure 0 = density o = density relative to sea level 
u = coefficient of viscosity » = coefficient of kinematic viscosity 


The suffix ( )9 refers to sea level (h = 0) 


| a p | ux 107 | » x 104 
£ aC es Tey eee ws Lp Aeolian fof — A i 5 ft ot shig? jest] 
“9 Sten Ps Po sq. ft. bars | Con Beers ft. sec. | sec. 
0 | 15-00 | 1-0000 |1117 | 761-6 | 1-0000 |2116-2 |1013-2 | 1-0000 | 0-002378 | 1-0000 | 3-719 | 1-564 
1,000 |' 13-02 | 0-9966 |1113 | 759-0 | 0-9644 |2040-9 | 977-1 | 0-9711 ; 0-002310 | 0-9854 | 3-699 | 1-602 
2,000 | 11-04 | 0-9931 |1109 | 756-3 | 0-9298 |1967-7 | 942-1 | 0-9428 | 0-002242 | 0-9710 | 3-679 | 1-641 
3,000 9-06 | 0-9896 |1105 | 753-7 | 0:8963 |1896-7 | 908-1 | 0-9152 | 0-002177_ | 0-9566 | 3-659 | 1-681 
4,000 7-08 | 0-9862 |1102 | 751-0 | 0-8637 |1827-7 | 875-1 | 0-8881 | 0-002112 | 0-9424 | 3-639 | 1-723 
5,000 5-10 | 0-9827 |1098 | 748-4 | 0-8321 |1760-8 | 843-0 | 0-8617 | 0-002049 | 0-9283 | 3-618 | 1-766 
6,000 3-12 | 0-9792 |1094 | 745-7 | 0-8014 |1696-0 | 812-0 | 0-8359 | 0-001988 | 0-9143 | 3-598 | 1-810 
7,000 |+ 1-14 | 0-9756 |1090 | 743-0 | 0-7717 |1633-0 | 781-8 | 0-8107 | 0-001928 | 0-9004 | 3-577 | 1-855 
8,000 |— 0-84 | 0-9721 |1086 | 740-4 | 0-7428 |1571-9 | 752-6 | 0-7860 | 0-001869 | 0-8866 | 3-557 | 1-903 
9,000 |— 2-82 | 0:9686 |1082 | 737-7 | 0-7148 |1512-8 | 724-3 | 0-7620 | 0-001812 | 0-8729 | 3-536 | 1-951 
10,000 |— 4-80 | 0-9650 |1078 | 734-9 | 0-6877 |1455-4 | 696-8 | 0-7385 | 0-001756 | 0-8594 | 3-515 | 2-002 
11,000 |— 6-78 | 0-9614 |1074 | 732-2 | 0-6614 |1399-8 | 670-2 | 0-7156 | 0-001702 | 0-8459 | 3-495 | 2-054 
12,000 |— 8-76 | 0-:9579 |1070 | 729-5 | 0-6360 |1345-9 | 644-4 | 0-6932 | 0-001649 | 0-8326 | 3-474 | 2-107 
13,000 | —10-74 | 0-9543 |1066 | 726-8 | 0-6113 |1293-7 | 619-4 | 0-6714 | 0-001597 | 0-8194 | 3-453 | 2-163 
14,000 |—12-72 | 0-9507 |1062 | 724-0 | 0-5875 |1243-2 | 595-2 | 0-6500 | 0-001546 | 0-8063 | 3-432 | 2-220 
15,000 |—14-70 | 0-9470 |1058 | 721-2 | 0-5644 |1194-3 | 571-8 | 0-6292 | 0-001497 | 0-7933 | 3-411 | 2-280 
16,000 |—16-68 | 0-9434 |1054 | 718-5 | 0-5420 |1147-0 | 549-1 | 0-6090 | 0-001448 | 0-7804 | 3-390 | 2-341 
17,000 |—18-66 | 0-9397 |1050 | 715-7 | 0-5203 [1101-1 | 527-2 | 0-5892 | 0-001401 | 0-7676 | 3-369 | 2-404 
18,000 | —20-64 | 0-9361 |1046 | 712-9 | 0-4994 |1056-9 | 506-0 | 0-5699 | 0-001355 | 0-7549 | 3-347 | 2-470 
19,000 | —22-62 | 0-9324 |1041 | 710-1 | 0-4792 |1014-0 | 485-5 | 0-5511 | 0-001311 | 0-7424 | 3-326 | 2-538 
20,000 | —24-60 | 0-9287 |1037 | 707-3 | 0-4596 | 972-6 | 465-6 | 0-5328 | 0-001267 | 0-7299 | 3-305 | 2-608 
21,000 | —26-58 | 0-9250 |1033 | 704-5 | 0-4406 | 932-5 | 446-4 | 0-5150 | 0-001225 | 0-7176 | 3-283 | 2-681 
22,000 |—28-56 | 0-9213 |1029 | 701-6 | 0-4223 | 893-8 | 427-9 | 0-4976 | 0-001183 | 0-7054 | 3-262 | 2-757 
23,000 | —30-54 | 0-9175 |1025 | 698-8 | 0-4047 | 856-4 | 410-0 | 0-4807 | 0-001143 | 0-6933 | 3-240 | 2-834 
24,000 | —32-52 | 0-9138 |1021 | 695-9 | 0-3876 | 820-3 | 392-7 | 0-4642 | 0-001104 | 0-6813 | 3-218 | 2-915 
25,000 |—34-50 | 0-9100 |1017 | 693-1 | 0-3711 | 785-3 | 376-0 | 0-4481 | 0-001066 | 0-6694 | 3-196 | 2-999 
26,000 |—36-48 | 0-9062 |1012 | 690-2 | 0-3552 | 751-7 | 359-9 | 0-4325 | 0-001029 | 0-6576 | 3-174 | 3-087 
27,000 | —38-46 | 0-9024 |1008 | 687-3 | 0-3399 | 719-2 | 344-3 | 0:4173 | 0-000993 | 0-6460 | 3-153 | 3-177 
28,000 |—40-44 | 0-8986 |1004 | 684-4 | 0-3251 | 687-9 | 329-3 | 0-4025 | 0-:000957 | 0-6345 | 3-130 | 3-270 
29,000 | —42-42 | 0-8948 | 999 | 681-5 | 0-3108 | 657-6 | 314-9 | 0-3882 | 0-000923 | 0-6230 | 3-108 | 3-367 
30,000 | —44-40 | 0-8909 | 995 | 678-5 | 0-2970 | 628-5 | 300-9 | 0-3741 | 0-000890 | 0-6116 | 3-086 | 3-469 
31,000 | —46-38 | 0-8871 | 991 | 675-6 | 0-2837 | 600-4 | 287-5 | 0:3606 | 0-000858 | 0-6005 | 3-064 | 3-573 
32,000 | —48-36 | 0-8832 | 987 | 672-6 | 0-2709 | 573-3 | 274-5 | 0-3473 | 0-000826 | 0-5894 | 3-041 | 3-682 
33,000 | —50-34 | 0-8793 | 982 | 669-7 | 0-2586 | 547-3 | 262-0 | 0-3345 | 0-000796 | 0-5784 | 3-019 | 3-795 
34,000 | —52-32 | 0-8754 | 978 | 666-7 | 0-2467 | 522-2 | 250-0 | 0-3220 | 0-000766 | 0-5675 | 2-997 | 3-913 
35,000 | 54-30 | 0-8714 | 973 | 663-7 | 0-2353 | 498-0 | 238-4 | 0-3099 | 0-000737 | 0-5567 | 2-974 | 4-036 
36,000 |—56-28 | 0-8675 | 969 | 660-7 | 0-2244 | 474-8 | 227-3 | 0-2981 | 0-000709 | 0-5460 | 2-951 | 4-159 
36,090 |—56-5 | 0-8670 | 968 | 660-3 | 0-2234 | 472-7 | 226-3 | 0°2972 | 0-000707 | 0-5451 | 2-948 | 4-172 
36,090 |—56-5 | 0-8670 | 968 | 660-3 | 0-2234 | 472-7 | 226-3 | 0-2972 | 0-0007068 | 0-5451 | 2-948 | 4-172 
37,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-2138 | 452-5 | 216-7 | 0-2845 | 0-0006766 | 0-5334 | 2-948 | 4-358 
38,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-2038 | 431-2 | 206-5 | 0-2711 | 0-0006448 | 0-5207 | 2-948 | 4-573 
39,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1942 | 411-0 | 196-8 | 0-2584 | 0-0006145 | 0-5083 | 2-948 | 4-799 
40,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1851 | 391-8 187-6 | 0-2463 | 0-0005857 | 0-4963 | 2-948 | 5-034 
41,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1764 | 373-4 178-8  0-2347  0-0005582 | 0-4845 | 2-948 5-282 
42,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1681 | 355-8 | 170-4 | 0-2237 | 0-0005320 | 0-4730 | 2-948 5-542 
43,000 |—56-5 | 0:8670 | 968 | 660-3 | 0-1603 | 339-1 | 162-4 | 0-2132 | 0-0005071 | 0-4617 | 2-948 | 5-815 
44,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1527 | 323-2 | 154-8 | 0-2032 | 0-0004833 | 0-4508 | 2-948 | 6-102 
45,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1456 308-0 | 147-5 | 0-1936 | 0-0004605 | 0-4400 | 2-948 6-403 
46,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1387 | 293-6 | 140-6 | 0-1846 | 0-0004390 | 0:4296 | 2-948 | 6-718 
47,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1322 | 279-8 | 134-0 | 0-1759 | 0-0004184 | 0-4194 | 2-948 | 7-049 
48,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1260 | 266-6 | 127-7 | 0-1676 | 0-0003987 | 0-4094 | 2-948 | 7-396 
49,000 |—56-5 | 0-8670 | 968 | 660-3 | 0-1201 254-1 | 121-7 | 0-1598 | 0-0003800 | 0-3997 | 2-948 | 7-761 
50,000 |—56-5 | 0-8670 968 660-3 0-1145 242-2 116-0 0-1523 0-0003622 | 0-3902 2-948 | 8.141 


(55239) 


We. 


/ 
: 


e) 10,000 20,000 h 30,000 36,090 40,000 50,000 


Fic. 1.—Properties of standard atmosphere. 
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Fic. 2.—Variation of viscosity with height. 
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SYSTEM OF AXES 


Symbol | 
ik Designation he es : 
Xes Positive longitudinal lateral normal 
Aibeets forward starboard downward 
: irection 
Force Symbol x yi Z 
M t Symbol L M N 
ea Designation rolling pitching yawing 
Rotation Sue é é ’ 
Volpe | Linear U v w 
y | Angular p q | Y 
Moment of Inertia A B | C 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 

The symbols for the control angles are :— 

€ aileron angle 

n elevator angle 
ny tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


Mass 

Time 

Resultant linear velocity 

Resultant angular velocity 

Density, o relative density 

Kinematic coefficient of viscosity 

Reynolds number, R = /V/» (where / is a suitable linear dimension) 
Normal temperature and pressure for aeronautical work are 15° C 
and 760 mi. 
For air under these ¢ p = 0-002378 slug/cu. ft. 

conditions | » = 1:56 x 10-4 ft.?/sec. 

The slug is taken to be 32-2 lb.-mass. 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio; A = 62/S 

Lift, with coefficient C, = L/4eV?S 

Drag, with coefficient Cp = D/4eV2S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/teV2bS 

Pitching moment, with coefficient C,, = M/deV?cS 

Yawing moment, with coefficient C, = N/}pV7bS 


2. AIRSCREWS 


Revolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient ky = T/om?D4 
Torque, with coefficient ke = Q/on?D> 
Efficiency, 7 = TV/P = Jky/2ak, 


Drag of Circular Cylinders in Flight 
By 
A. V. Strepuens, M.A. 


COMMUNICATED BY THE DIRECTOR OF SCIENTIFIC RESEARCH, MINISTRY OF AIRCRAFT PRODUCTION. 


Reports and Memoranda No. 1892 
gth October, 1940* 


Summary—Reasons for enquiry.—Measurements of the drag of a circular cylinder in flight were required. 


Range of investigation —Measurements of the drag of smooth cylinders 2 in. and 3 in. in diameter of effective length 
5 ft., were made made over a speed range of 100 to 200 m.p.h. in non-turbulent air. The effects of roughening the 
surface were investigated. 

Conclusions—The sharp fall in the drag coefficient characteristic of circular cylinders occurred at approximately 
the same Reynolds number in flight as had previously been found in low turbulence wind tunnels. Roughening the 
cylinder caused this fall to occur at a lower Reynolds number, but the lowest drag coefficient then reached was greater 
than that of the smooth cylinder. 

The following values for the greatest drag per ft. run of smooth and rough cylinders at indicated speeds up to 
250 m.p.h. at 15,000 ft. illustrate the magnitude of these effects. 


Drag per ft. run, lb. 
Diameter, in. 


Smooth. Rough. 
1 14-0 14-0 
13 21*0 12 <7 
. 24-0 12-7 
24 19-0 16-0 
3 14-7 19-2 


The roughness which produced the above results was such as to show the greatest reduction in drag. On a cylinder 
2 in. in diameter it consisted of sharp particles of abrasive approximately 0-001 in. high. It should be emphasised 
that the above results apply only to cylinders in non-turbulent air. 


Introduction.—Much data has been collected on the drag of circular cylinders in wind tunnels » ”. 
In all tests the drag coefficient is found to maintain a steady value of approximately 1-2 up to 
a certain Reynolds number and then to fall progressively to a much lower value. The lower the 
turbulence and the smoother the cylinder the higher does this Reynolds number become. _ It is 
now well established* that no turbulence of an intensity and scale likely to promote transition 
is to be found in the free atmosphere. It therefore seemed possible that the sudden fall in the 
drag coefficient might be delayed to a considerably higher Reynolds number than had been 
found in these experiments, in the case of a cylinder carried on an aeroplane. In fact on the 
basis of the wind tunnel data it was found impossible to predict the drag of a particular cylinder 
some 2 in. in diameter attached to a high performance aeroplane without a possible error of even 
100 per cent. A simple experiment was therefore conducted in flight to provide further data. 


* R.A.E. Report, May, 1940. 
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Method of experiment.—A steel tube 7 ft. long and 2 in. in diameter was mounted so as to 
project vertically out of the gun turret of a Blenheim aeroplane for a distance of 5 ft. It was 
thus well clear of the slipstreams. The tube was rigidly attached to the gun mounting and 
arrangements were made to determine its resistance by measuring the deflection of the tube itself 
at a point 18 in. from the fixed end. 


For this purpose an Ames dial was mounted on a stanchion based upon a stout steel plate welded 
to the bottom of the cylinder. This arrangement is shown in Fig. 1. The apparatus was calibrated 
on the ground by applying known loads to the centre of the protruding part of the cylinder and 
noting the corresponding deflections. The drag of the cylinder in flight could then be obtained 
by reading the Ames dial. Strictly speaking this method involved the assumption that the 
resultant drag occurred at the centre of the exposed cylinder. This could not be guaranteed 
because the air velocity along the cylinder was probably not exactly constant. There was also a 
probability of end effect. In point of fact neither of these uncertainties vitiated the experiments, 
because the cylinder tested was of the same general proportions as that to which it was proposed 
to apply the results. It should be noted that the deflection of the cylinder was taken at a point 
relatively near the fixed end so that the deflection was in effect a measure of the torque, due to 
the resultant wind force on the cylinder, about an axis through the fixed end of the cylinder. 
Beam theory shows that even if the resultant force occurred at a point 6 in. above the centre of 
the protruding part of the cylinder the Ames dial would indicate its torque about the end of the 
cylinder correct to 1} per cent. The use of this crude and convenient method of measurement 
was justified by the results obtained. The maximum and minimum drag coefficients of a smooth 
cylinder agreed reasonably well with the results of previous measurements in wind tunnels, so 
there was no reason to doubt the intermediate values which were of primary interest in the 
present tests. 


The inclination of the cylinder to the vertical in horizontal flight ranged from 1° at an indicated 
speed of 220 m.p.h. to 7$° at 100 m.p.h. Due to downwash effects it is probable that the cylinder 
was approximately perpendicular to the local air stream at all speeds. Owing to the inclination 
of the cylinder to the vertical, small corrections were made to the measured drag to allow for the 
component of the weight of the cylinder in this direction. 


The local velocity at a point 1 ft. in front of the centre of the exposed part of the cylinder was 
measured by means of a swivelling pitot-static head attached to the cylinder itself during one 
flight. Having established the relation between local indicated speed and ordinary A.S.I. reading 
all drag measurements were made by reference to the latter, with the swivelling pitot removed. 


The Reynolds number was calculated on the basis of the local speed at the cylinder. Air 
density and viscosity were obtained from readings of aneroid and air thermometer. The range 
of Reynolds number covered by the experiments was extended by fitting a cylinder 3 in. in 
diameter for the outer five feet of the 2-in. cylinder for some of the tests. 


Results—Smooth cylinder.—Data obtained on the drag of smooth cylinders is given in Fig. 2. 
It will be seen that the nature of the surface was varied during the tests. For the first flights 
with the 2-in. cylinder the steel was simply polished with fine emery paper and left bare. When 
it became apparent that the results were highly sensitive to the nature of the surface the tests 
were repeated with the 3-in. cylinder in a very smooth condition. The surface was filled and 
given a number of coats of cellulose paint in accordance with modern practice for high-grade 
coachwork. A high polish resulted and every confidence could be felt that the cylinder was then 
aerodynamically smcoth. All the results were probably affected to some extent by vibration of 
the cylinder. This was most noticeable on the larger cylinder at the lower speeds. In drawing a 
curve through these points the results for the 2-in. cylinder have therefore been given more 
weight than those for the 3-in. cylinder at Reynolds numbers around 3 x 10°. The curve has also 
been drawn through the set of points giving the higher drag for the 2-in. cylinder, as it is con- 
sidered that the cylinder must have been smoother on this occasion than for the previous flight. 


3 
It had been repolished between the two flights to remove traces of rust. The results of experiments 


in the C.A.T. by Mr. E. F. Relf are given for comparison, and in order to provide a basis for 
extrapolating the flight results to higher Reynolds numbers than were covered by the test. 


Drag of rough cylinder.—During the tests it was noticed that the drag of the 2-in. cylinder 
was greatly reduced when flying through cloud. It was therefore decided to make tests of 
artificially roughened cylinders to see whether practical use could be made of this phenomenon. 
As is usual in experiments of this nature difficulty was experienced in specifying the nature of 
the surfaces tested. The present tests were mainly exploratory and the surfaces will therefore 
be related merely to the processes by which they were produced. 


1. Extreme roughness.—The cylinder was sprinkled with fine sand while the paint was still 
wet and then re-sprayed with paint. The resulting surface resembled a coarse sand paper. 


2. Moderate roughness.—The cylinder was sprayed while wet with a fine silicon abrasive. 
The surface was examined with a microscope and found to consist of closely packed sharp cornered 
particles averaging perhaps 0-001 in. interspersed with isolated particles of 0-003 in. to 0-005 in. 


3. Slight roughness.—The surface was prepared as for (2) but rubbed down with fine sand paper. 
This removed the larger isolated particles, leaving the equivalent of a very fine emery paper. 


4. Stove enamelled.—The cylinder was polished with emery and given one coat of glossy stove 
enamel. Under the microscope some of the scratches due to polishing appeared and the surface 
was found to be composed of “ orange peel’’ markings. 


The curves of drag coefficient against Reynolds number are plotted in Fig. 3 for the surface in 
these four conditions. The curve for smooth cylinders from Fig. 2 is given for purposes of com- 
parison. It will be seen that roughening has two distinct effects, both of which increase 
progressively with the degree of roughness over the range covered. Firstly, the fall in drag co- 
efficient occurs at lower Reynolds numbers and, secondly, the drag coefficient reached after the 
critical change becomes progressively increased. Even though the degree of roughness of the 
cylinders tested has not been precisely defined, the curves taken together give a good idea of the 
lowest drag coefficient that can be obtained by roughening at any Keynolds number between 
210° and dx 10°. 


Application of results —Some consideration is required to apply the results obtained to design 
a cylinder giving the least drag in a practical case. It will be readily seen that over limited 
~ ranges of Reynolds number the drag of a cylinder can be reduced— 
| (a) by increasing the diameter, 
(b) by roughening the surface. 


Both these expedients will, however, increase the drag at Keynolds numbers outside these 
ranges. The maximum effects obtainable by either method -are large and seem to warrant 
investigation. In approaching the problem it will be assumed that it is the greatest drag which 
can be encountered under any probable circumstances of flight which will be of primary concern 
to the designer. This greatest drag will not necessarily occur at the greatest indicated speed to 
be reached. 


Effect of diameter on maximum drag of cylinder up to a particular speed.—In Figures 4 (a) and 
4 (b) curves are given for the actual drag per foot run of smooth cylinders of diameter 1 in. to 
3 in. at indicated speeds up to 350 m.p.h. Values are computed for sea level and an altitude of 
15,000 ft. These curves are based entirely on the single experimental curve of Fig. 2 but present 
the results in a form more suitable for the present purpose. From these curves it is possible to 
read off the greatest drag reached up to a given speed at a given altitude. Values of greatest 
drag have been plotted against diameter in Fig. 5 for indicated speeds up to 250, 300 and 
350 m.p.h, 


+ 


The possibilities of reducing the drag of existing cylinders by increasing the diameter can be 
easily assessed from this diagram. It should, however, be stressed that these results apply only 
to very smooth cylindrical sections in an airstream free from turbulence. For a cylinder in a 
slipstream it would probably not be advantageous to increase the diameter. 


Effect of roughness on drag of cylinder.—From the curves of Fig. 3 it appears that the surface 
designated as “ slight roughness” gives the greatest reduction of drag obtainable by roughening 
at a Reynolds number of 310°. At this Reynolds number the drag of a smooth cylinder is 
reduced by 52 per cent. and some reduction is obtained between 1-6 10° and 3:9 10°. On the 
other hand at 4-5 10° the drag is increased by 80 per cent. by this degree of roughness. In 
order to simplify the problem this roughness will be adopted for purposes of computation. It 
should be remembered that in applying the results to cylinders of size different from the test 
cylinder the scale of the roughness should be considered to be increased in proportion to the 
diameter of the cylinder. 


Taking the curve for slight roughness in Fig. 3, values of the actual drag of cylinders of various 
diameters have been computed for indicated speeds up to 350 m.p.h. and altitudes of 0 ft. and 
15,000 ft. In all cases it was found that the curves of drag against speed did not reach maxima 
at any point in the speed range. In consequence the greatest drag encountered up to a given 
speed is merely the drag at that speed. The values of the actual drags per foot run at indicated 
speeds of 250, 300 and 350 m.p.h. have been plotted against cylinder diameter in Fig. 6 (a) and 
6 (b), for comparison with the corresponding curves from Fig. 5. These figures indicate the 
advantage that can be obtained by “ roughening slightly’ if the diameter of the cylinder and 
the maximum indicated speed are known. The curves also set limits to the maximum drag that 
may be expected for ordinary commercial finishes. 
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Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. ; 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 

The symbols for the control angles are :-- | 

€ aileron angle 

n elevator angle 
n; tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


1. GENERAL 
m Mass 
t Time 
V Resultant linear velocity 
Q Resultant angular velocity - 
e Density, o relative density 


Kinematic coefficient of viscosity 

Reynolds number, R = /V/»*(where / is a suitable linear dimension) 
Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. 
For air under these ¢ 9 = 0:002378 slug/cu. ft. 

conditions \ » = 1:56 x 10-4 ft.?/sec. 

The slug is taken to be 32:2 lb.-mass. 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio, A = R/S 

Lift, with coefficient C, = L/4eV°S 

Drag, with coefficient Cp = D/4pV?S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/4eV?0S 

Pitching moment, with coefficient C,, = M/4eV?cS 

Yawing moment, with coefficient C, = N/$eV70S 


2, AIRSCREWS 


W 


i (OB ay Be © Hey i 


Revolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient ky = T/pn?D? 
Torque, with coefficient ke = Q/oen?D® 
Efficiency, 7 = TV/P = Jky/2akg 
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Summary.—This report gives a method of transforming a cascade of straight blades into a 
cascade of aerofoils. As the general flow past a cascade of straight blades is known! 2 the 


general flow past the cascade of derived aerofoils can be obtained. An illustrative example is 
given. 


1. Introduction.—lf any particular flow past a series of oval bodies on an axis of a complex 
plane is known, a transformation can be found which converts the ovals into a cascade of straight 
blades, since the corresponding flow past the latter is known. Moreover, by applying the reverse 
transformation to the general flow round the straight blades, the general flow round the ovals 
is obtained. The direct transformation can also be applied to an offset series of ovals, which 
then transform into a cascade of aerofoils: the process corresponds to the usual transformation 
of a single circle into an aerofoil. This method has been used before by one of the writers’. 
The present report applies the same method to a series of ovals which are in a physical sense 
much simpler than those previously used and which yield aerofoils approximating more closely 


__to those met with in practice. An illustrative example is given. 


2. Outline of argument.—The method employed may be briefly summarised as follows :— 


(i) The normal flow past a series of doublets on the imaginary axis is known (it has been 
given by Lamb‘) and can be interpreted as the normal flow past a series of oval bodies. 


(ii) The general flow past a series of laminae lying along the imaginary axis is known ; 
in particular the normal flow is known, and a transformation can therefore be obtained which 
converts the laminae into the ovals discussed under (i). 


(iii) Application of this transformation to the general flow round the laminae gives the 
general flow round the ovals. 


(iv) A particular case of the general flow round the ovals is that for which the flow at 
infinity is inclined to the cascade axis, but for which there is no circulation. In this case 


the ovals can immediately be transformed into a cascade of flat plates parallel to the direction 
* of flow at infinity. 


(v) Application of the transformation of (iv) to an offset series of ovals yields a cascade 
of aerofoils. If the offset ovals pass through the points corresponding to the trailing edges 


of the inclined flat plates (which are stagnation points in the flow round the transform ovals) 
the aerofoils will have cusps at their trailing edges. 


2 


3. General Flow past a Cascade of Ovals.—The flow due to a series of uniform doublets lying 
along the imaginary axis in the /-plane with a period az is known to be# 
WMG coe, i 


where C is a parameter defining the strengths of the doublets. If a uniform flow of velocity U 
parallel to the real axis is added the flow becomes 


w= Ul + C.coth J, 
== U4) + sinh*6coth 7} =. e e ae ry + (1) 


where for convenience C is replaced by U sinh?8. Equation (1) gives the velocity potential and 
and stream function respectively as 


eS sinh?8 sinh 2m } 
0 Feu lew ts cosh 2m — cos 2n} ’ 


sinh? sin 2 | 
cosh 2m — cos 2n J ; 


v= U |n _ 
The streamline | = 0 evidently consists partly of the axis 7 = 0 and partly of the closed oval 


me AY 
cosh 2m = cos 2n 4 aaa rk cz. 4) fe Pa fe (2) 


This oval has its centre at the origin, and repeats at intervals of a along the imaginary axis. 
The lengths of the semi-diameters are defined by the equations 


m =f, 
nm tan n = sinh?8. 
Equation (1) may therefore be regarded as giving the normal flow past a series of oval bodies, 


of shape given by (2), immersed in the fluid. The oval corresponding to B = 1-0 is shown in 
bi tau & 


Differentiation of equation (1) gives the velocity components of the flow past the ovals in 
the form 
EP GE VAG todd sinh*6 | 
a Ger sinh’ |’ a. Pe oy < xi a. (3) 


so that the stagnation points of the flow are given by / = + $+ ria and therefore, as is other- 
wise obvious, lie at the ends of the major axes of the ovals. 


Now consider the flow in the ¢-plane given by 


wr U-coshesisec a Cosi) seam om oe y a es, (4) 
or 
dw U sinh Z 
ae Tas ais HM t. se a Sg rr (5) 
SG ,/sinh?t + sin? a 


It represents? ? the flow normal to a series of laminae of length 2a distributed along the 7-axis 
with a period of z. Equation (5) shows that there are stagnation points when sinh f vanishes, 
and that the velocity of flow at infinity is U, as in the case of equation (1). Elimination of w 
between (1) and (4) therefore yields a transformation which converts the laminae of (4) into the 
ovals of (1) ; the transformation is 


cosh-1 (sec a cosh C) = / + sinh? 6 coth J. ie ie oy oF ae ee. (6) 


In addition, the stagnation points in the two flows must correspond; accordingly, we put 
cosh ¢ = 1, / = + 6 in (6), which becomes 


sec a = cosh { 8 + sinh* 6 coth 8}, a ie 3 fe a (7) 


& 


wr 
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a relation which connects the two parameters a and 8. It will be convenient in what follows to 
employ the two auxiliary quantities 


y = 8+ sinh?6 coth 6, id Lé, gf Ry iy if (8) 
A = 1 + sinh?6 coth J. cs i ge eg ‘ ih (9) 

The transformation (6), connecting the /-plane ovals and the @-plane laminae, then becomes 
Biot 7 cost © =eosit A.7* bs a i. on nf ae A, (10) 


It may be noted that (7) and (8) show a to be the gudermannian of y, so that the circular 
functions of a, 1f required, are given by 
Siting seotanilir ye: 
COS S=(SeChiy F's): bt. = i, be, oe uN sarong E} 
tan a = sinh y 


The general flow past the ¢-plane laminae, of which (5) is a particular case, is 
dw Usinh¢+72W cosh Z 


— = : : 9 
at a/ sinh? ¢ va one e uv. sor as = re ack (12) 


Here V is the velocity parallel to the y-axis, and there is circulation 27W round each lamina. 
We now apply the transformation (10) to the flow (12). The laminae transform as before into 
the /-plane ovals, and the general flow round the ovals becomes 


sinh2 / 


(13) 


dw [u EW. cosh, ha sinh A 1 sinh? 6 
a in } 


When W = V = 0, this agrees with equation (3). 


4. Transformation of the Ovals into Inclined Flat Plates and Aerofotls—As a particular case of 
(13), put W = 0 and V= Utano. The equation then becomes 


dw a tan o sinh A ipa sinh? 6 
eae pemsaeerat | 

/ sinh? \ — sinh? v ea. sinh? / 
{n this case the flow past the ovals is inclined at an angle c to the real axis at large distances on 


either side of the cascade. Now imagine the flow pattern of (14) to be replotted in the z-plane, 
with the streamlines all straightened out. The flow in the z-plane is obviously 


dw 

Se ecuatlows P talc). ft: i Ue Y ei; oy -. — (1S) 
Since the flow in the z-plane is uniform and rectilinear throughout, it is evident that the ovals 
have transformed into flat plates lying along the lines of flow. 


(14) 


Elimination of w between (14) and (15) gives the transformation connecting the /- and z-planes 


as 
dz ae {cos F woe o vane Kast | . a (16) 
dl r/ sinh?, — sinh?y J sinh*/ 
or 
z = e'°{) cos o — i sin o cosh—! (sech y cosh 4) }... = onl) 


At first sight, it would appear that (16) vanishes for / = 8. However, when/ = 6, A = yalso, 
and it is found that in this case dz/dl is neither zero nor infinite. Accordingly, the only infinities 
of (16) are for sinh / = 0; these points are the centres of the ovals. The zeros of (16) are given by 

7 sino sinh A 
cos 6 = SG, 
a/ sinh?) — sinh?y 


or 
sinhA=-+coso sinhy. .. meg Be: "¢ Me hie eitiln «Ohee) 


(56861) At 
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It is evident that when (16) vanishes, (14) also vanishes, and the singularities (18) therefore 
coincide with the stagnation points in the flow (14). These correspond to the leading and trailing 
edges of the flat plates in the z-plane. By substitution from (18) in (17), these points are given as 


z= +e? {cos o sinh (cos o sinh y) + sin osin“ (sin otanhy)} .. ey guts) 


Equation (17) is thus a transformation from the /-plane to the z-plane which converts the 
cascade of ovals into a cascade of straight blades inclined at an angle o to the real axis. Now, 
in a manner analogous to the transformation of a circle into Joukowski aerofoils, a cascade of 
offset ovals can be transformed by (17) into a cascade of curved aerofoils. The process is con- 
veniently illustrated by the example of § 5. It may be noted that, for computational purposes, 
it is convenient to draw the transformed aerofoil in the Z-plane, where Z = ze~’*, since a factor e?? 
is thereby avoided. The Z-plane axes are rotated through o relative to those of the z-plane ; 
and the transformation (17) is simply 


Z = cose — 1 sino cosh-} (sech y cosh 4). = fe. «se @A(20) 


In the application of the transformation to offset ovals it sometimes happens that the resultant 
Z-plane contours cross themselves. Usually this is immediately apparent, but occasionally such 
a cross-over may occur close to a cusp. Care must be taken to determine whether there is a 
cross-over or not, since it is obviously impossible to associate a fluid flow with a contour which 
crosses itself. 


5. Illustrative example.—Take the series of ovals defined by 8 = 1-0, a member of which is 
shown in Fig. 1. We shall transform this into a cascade of straight blades at 30° to the real 
axis, so that o = z/6. We then have the following data :— 


Cesare) 
i ) 
o= 2/6, } (21) 
and substitution in (8) gives 
pe eA OI ESS oe be 8 x: re 7 ay ee A's) 

The stagnation points (18) now give 

eee / Laas re es ae e . Me i eee 
and substitution for ’ from (9) gives an equation for / of which the solution is found to be* 

+1 = 0-943 + 0-3237. (5 a F ae ¥, w (24) 


Moreover, substitution of the numerical value of y in equation (19), and replacement of z by Ze, 
gives the ends of the Z-plane flat plates as 


Le 13. ne es = S e A me Se ele 
The Z-plane flat plate is shown in Fig. 2. 


Now consider the ovals given by 8’ = 1-1. _ If these ovals are offset so that they pass through 
one of the points (24) (and the corresponding points along the series) and so that they enclose the 
other zeros and the infinities, they will transform into a cascade of aerofoils with cusps at their 
trailing edges. In Fig. 1 the oval corresponding to 6’ = 1-1 has been displaced so that its centre 
is at the point 


1 = 0-139 — 0-131 2, fs ae “4 ae a: <i 4+ (28) 
and it then passes through one of the singularities (24), namely 
L= — 0-943 — 0-323 7. A, & = a ea a GHZT) 


* The solution of this equation may be obtained by finding an approximate value of / and applying Newton’s method 
for a closer approximation. The co-ordinates are of course related by equation (2). 


x 


fe) 


There is no difficulty in the computation of the Z-plane transform ; any point / on the displaced 
oval is substituted in (9) with 6 = 1-0 to find A, and this is then substituted in (20) to find Z. 
The Atlas Supplement to Kennelly’s Tables of Complex Hyperbolic and Circular Functions® was 
found most convenient for this purpose. 


6. Properties of the Derived Cascade, and Determination of an Equivalent Cascade of Flat Plates.— 
The flow at infinity is the same for both the @’ ovals and the derived aerofoils ; accordingly the 
properties of the cascade can be determined from those of the ovals. If the Kutta-Joukowski 
hypothesis is assumed, the velocity of flow at the aerofoil trailing edges is finite ; and therefore 
there is always a stagnation point at each of the corresponding points on the @’ ovals. 


Now in Fig. 1 we may regard the 8’ oval as being transformed with respect to the 6 oval : 
the transformation has a singularity at one of the points of intersection of the two ovals, and 
this is the stagnation point of the flow round the 8’ oval. Suppose we now select a different 
8 oval (still with its centre at the origin) and displace the original 8’ oval, and the flow pattern 
round it, so that the stagnation point on the 8’ oval lies at the new intersection of the two ovals. 
With this point as a singularity we may transform the 8’ oval with respect to the new 8 oval, 
and shall obtain a cascade which has the same properties as the 8’ oval and, therefore, as the 
original cascade. An infinite number of cascades having the same properties may be obtained 
in this way. In particular, when the 8’ oval coincides with the 6 oval, so that its centre is now 
at the origin, the process of transformation yields a cascade of flat plates, which may be described 
as equivalent to the original cascade. 


From (26) and (27) the stagnation point, when the 8’ oval is not displaced, is 
1 = — 1-082 — 0-1927. — a e. - Zz - dent 2a) 


Accordingly, this point and its reflexion in the origin are the singularities of the new trans- 
formation ; and they are defined by equation (18). In this equation, y is defined by (8) with 
6 = 1-1 and A by (9) and (28). Accordingly (18) gives o, and when these values of y, 4, and o 
are inserted in (19) the half chords of the flat plates are derived. The results are 


ies Ole 

N= 282 | 

Goa -S0SsTadian (29) 
alee 

Li Oe, 


The equivalent cascade of flat plates is thus inclined at 17-4° to the real axis, and has a gap- 
chord ratio 
S 7 


c 6-444 
The properties of cascades of this type are given in Reference 1, and define the properties of the 


curved aerofoils of Fig. 2. It may be remarked that the angle of the flat plates to the real axis 
(17-4°) is the angle of flow for no lift on the cascade of Fig. 2. 


= 0-488. ae a eee eet) 


7. The writers wish to express their thanks to Messrs.Metropolitan-Vickers Electrical Co., Ltd., 
and to the Aeronautical Research Committee for permission to publish this paper. 
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SYSTEM OF AXES 


Symbol =|, y : 
~ Axes tee longitudinal lateral normal 
recon, forward starboard downward 
Force- Symbol x an Z 
Symbol L M N 
Mone) Designation rolling pitching yawing 
i ; Angle of 
| Rotation Mee ae . d 
Telopt Linear u v w 
y Angular Bova q f 


Moment of Inertia | A | B C 
| | 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative’ moment about the corresponding axis. 


The symbols for the control angles are :— 
2 é aileron angle 
,) 3 n elevator angle 
ny tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


1, GENERAL 


m Mass 
i Time 
V Resultant linear velocity 


Q Resultant angular velocity 
0 Density, o relative density 
v Kinematic coefficient of viscosity % 
R Reynolds number, R = /V/» (where / is a suitable linear dimension) 

Normal temperature and pressure for aeronautical work are 15° C 

and 760 mm. 

For air under these ¢ 9 = 0:002378 slug/cu. ft. 

conditions | v = 1-56 x 10-4 ft.?/sec. 
The slug is taken to be 32:2 lb.-mass. 


OL Angle of incidence 

é Angle of downwash 

3 Area | a 

b Span tb 
C Chord x 
A Aspect ratio, A = 0?/S 

di, Lift, with coefficient C, = L/4eV?S 

D Drag, with coefficient Cp = D/ZeV*S 

y Gliding angle, tan y = D/L 

rE, Rolling moment, with coefficient C,; = L/4eV7bS 

M Pitching moment, with coefficient C,, = M/4eV?cS 

N  Yawing moment, with coefficient C, = N/$eV20S 


2. AIRSCREWS 


nt Revolutions per second be 
D Diameter 

J V/nD 

P Power 

-. Thrust, with coefficient ky = T/pn?D4 

O Torque, with coefficient kg = Q/on?D° 

7 Efficiency, 7 = TV/P = Jkr/22ka 
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Compression Tests on Curved Plates of 
Thin Sheet Duralumin 
By | 
Bie lwcGox, M.A AND. W.-]. CLensnaw, B.Sc., A.Inst.P., 
OF THE ENGINEERING DEPARTMENT, N.P.L. 


Reports and Memoranda No. 1894 
28th November, 1941 


Summary—(a) Purpose of investigation.—The results of compression tests previously carried out on flat plates showed? 
that the stress-strain relationship of such plates could be represented approximately by two straight lines. The present 
tests were undertaken in order to investigate the stress-strain relationship for curved plates and to determine whether 
in this case a similar simplification of the stress-strain relationship was permissible. 


(0) Range of investigation.—Tests have been made on plates of duralumin having the following dimensions :— 


0-067 in. thick 

0-038 - 

0-030 PE curved to 18 in. radius 

0-024 . 

0-017 “4 4, 6, 8, 10 and 12 in. wide (circumferentially). 


0-067 __,, 
G033: 40"), 


All the plates were 24 in. in length parallel to the axis of curvature and all the edges were so held that they could 
be considered clamped. Compressive loads were applied to the curved edges of the plates (parallel to the length) and 
the loading conditions were such that the planes containing these curved edges were maintained parallel. The 
displacements in the direction of the load were measured over the whole plate (including the end fixings) and an amplified 
record was obtained by means of a hydraulic and mechanical transmission. The use of this transmission in conjunction 
with the automatic load. recorder fitted to the testing machine enabled load- deformation curves to be drawn 
autographically. 

(c) Conclustons.—Curved plates subjected to compression parallel to the axis of curvature deform in a similar manner 
to flat plates. Below the buckling load the plate behaves as a solid test piece, the compressive stress. being uniform 
over the whole plate. After buckling the central part of the plate is stressed less heavily than the edges (which are 
constrained to remain straight) and, since the strain is proportional to the stress in the edges, the stiffness (rate of 
increase of load/rate of increase of strain) of the plate is reduced. 

In contrast to flat plates, buckling of curved plates may result not only in a reduction of stiffness but also in an 
actual reduction of load. The reduction is greatest in thick, wide plates curved to a small radius ; but the dimensions 
of plates for which the effect is serious probably le outside the range of normal aircraft engineering. The elastic 
behaviour of the thinner, narrower plates may be represented fairly closely by two straight lines. The buckling stress 
of a curved plate is greater than that of an equivalent flat plate ; but the reduction of stiffness after buckling is also 
greater. 

No theoretical analysis of the buckling and behaviour after buckling of curved plates with clamped edges is available ; 
but on the basis of the experimental results, empirical formulae for the buckling stress and for the reduction of stiffness 
after buckling have been developed. In choosing these formulae, use has been made of all the theoretical results 
available 1» 2. 


seurved to 36 in. radius ] 


Introduction.—In R. & M. 1554 an approximate theory was developed describing the buckling 
and the behaviour after bucklingt of a flat rectangular plate of thin sheet material subjected to 
compression parallel to its longer edges. It was shown that, owing to the formation of waves at 
buckling and the consequent reduction of stress in the centre of the sheet, the effective stiffness 


* This work was carried out during 1934. 
+ The edges of the plate being supported, buckling need_not lead at once to complete collapse. 
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of the plate after buckling should be reduced ; that the ratio of the stiffness immediately after 
buckling to the stiffness before buckling should be 4 if the edges of the plate were simply supported 
and about 4 if the edges were clamped ; that a further decrease in stiffness should occur as the 
load was increased above the buckling load ; and that, if total collapse occurred at some definite 
value of the strain (i.e., edge stress), the ultimate load carried P,, should be given by the equation 


P,, = ¢(At + Ba), 


A and B being constants for any given material with given edge conditions, and the term A? 
being considerably larger than the term Bd. The results of tests carried out in America® on . 
plates with simply supported edges were shown to agree reasonably well with the theoretical 
conclusions, and the results of a series of tests carried out at the N.P.L.1 on plates with clamped 
edges also confirmed the main theoretical conclusions. 


Although for some purposes, formulae giving the ultimate load a plate can carry may 
be sufficient, such formulae require the specification. of the value of the limiting edge 
stress. If the value of this stress is in doubt, it may be more convenient to refer directly to the 
stress-strain relationship of the plate. Thus methods can be developed for determining from 
a simplified stress-strain diagram for the plate in compression the behaviour of stiffened tubes 
and quadrants in bending.4 Treatment by this method has the advantage that the theoretical 
prediction of the behaviour can be compared with the experimental results not only with regard 
to the value of the collapsing load but also with regard to the deformation at lower loads. 


The applicability of this method depends, of course, upon the accuracy with which the stress- 
strain relationship of the plate can be represented by a simple diagram ; but since in any practical 
construction the stiffeners will carry the major portion of the load, considerable inaccuracy in 
the representation of the behaviour of the plate may cause relatively little error in the estimate 
of the behaviour of the complete structure. Thus, for flat plates the simplified stress-strain 
diagram consisting of two straight lines is probably sufficiently accurate for all practical purposes ; 
but, before a similar diagram could be assumed to represent the behaviour of curved plates it 
became necessary to investigate whether for these plates this assumption was justified. 


Theoretical Considerations.— The only available theory relating directly to the behaviour of 
curved plates in compression is the calculation by the Ritz energy method of the buckling stress 
of a plate with simply supported edges.? In addition.as limiting cases we have the approximate 
results of R. & M. 1554 for the buckling (and the behaviour after buckling) of flat plates and the 
exact calculation of the buckling stress of a complete circular tube in compression given by 
Professor Southwell’s general theory of elastic stability.® This latter result is incorporated as 
a limiting case in the results of reference 2. 


In practical stiffened constructions, the edge conditions of the individual plates approximate 
more closely to clamping than to simple support. Since the limiting condition of the complete 
tube is presumably independent of edge conditions, it would be reasonable to assume that the 
‘buckling of a plate with clamped edges might be defined by the formula of reference 2 with 


Ty? 


22 
the term = replaced by 


5 Comparison 


of the experimental results with this amended formula is made below. 


With regard to the behaviour after buckling, it might be anticipated that, since the develop- 
ment of the buckled form must to a large extent destroy the essential curvature of the plate, 
this curvature should have relatively little effect. This does not imply that the reduction of 
stiffmess after buckling should be independent of curvature; but rather that the reduction 
should be slightly greater for curved than for flat plates, so that, in conjunction with the higher 
buckling stress for curved plates, the mean stress/maximum stress ratio at maximum stresses 
well above the buckling stress should be effectively the same for curved as for flat plates. 
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From the description of the experimental results below it will be apparent that the buckling 
of curved plates is rather different in type from the buckling of flat plates. The effect of this 
difference on the effective description of the behaviour of curved plates by simplified stress-strain 
diagrams will be discussed after description of the experimental results. 


Description of Apparatus.—The apparatus used for supporting the panels in the present tests 


2-is shown in Fig. 1. The top and bottom (curved) edges of the plate to be tested were clamped by 
4 curved cover plates to the curved faces of angle pieces bolted to the top and bottom plattens. 
x The longer (straight) edges of the plate were guided in slots 3 in. deep, cut in the edges of a series 


of strips supported upon two pairs of pillars. These strips were of 16 S.W.G. duralumin and were 
spaced along the supporting pillars by rubber washers } in. thick. The whole ‘“‘ comb ”’ of strips 


© on a pair of supporting pillars formed one unit which could be turned back to front or end for end. 
= Two slots of different widths were cut in each edge of each strip slightly off the centre line, so that 
* the four possible positions of the ‘combs ” brought four different widths of slot into alignment 
.. with the top and bottom angles. The lower ends of the supporting pillars of the combs weré 
= supported in cylindrical cups fixed to the upper face of the lower platten, and their upper ends 


were braced laterally by bars connecting one pillar of each pair to the corresponding pillar of the 
other comb. A number of different positions for the cups supporting the pillars were provided | 
to allow plates 4, 6, 8, 10 and 12 in. in width to be tested. 


All the tests were carried out in the 25-ton Denison testing machine installed in the Engineering 
Department of the National Physical Laboratory. In this machine the upper and lower straining 
heads are constrained to remain parallel, any provision for axial loading being made independently. 
In the present tests the plate apparatus was supported on the lower head of the testing machine 
and the load was applied to the upper platten through a cylindrical seating, 3 in. in diameter, 
of which the axis was perpendicular to the plane of the centre of the plate. This seating accom- 
modated any lack of parallelism between the upper platten and the upper head of the testing 
machine in the lateral direction, while the slight flexibility of the upper platten on the plate 
accommodated any difference in the perpendicular direction. So soon as any appreciable load 
was applied both these adjustments became inoperative and the plate was subjected to uniform 
compressive strain. 


The apparatus used for recording the strain of the plates is shown in Figs. 1 and 2. To the 
~ lower head of the testing machine, at the four corners of the lower platten of the plate apparatus, 
“were fixed four cylinders containing oil, each cylinder being # in. in diameter. These cylinders 
‘were closed by pistons, to which the motion of the top platten of the plate apparatus was 
transmitted by ball ended rods bearing against conical depressions in the tops of the pistons and 
~in the under side of the upper platten. All four cylinders were connected by stout copper pipes 


~to a single recording cylinder, about + in. in diameter, the bore of this cylinder being straight 


‘and true for a length of 4 in. The motion of the piston closing this cylinder was thus about 


~40 times the relative motion of the top’and bottom plattens of the plate apparatus. 


_ The Denison testing machine used in the present tests is fitted with an autographic recorder. 
« The load is weighed through a multiple-lever system by a spring, two springs, one of 25 tons 


sand the other 5 tons capacity, being provided. The motion of the beam to which the spring is 


“attached, rotates, by means of a system of pulleys and wires, a drum upon which a piece of . 


«graduated paper is fixed. A bar carrying a pencil is traversed axially over this drum on two 


pulleys, one of which is actuated by the strain recording mechanism. In the present tests, 
the pulley moving the pencil was rotated by a wire connected at one end to the piston of the 
recording cylinder, passing round the pulley and carrying a weight of about one pound at the 
free end, a mechanical magnification of 24 to 1 being obtained. The motion of the pencil was 
therefore about 100 times the relative motion of the upper and lower plattens of the plate 


“apparatus. : 
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The method of supporting the plate* proved very satisfactory throughout the tests. The 
support afforded to the edges of the plate. was sufficient to prevent buckling at the edge, even 
of plates 0-018 in. in thickness stressed to 12 tons/sq. in. In most cases the plates slipped quite 
freely through the slots in the strips of the comb; but, even where, owing to severe local 
buckling, the plate seized in the slots, the rubber spacing washers prevented any considerable 
variation of stress being set up. On the other hand, the appearance of the buckled plates 
showed quite clearly that the edge support was sufficient to maintain the edges of even the thickest 
plates (16 S.W.G.) in the original form. 


The strain recording apparatus also worked satisfactorily. The lost motion was almost 
negligible and there was no other source of appreciable hysteresis. No. particular provision had 
been made for preventing oil from leaking past the pistons; in fact, although the pistons all 
fitted the cylinders fairly closely, reliance was placed upon the leakage (under the pressure of 
about 20 Ib./sq. in.) to maintain effective lubrication. Owing to this leakage there was a slow 
movement of the strain recorder with time ; but this movement was far too slow to have any 
measurable effect in the period of a normal test. 


The load recording device unfortunately did not work so satisfactorily. For the thicker plates, 
for which loads as high as 5 tons were in certain cases required, it was possible to use either the 
5-ton spring of the recorder with a low mechanical magnification or the 25-ton spring with a 
correspondingly higher mechanical magnification—either system giving a good positive drive 
to the recording drum ; but, for the thinner plates, for which much smaller loads were required, 
it was necessary, in order to obtain a record of reasonable size, to use both the light spring and a” 
high magnification. In such cases the actuating torque on the recording drum was very small 
and great care had to be exercised to reduce the effect of friction in retarding the-motion of the 
drum. In spite of all precautions, however, the lost motion due to friction was in some cases 
an appreciable part of the total motion. The accuracy of the results of tests given below is for 
this reason lower for-thin plates than for thick ones. 


Materials Tested.—Part of the material used in the present tests was the remainder of the 
duralumin sheet (GOK) previously used in the tests on flat plates described in reference 1. 
A further supply of duralumin sheet (HKZ) was subsequently obtained. The surface of the 
material HKZ was anodically treated, that of the material GOK was untreated. 


Method of Test.—Preliminary tests on fairly thick wide plates showed that the behaviour of a 
plate on first loading might be appreciably different from its behaviour after it had once been 
buckled. For this reason, except for the first few plates, each plate was loaded and unloaded 
three times. In the first test the plate was loaded until it buckled, in the second and third it 
was loaded up to an edge stress of about 10-12 tons/sq. in. The load was in most cases applied 
at a rate of about 2-3 tons/sq. in. edge stress per minute, a complete test on one plate taking 
about 45 minutes. 


Results of Tests.— For each plate tested a complete record of total load against overall strain 
was obtained autographically. Unfortunately, it is impracticable to reproduce all these diagrams ; 
but typical diagrams together with sketches of the buckled forms are shown in Figs. 3-14. 


Fig. 3 shows the results of tests on the plate GOKICB1. On first loading this plate buckled 
with considerable noise into one central wave at a load of 5:5 tons (stress 7-6 tons/sq. in.). (So 
violent an oscillation of the recorder was set up that the drum slipped; but the amount of slip 
was subsequently. determined from the final position of the drum after the plate had been 
unloaded.) Immediately upon-buckling, the strain increased to that corresponding to a load 
in the unbuckled state of 6-7 tons (stress 9-3 tons/sq. in.), while the actual load carried fell to 
2-8 tons. On unloading, the load decreased to just over 1 ton at a strain corresponding to 


* The method is a slight modification of that developed at the R.A.E. by the late Mr. Gerard for supporting the 
edges of a semi-circular tube when compressed parallel to its axis. 
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a load in the unbuckled state of 2-5 tons; then the buckle suddenly disappeared, the load 
increasing to about 1:8 tons. During the remainder of the unloading the plate behaved as a 
solid test piece. On reloading the buckle began to develop at a load of 2:1 tons (stress 2-9 tons/sq. 
in.) ; but did not assume its final form until the edge stress was about 4:5 tons/sq. in. Then 
the load, which had remained practically constant since buckling began, fell suddenly to about 
1-5 tons. Upon further loading the mean stress increased at about half the rate of the edge 
stress (the strain). Upon unloading an increase of load from 0-5 to 1-0 ton as the plate returned 
to the unbuckled form was again observed. The behaviour of the plate in the third loading was 
very similar to its behaviour in the second loading. . Buckling began at a slightly lower load 
(1-8 tons) but again the buckle took some time to develop. Once again, an increase of load as 
the buckle disappeared was observed. 


The behaviour of the plate GOKICBI has been described in some detail because it is typical 
of the behaviour of many other plates of similar dimensions.. As a contrast to this type of 
record, Fig. 4 shows the results of tests on the plate GOKICE1. In this narrow plate, the process 
of buckling was similar to that of a flat plate. The reduction of stiffness at buckling was well- 
marked ; but there was no actual reduction of load, and the reduction of buckling load on 
reloading was relatively small. In general, the sudden changes of load became less marked 
as the width of plate was decreased, as the thickness was decreased or as the radius ot curvature 
was increased ; but the variation of the magnitude of the changes with width was much greater 
for thick than for thin plates. 


For comparison, in particular with Fig. 3, Fig. 5 shows the results of the tests on the plate 
HKZ2CB1. In this diagram the points chiefly to be remarked are the sudden fall of load (due to 
a new buckle) that occurred at an edge stress of about 12 tons/sq. in. and the nature of the 
buckled form finally developed. For plates of which the “ bulge’’ was appreciable this form 
of buckling appeared to be the normal one. Narrower plates usually buckled into one or more 
central waves in the manner of flat plates; but the occurrence of central waves in wide plates 
(as for instance in GOKICBI, Fig. 3) may be regarded as abnormal. Another fairly good example 
of the “ side wave ’’ method of buckling is shown in Fig. 6 (HKZ2CC2). 


A diagram obtained in tests on the plate HKZ3CBI1 is shown in Fig. 7. This thickness of 
material (0-030 in.) was fitted into the slots used for the thicker materials HKZ2 and GQK2 
(0-038 in. thick). In consequence, the edge support was not so good as in the other tests and 
there was a tendency (indicated in Fig. 7) for the edges of the plate to buckle. This effect probably 
to some extent reduced all the loads carried by the plates of this thickness. 


Typical diagrams for still thinner plates are shown in Fig. 8 (GOQK8CB1), Fig. 9 (@QK4CB1) 
and Fig. 10 (GOK4CD1). As might be expected, plates curved to a larger radius behaved more 
nearly like flat plates. Examples are shown in Fig. 11 (HKZ2DB1) and Fig. 12 (HKZ1DC1). 


Analysis of Results —-From the diagrams, Figs. 3-12, it will be apparent that theoretical 
description of the behaviour of the plates is likely to be difficult. The fact that the buckled 
form is developed suddenly with simultaneous increase of strain and decrease of load suggests 
that the unbuckled state may be in unstable equilibrium at loads well below that which actually 
causes the plate to buckle. Moreover, in the thicker plates, the sudden deflections due to 
_ buckling are very large and may well.cause permanent deformation ; so that, in subsequent 
_tests, the behaviour of-the plate may be affected by damage done in the first loading. 


The only figure ‘therefore that is likely to have any real theoretical significance is the value 
of the primary buckling stress on first loading. The values of this stress deduced from. the 
experimental results are given in the sixth column of Table 1 and in the ue column are given the 
values calculated from the formulae :— 


: [42 /4 {ia | i 
Buckling Stress = E ape +  /9— + — } ue ri SF (1) 


E 
1.— o? 


E = Young’s modulus | 


where EF’ = 


- for the material 


o = Poisson’s ratio | 

¢ = thickness 

6 = width (circumferentially) of the plate. 
y = radius of curvature 


Apart from the comparison of the values of the primary buckling stresses, there remains the 
practical problem whether it be possible to devise empirical formulae that shall describe the 
behaviour of the plates with accuracy sufficient for the ordinary purposes of design. 


The first question that arises in this problem is whether attention should be concentrated upon 
first loading or upon some subsequent loading. In some cases (cf., for instance Figs. 4, 8 and 9) 
the choice makes but little difference ; but in others (cf., for instance Figs. 3 and 5) the difference 
is considerable. This question should be decided by consideration of the severity of the loading 
condition which, in practical constructions, would cause the constituent plates to buckle. If 
this loading were so severe that it was never likely to occur in practice, then first loading might 
be adopted. If it represented ordinarily severe service conditions then perhaps second loading 
might be used ; whilst, if normal conditions carried the plates well past buckling, third or some 
subsequent loading might be more suitable. 


It may be said at once that, if it is required to represent the behaviour of a plate after an 
indefinite number of fairly severe loadings (up to 10 tons/sq. in. edge stress, say), the best course 
is to assume simply that the load is proportional to the strain, the stiffness being, of course, only 
a fraction of the elastic stiffness. This assumption is equivalent to that made by Sechler and 
others that the plate is reducible to an “ equivalent width ”’ ; but it does not necessitate Sechler’s 
second assumption that the equivalent width is independent of the actual width. 


On the other hand, edge stresses of the order of 10 tons/sq. in. are not likely to be set up 
frequently under working conditions, and, if plates are designed on the “ chord stiffness ”’ from 
zero to 10 tons/sq. in. edge stress, the strength at lower edge stresses may in many cases be 
seriously underestimated. Figs. 13 and 14 show load-deformation curves obtained in repeated 
loading firstly up to about 5 tons/sq. in. edge stress (curves marked 1, 2 and 3) secondly, up to 
10 tons/sq. in. (curves marked 4 and 5), and thirdly up to about 13 tons/sq. in. (curves marked 
6 and 7). It will be seen that the repetition over each range is fairly good and that the behaviour 
at values of the edge stress in the neighbourhood of 5 tons/sq. in. is represented better by the 
empirical formula than by the “ reduced stiffness ”’ line. 

It is suggested therefore that second loading is in general the best criterion. This indicates 
1s oe of a plate that has been buckled but which has not been stressed very fat above 
uckling. 7 

Figs. 3-14 show that the behaviour of many plates may be represented fairly accurately by. 
two straight lines, i.e., by the elastic line up to the (nominal) buckling stress f, and thereafter 
by a-line of which the slope (stress/strain) is 1/nth of the slope of the elastic line. Results such 
as those shown in Figs. 3 and 5 do not, of course, conform very closely to this simplified diagram ; 
but this type of diagram is obtained only from plates lying somewhat outside the range of normal 
practice in aircraft construction. 


After trial of several different expressions for f, and », the following were adopted :—. 
[ats t 
= OY 6 were : EY ee ee ‘ese one ee ee ee 
Io | i 00:21 3 f (2) 
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Although, in the choice of these formulae, some attention has been paid to theoretical con- 
siderations, the formulae should be regarded as founded entirely on the experimental results. 
They represent the behaviour of plates of the dimensions normally employed in aircraft 
construction with reasonable accuracy ; but outside that range they are not applicable. With 
regard to the formula for f, greater accuracy over a wider range could be obtained by substituting 
for (2) an expression similar to (1) with a lower coefficient of t/r (this adjustment might perhaps 
be regarded as an allowance for the effect of initial waves or waves formed in first loading) ; 
but the gain in accuracy in the practical range is considered to be outweighed by the greater 
complexity of the formula. 


The expression for 7 has an even less rational basis ; certainly increases with the ratio b/f, 
and the increase is greater, the smaller the value of 7. Otherwise the sole: justification for this 
expression is that it enables the experimental results to be reproduced with reasonable accuracy. 
The form of the expression for 1 suggests that ~ can be less than 2, a result at variance with the 
conclusions of reference 1. It will be found, however, that the values of b/t that make 1 < 2, 
also make f, > about 6 tons/sq. in., so that the exact slope of the line above f, makes very little 
difference to the total load at-reasonable values of the edge stress. 


The simplified diagrams obtained by using the values of /, and m given by (2) and (3) are 
shown on Figs. 3-14, and in Table 1 the experimental values of the mean stresses at buckling and 
at 5 tons/sq. in. and 10 tons/sq. in. edge stress are compared with the values given by the simplified 
diagram for all the plates tested. The mean stresses quoted in Table 1 are the mean stresses in 
the free width of plate, the load carried by the 2-in. widths of panel actually supported by the 
combs having been subtracted. The loads recorded on Figs. 3-14 are, of course, the total loads. 


The range over which the formulae (2) and (3) may be regarded as valid is rather difficult to 
define, since it includes both flat plates and plates for which the first term of (2) is negligible in 
comparison with the second. On the whole, it appears reasonable to limit the application of the 
formulae to plates for which the value of /, given by (2) is less than 5 tons/sq. in. It is improbable 
that any practical form of construction would have a higher value of f,, and Figs. 3-14 and 
Table 1 show that with this limitation the agreement between the estimated and the actual 
stresses 1s reasonably close. 
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Application of Results.—The empirica: formulae (2) and (3) for f, and. respectively are 
suggested for use in analysis on the lines indicated in reference 4. It may, however, be considered 
that the assumption that the load carried by the panel is proportional merely to f 


b ; sete 
(eq uivalent to the assumptions f/f, = 0,” « j ) is easier to apply and is equally well justified by the 


experimental results. In the last two columns of Table I are listed the ratios of the effective 
widths of the panels to their thicknesses at 10 tons/sq. in. edge stress according to both the 
theoretical and experimental results. According to the empirical formulae (2) and (3) there should 
be aslight increase in this ratio with increasing free width of panel ; the experimental evidence for 
this variation with width is ‘very slight. On the other hand the formulae (2) and (3) indicate 
a slight increase in the effective width with increasing curvature, and in this respect the experi- 
mental results agree better with the empirical formula than with the assumption that the free 
width is proportional merely to the thickness. 


So far as the values of effective width at 10 tons/sq. in. are concerned, therefore, there appears 
to be little to choose between the two assumptions. At lower values of the edge stress, the 
empirical formulae (2) and (3) will normally give better results. It may be suggested that this 
superiority is of little importance, on the grounds that the conditions at failure are alone of 
interest ; but it must be remembered that in many cases, e.g. tubes tested in bending, the 
limiting stress is developed only locally. In such cases, underestimate of the load carried at 
lower stresses wifl lead to underestimate of the collapsing load. 


Lastly, it is necessary to consider whether the method of reference 4 or the “ reduced effective 
width ’”’ method is easier to apply to practical structures. The chief difficulty of the second 
method consists in determining which parts of the structure are subjected to compression. Thus 
in estimating the strength of tubes in bending, it is necessary to find the neutral axis by a trial 
and error method. This difficulty is avoided by the method of reference 4. In fact, in the case 
of tubes of circular section if the values of /, and ” and of the ratio of the weight of the stringers 
to the weight of the sheet be given, the analysis of reference 4 enables the bending moment 
corresponding to any value of the limiting edge stress to be written down at once. Similar analysis 
could, of course, be worked out for other cross sections under other types of loading. Finally, it 
may be remarked that the method of reference 4 may be used in conjunction with the “ reduced 


; Sb SD 
effective width ’’ assumption merely by inserting the values f, = 0,1 = C ;where C is a constant. 
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AERODYNAMIC SYMBOLS 


1, GENERAL 
m Mass 
t Time 
Vv Resultant linear velocity 
Q Resultant angular velocity 
o Density, o relative density 
v Kinematic coefficient of viscosity | . 
R Reynolds number, R = /V/» (where / is a suitable linear dimension) 
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Normal temperature and pressure for aeronautical work are 15° C 
and 760 mm. 
For air under these ¢ p = 0-002378 slug/cu. ft. 
haan = OG: 11077 {t.*/sec, 

The slug is taken to be 32:2 Ib.-mass. a 

Angle of incidence 

Angle of downwash 

Area 

Span 

Chord 

Aspect ratio, A = 5?/S 

Lift, with coefficient C, = L/4eV?S 

Drag, with coefficient Cp = D/$pV?S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/$eV7bS 

Pitching moment, with coefficient C,, = M/4pV?cS 

Yawing moment, with coefficient C, = N/4}pV7dS 


2. AIRSCREWS 


SOrn WS oO:s 


Revolutions per second 

Diameter 

Power 

Thrust, with coefficient k, = T/on?D4 
Torque, with coefficient kg = Q/on?D*® 
Efficiency, 7 = TV/P = Jky/2ak, 
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Summary.—The paper deals with transition from laminar to turbulent flow in the boundary 
layer of a body caused by turbulence in the free stream. A non-dimensional number 
representative of the conditions in the boundary layer at transition is derived, and it is shown 


that this number should tend to have a definite value when the transition is caused solely by 
isotropic turbulence. This number, denoted by &, is 


TC OROmean| 


where v is the momentum thickness of the boundary layer, 


nA m@ is the root-mean-square value of the downstream component of the turbulent 
velocity u;, 


v is the kinematic viscosity, 


L is the Taylor scale of turbulence given by in R,dl, where R, is the correlation 
0 


coefficient for u-velocity fluctuations at two points separated by a lateral 
distance J, 


and = 4, is the mean tangential velocity at the outer edge of the boundary layer. 

Values of & are calculated from the results of experiments made on flat plates, an elliptic 
cylinder and a long streamline body of revolution in isotropic turbulent streams behind square- 
The average value of & obtained is 0-985, and of the six values calculated, four 


differ from the average value by amounts within + 6 per cent., one is 28 per cent. greater and the 
other is 29 per cent. smaller. 


Calculations of & are also made for two flat plates and thirteen aerofoils in wind-tunnel 
streams. The conditions of flow in these streams are diverse and. differ from those taken in 
the derivation of &. A representative value of ~, for these models appears to be 0.49. 


1. Theory.—It has been suggested by Professor Taylor!!® that the way in which turbulence 
is most likely to effect a laminar boundary layer is through the action of local pressure gradients 
which necessarily accompany turbulent flow, and that these gradients when sufficiently great 
and directed against the flow will cause a reversal of flow and a separation from the surface : and 


~ also that pressure gradients smaller than those necessary for the separation of a laminar boundary 


layer in a steady stream may be capable of causing the boundary layer to become turbulent. 


2 


The relation for the root-mean-square value of the fluctuating pressure gradients in the 


downstream direction wh (dp/dx)?, obtained by Professor Taylor! for isotropic turbulence 
produced by passing a stream through a rectangular grid of Mesh M is 


ae EL 
vEneit (Malt! 
v 


where 
M,/ ie 
> 60, 
Vv 
d is a length which determines the scale of the eddies responsible for dissipation of 
energy by viscous forces, 
An is a length which can be determined from diffusion experiments, 


wu; is the root-mean-square value of the downstream component of turbulent 


velocity 1, 
££ -j—1 
‘(75) 
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U, is the mean velocity in the free stream 


Aetak | ws is the downstream distance behind the grid. 


Professor Taylor then proceeds to consider the effect of wh (0p/0x)? on separation for the laminar 
boundary layer of a flat plate and that for a sphere, and, on the assumption that A/A, and A 


are constants, obtains a relation for the condition that the boundary layer on a flat plate will 
become turbulent, and the well-known relation 


ma! AE) 


for the critical Reynolds number of a sphere. 


2. In a recent paper® the writer and Dr. Preston adopted the idea that transition can be caused 
by a separation due to the adverse pressure gradients associated with turbulence in the free 
stream, and obtained a non-dimensicnal number whose value at the point of transition in the 
boundary layer of a body of any form in an isotropic turbulent stream should tend to be constant, 
when the transition is caused solely by turbulence in the free stream. According to the KArman- 
Pohlhausen theory it is to be expected that the instantaneous distribution of velocity in a 


cross-section of the boundary layer would depend on the local value of — = sg and that the 
1 
S Op\- ae 
non-dimensional number — Te (& should tend to have a definite value at transition, 
1 s 


where § is the boundary layer thickness, w, is the mean tangential velocity at the outer edge of 


3 


the layer, vis the viscosity, p is the instantaneous pressure, and s is the distance along the surface. 
On substitution of the relation 


(Zin _ V2dpu? al 
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Se [rep 
the number — — &) becomes 
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which was written, in the earlier paper?, in the form 


pence) (LE) elt 


where uw, is the maximum value of u,, 


bol 


ye ra, 


vis the kinematic viscosity, 


vis the momentum thickness, 


given by i Ca ( 1 — a) dy, 


and « _ isthe mean tangential component of the velocity in the boundary layer at a normal 
distance y from the surface. 


The number was then reduced to the simpler form 


) Gp GF 


obtained by taking the (2)'" root, and ignoring the factor 
@as (2) (va U; ie 5 
Ady Uv Uy 


It was shown that values of (vw-;/v) (v/M)# (v/w,v)? at transition, calculated from measure- 
ments taken for a long streamline body of revolution in water streams behind square-mesh 
grids and also for a flat plate in isotropic turbulent streams, from measurements made by Hall 
and Hislop®, had the same order of magnitude. 


which was taken to be constant. 


3. It is convenient in the present paper to introduce the Taylor scale of turbulence, L, 
wheres) { * R,dl, and R, is the correlation coefficient for u-velocity fluctuations at two 


0 
points separated by a lateral distance 7. The number whose value should, on the assumptions 
made, be constant at transition can then be written in the form 


eta eG | 
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The condition that transition is caused solely by turbulence holds when the value of A, that 
is (8?/v) du,/ds, is zero. The value of 8/v calculated from a solution of the laminar boundary 
layer equations for a known pressure distribution is therefore constant. The value given by the 
Karman-Pohlhausen theory is 8:51. A representative relation for L/M obtained by Dryden 
and others (Table III) from measurements made in the turbulent streams behind square-mesh 
grids, M = 0-25, 0:5, 1-0 and 3-25 in., is L/M = 0: 147 + 0-0025 (x/M). The relation for A 
obtained by Taylor !” (p. 441) is 


Uh a aye 
( ) dx MA?’ 
which can be written 
a) y 5 
a (x/M) - APT 


Measurements of U,/ wh 12 obtained from the experiments of Dryden and others!4 (Table IV), are 
plotted against x/M in Fig. 1. Included in the figure is a curve given by the relation 


The measurements lie close to this curve so that A = 5 [0-147 + 0-0025 x/M;*?. The value of 
(L/M)'/A is therefore 0-2. Experimental values of 4,/A obtained for turbulent flows behind 
square-mesh grids are 0:5 (Taylor!’, p. 478) and 1-6 (Hall®, p. 13), and theoretical values obtained 
for definite eddy systems are 0-7 (Taylor?, p. 13) and 0-8 to 1 (Green, p. 12). Anothef experi- 
mental value, which lies between Green’s theoretical values, 0-8, and 1-0, has been obtained by 
Hall, since the publication of his paper?. It is a matter of some difficulty to obtain a good value 
for Ay/A, and no experiments have been made to determine its dependence on the turbulence 


number L Al a Taking a representative value A,/A = 1, and substituting (8/v)? = 72-5 and 


(L/M)#/A = o we get 
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The analysis therefore leads to the conclusion that the value of 


vas [(X2) @ GI 


denoted by @, (see footnote), should tend to have a constant value at the point of transition 
in the boundary layer of a body in a turbulent stream behind a square-mesh grid, when the 
transition is caused solely by turbulence. 


1 2 
The symbol & was used for the non-dimensional number ( (i) (en in the earlier paper?. 
f 1 


4. The Taylor relation 
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was obtained by the substitution of 
rN V 
biel (ey. set 
M M ty, uy | 
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(see Ref. 17, p. 440). From Dryden’s experiments we have obtained the relation A = 5(L/M), 
and on substitution in the relation 


bole 


in the relation 


we obtain 


A 3 

oder a 

_A relation of this form can be obtained from considerations of similarity (see Taylor, Ref. 18, 
p. 299). For similar flow systems the Reynolds stresses are proportional to 9 u? (uw, = typical 
velocity), and the rate of dissipation per unit volume is proportional to o ( ed ue) / Loti atypical 


length). For isotropic turbulence the rate of dissipation due to viscosity is 15 ww?/22. The two 
ways of regarding dissipation are consistent if 
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where Bisa constant. It should be noted that the values of U,/ 1: uz measured by Dryden at a 
given value of x/M behind a grid were independent of U, but they were not independent of M. 
The conditions for similarity were therefore not fully satisfied. The relation A = 5(L/M)! given 
in section 3 was obtained for mean curves of U,/ Nie u? and L/M against x/M, from Dryden’s results 
for M = 0-25, 0-5, 1-0 and 3-25 in. 


5. Values of & at transition (A = 0) calculated for a large number of models of different 


size and shape and tested in turbulent streams for which the values of A ue and of L differ 
appreciably are given in Tables 1 and 2. The experiments analysed fall into two groups: 
I, those made on models in streams behind square-mesh grids, and II, those made on models in 
wind-tunnel streams. The conditions of the experiments of Group IJ do not satisfy those 
assumed in the derivation of &, but the calculations were considered to be worth while since 
they show the probable variation in & to be expected and whether a reasonably representative 
value of &€ can be obtained. 
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The results of twenty experiments have been analysed, five in Group I and fifteen in Group II. 
Descriptions of the models tested and relevant details of the experiments made on them are given 
below. 


Group I. Models in Streams Behind Square-Mesh Grids. (1) Two Flat Plates. Cambridge?.— 
One plate was of smooth plate glass, three-eighths of an inch in thickness, and had an elliptical 
leading edge. The other plate was of rolled duralumin of thickness one-eighth of an inch, and 
had a sharp leading edge. Both plates had excellent surfaces. M =—0-5 and 1 in. Position of 
transition point obtained from surface-tube observations. 


(2) Flat Plate. National Bureau of Standards, U.S.A.4.—Polished aluminium plate, 0-125 in. 
thick, with sharp bevelled leading edge. M = 0:5 and 1in. Position of transition point obtained 
from inspection of oscillograph records. 


(3) Elliptic Cylinder. National Bureau of Standards, U.S.A.°.—The cylinder was made of 
wood, and had a smooth polished surface. Major and minor axes 11-78 and 3-98 in., respec- 
tively. M = Lin. Position of transition point deduced from change in shape of velocity distribution 
in boundary layer. 


(4) Long Streamline Body of Revolution. Small Water Tunnel, N.P.L.?.—Brass model with 
ellipsoidal nose. Polished surface. M = 0-5 in. Position of transition point obtained from 
behaviour of filament bands in boundary layer. 


(5) Long Streamline Body of Revolution. Water Tunnel, N.P.L.2,—Same as (4), except 
Me. -2o%in. 
Group II. Models in Wind-Tunnel Streams.—(6) Flat Plate. National Bureau of Standards, 


U.S.A.4.— Polished aluminium plate, 0-125 in. thick, with sharp bevelled leading edge. Position 
-of transition point obtained from an inspection of oscillograph records. 


(7) Flat Plate. Compressed Aty Tunnel®.—Aluminium plate 0-125 in. thick having an elliptical 
leading edge. 


(8) Symmetrical Aerofoil. D.V.L. 5 x 7 Metre Tunnel’.—Composite model with a marble 
surface finish. Chord = 10-5 ft. Maximum thickness = 0: 127 chord, 0° incidence. 


(9) Symmetrical Aerofoil. D.V.L. 5 X 7 Metre Tumnnel’—Same as (8) except maximum 
thickness = 0-206 chord. 


(10) Aerofoil N.A.C.A.23012. Variable Density Tunnel®.—Duralumin model. Chord = 5 in. 
Maximum thickness = 0-12 chord. Incidence for minimum profile drag. 


(11) Aerofoil N.A.C.A.0012. Variable Density Tunnel®.—Duralumin model. Symmetrical 
section. Chord = 5in. Maximum thickness = 0-12 chord. 0° incidence. 


(12) Aerofou N.A.C.A.0012. Full Scale Wind Tunnel N.A.C.A.°.—Model of steel-spar 
construction with ribs spaced at 12 in. intervals. The covering was + in. aluminium sheet. 
Glossy waxlike finish. Chord 6 ft. Maximum thickness = 0-12 chord. 


(13) Aerofoil N.A.C.A.0012. Compressed Airy Tunnel©,—Steel model chromium plated and 
hand-finished surface. Chord = 8 in. Maximum thickness = 0:12 chord. 0° incidence. 


(14) Jowkowskt Aerofoil. Compressed Air Tunnel4.—The model was made of cast aluminium 
and had an aerodynamic smooth surface up to Reynolds number 15 X 10°. Chord = 2 ft. 
Maximum thickness = 0-127 chord. Symmetrical section 0° incidence. Position of the point 
of transition determined by “‘ wires ”’ method. 


(15) Jowkowski Aerofoil. Compressed Air Tunnel4.—Same as (14) except maximum thickness 
= 0-236 chord. 


(16) Aerofoil A. Compressed Air Tunnel/®—Hiduminium model, hand finished. Symmetrical 
section. Chord = 7-84 in. Maximum thickness = 0-12 chord. 
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(17) Aerofoil B. Compressed Aiy Tunnel!®.—Hiduminium model, with good buffed surface. 
Chord = 8in. Maximum thickness = 0:12 chord. Incidence for which profile drag is a 
minimum. 


(18) Joukowski Aerofoil. No. 3, 7 ft. Tunnel, N.P.L.?.—Gunmetal model having a smooth 
polished surface. Chord = 39-7 in. Maximum thickness = 0:1507 chord. Transition point 
and transition region determined from surface-tube readings. 


(19) Aerofoiul C. 13 ft. x 9 ft. Tunnel, N.P.L.°.—Wooden model with hand-polished varnished 
surface. Chord = 70 in. Maximum thickness = 0-12 chord. Symmetrical section. 0° 
incidence. 


(20) Aerofoiul C. 13 ft. x 9 ft. Tunnel, N.P.L.°.—Same as (19), except chord = 42 in. 


6. For models, Nos. 1-6, 14, 15 and 18, the transition point was located experimentally, and 
for each of the other models (aerofoils) the position taken was that for which the profile drag 
calculated by the Squire-Young method!® for the known pressure distribution is equal to the 
measured drag. The experimental conditions considered are those for transition at or near 
the point on the surface where A = 0. Calculations of vw,/v and A were made bythe Karman- 
Pohlhausen method of solution of the laminar boundary layer equations for distributions of 
u,/U, determined either from measured or theoretical pressure coefficients, (f — f)/$eU5 (mean 
values), by the standard relation 


bade a E ima (p sank 
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Measured values of gf im and L were taken in the calculation of & for models 1, 2 and 3, and 


measured values of hs PH (taken to be (1/3-3)" of the value of w; measured with a fluid-motion 
microscope) and calculated values of L (from relations in Table 3, Ref. 14) for model 4 (and 5). 


ie £ 
For each of the other models a value of she GC) for the free stream, calculated from the 
0 


U 


for the known sphere number, Rorit. was taken. 


of 3 
relation log,, co CE @ | = 0-991 — 2-193 log.) [Rar X 1075], given in Ref. 15 
0 


7. Discussion of Results—The results of the calculations made for the models in turbulent 
streams behind square-mesh grids (Group I) are given in Table 1. Values of & are plotted 


L,/@ 
Vv 


the six values differ from this average value by amounts within +6 per cent., whilst the value 

for the flat plate (2a) is 28 per cent. greater and that for the long streamline body (5) is 29 per cent. 

smaller. The values of & for the flat plate (1) change from 1-07 to 0-90 of its mean value 

over the range vu,/v = 293 to 559 covered (Table 1). The values of & for the elliptic cylinder 

and the long streamline body do not change appreciably with A, over the small ranges covered 
M,/ ui | 

> 60, assumed in the derivation of the 


against the turbulence number in Fig. 2. The average value of & is 0-985. Four of 


(Table 1). It will be noticed that the condition 


relation for &,, is satisfied for each model. 
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TABLE 1—GROUP I 


as = Vit 1\3 
Model and Number i A a LV ia MVia E U, |) 
; oes : Lin ft. 
(1) Flat Plate 559 0 Od) [Soc 1 99 | 0-915 0-021 
(Cambridge) 493 0 0-5 | 32-6 98 | 0-980 0-025 
437 0 0:5 | 33°9 109 | 0-980 0-026 
360 0 O'S | 3l-4 435*7 111 1-015 +1-012 | 0-032 
328 0 1-0] 33-9 118 | 1-051 0-031 
322 0 1-0 | 41-7 148 | 1-060 0-032 
293 0 1-0 | 41-0 148 | 1-085 0-035 
(2) Flat Plate ee oe yen 0 0:5 186 518 | 1-260 0-069 
National Bureau of Stan- |(b) 470 0 1-0 115 338 | 0-965 0-027 
dards 
(3) Elliptic Cylinder es 217 0-4 | 1-0 348 1160 | 1-050 0-077 
National Bureau .of Stan- 235 0-1 | 1-0 327 1090 | 1-065 0-073 
dards 242 0 Lt) 315 1050 1-040 0-070 
205 —0-15) 1-0 294 980 | 1-020 0-065 
(4) Long Streamline Body of 210 32 Geto 0-045 bo-053 
Revolution. Water | 146 fl] 98/95 | 54 748 179 | 9-890 9-99 | 0-061 
Tunnel (N.P.L.) 168 0-3 0-5 47-1 196 | 0-920 0-054 
Estimated from values given 173 0 0:5 45-0 187 | 0-930 0-054 
above 
(5) Long Streamline Body of 310 0:6 | 0:25 28-0 76 | 0-700 0-028 
Revolution 260 0-3 | 0°25 34°5 108 | 0-700 0-032 
Estimated from values given 285 0 0-25 31-3 91 | 0-700 0-030 
above 


The results for the models in wind-tunnel streams, Group II, are given in Table 2. Values 
of & for these models and those for Group I are plotted against w,v/v in Fig. 3. The values for 
Group II are lower than those for Group I. A representative value of & for the models of 
Group II, except models 8, 9, 10, 11 and 12, see later, is 0-49. 


It is hardly to be expected that the values of &, for the models of Group II would lie closely 
on a smooth curve for they were obtained from experiments made under widely different 
conditions. It is necessary therefore to consider broadly whether transition on the models (both 
groups) is affected by conditions other than those due to turbulence, and whether the conditions 
in the turbulent streams differ from those taken in the derivation of the relation for &. The 


exactness of location of the transition points and the reliability of the values of 1), v, aif u? and L 
taken in the calculation of & must also be considered. 


8. Conditions other than those due to turbulence which affect transition are a gradient in 1“, 
2 

or more explicitly A (= 2 at} roughness due to small-scale excrescences, and surface waviness. 
Transition is assisted if #, falls in the direction of flow for, as is well known, transition may 
occur in a stream substantially free from turbulence when A reaches a limiting value, which 
according to Howarth!* is about —7:5. The value of &, at transition should therefore be greater 
for a positive than for a negative value ofA. The values of & plotted in Fig. 3 are for the 
condition A = 0, that is the condition taken in derivation of the relation for &,, except those 
for models 17, 19 and 20, for which A= 0-3, 0-6 and 0:6 respectively. Table 2 shows that for 
these models & changes slowly with A in the region near A = 0. 
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TABLE 2—GROUP II 


Sphere Ve (13 
Model and Number A " er D ae @ EL 
fori ey belie b. infti 

(6) Flat Plate. National Bureau of Standards. . b. 697 0 2-70 5 0-0132 0-560 
(7) Flat'Plate. (C.A.T.) ~. aye rie =e sr 448 0 2-24 eth 0-0173 0-614 
8) Symmetrical Aerofoil. (D.V.L.)  .. ae ne 1230 0 3°51* <8 0-0072 0-335 
(9) Symmetrical Aerofoil. (D.V.L.) ee ae ~ 1300 0 3-51* eon 0-0072 0-305 
Omelet N Ale Aloa0 12M. Dry Mune ae al E90 NT 0 15 | 8 | 0-0495 | 0-736 
(11) Acrofoil NA. CALOOIZ 9( V-Day 3. as ae 700 0 1-5 aie 0-0425 0-796 
(12) Aerofoil N.A.C.A.0012 (F.S.T.) - a .. | 1340 0 3°3 eed 0-0076 0-315 
(13) Aerofoil N.A.C.A.0012 (C.A.T.) y: am ae 780 0 2-24 aa 0-0173 0-371 
(14) Joukowski Aerofoil. (C.A.T.) 7 A? - 980 0 2°24 a 0-0173 0-512 
(15) Joukowski Aerofoil. (C.A.T.) Bes oe eek 00S. 0 2°24 Ge 0-0173 0-462 
(16) Aerofoil A. (C.A.T.) .. es Pe 7. Re 835 0 2°24 a 0-0173 0-468 
(17) Aerofoil B. (C.A:T.)  -. is ae ate ah ee ae 2-94 6 00173 mee 
(18) fone: Aerofoil. (N.P.L.) AP “e ee 370 0 2°10 Riga 0-0222 0-491 
(19) Aerofoil C. Neb eye C == 0. | S. ss Bar aia? 0-6? 8 13-55 ms 09-0092 | 0°465 
1205 1-4 3°25 3 0-488 

(20) Rei cle ave.) C = 42 La = 7: 7 955 0:6 5 3°55 Bey 0-0092 0-388 
co 1-4 3°25 | 3 | 0-412 


* Taken to be 0-91 the atmospheric value 3-85’. 


Fourteen of the models were made of metal, and of the others, models 3, 19 and 20 were made of 
wood, one of the Cambridge flat plates of glass and models 8 and 9 of marble (see section 5), 
All the models had polished surfaces, and they can be taken to be aerodynamically smooth 
at the Reynolds numbers of test. Surface waviness is also unlikely to affect the results. 


9. The reliability of a calculated value of & depends primarily on the reliability of the values 


of u,, v and VA wu? taken : and the reliability of wu, and v depends largely on the exactness of 


location of the transition point. The positions of the transition points for the flat plate 1 and 
the Joukowski Aerofoil 18 were determined by the surface-tube method, for the flat plate 2 
from inspection of oscillograph records, for the Joukowski Aerofoils 14 and 15 by the wires method, 
for the body of revolution 4 and 5 from the behaviour of filament bands in the boundary layer, 
and for the elliptic cylinder from the change in the shape of the velocity distribution in the 
boundary layer (see section 5). For the other models (Aerofoils) the transition points were taken 
to be those for which the profile drag calculated by the Squire-Young method!® for the known 
pressure distribution is equal to the measured drag. The relative precision of location obtained 
by these different methods when applied to bodies of different shape in different turbulent 
streams is not known: but for the same flow conditions the filament method gives results 
consistent with those given by the surface-tube method, and the wires method applied to 


i0 


Joukowski Aerofoils 14 and 15 in the C.A.T. gives results consistent with the Squire- Young 
method. For an abrupt transition, such as that which occurs on an aeroplane wing im flight, 
the position of the transition point obtained by the Squire-Young method for a thin wing 
at a low incidence is in good agreement with that measured by the surface-tube method. 
Transition on a body in a wind-tunnel stream having a high degree of turbulence is not 
however an abrupt process, and the position of the transition point predicted by the Squire- 
Young method occurs beyond the position given by a surface tube. The values of v and & 
will therefore be over-estimated. Thus, the value of & obtained for the measured point of 
transition on Joukowski Aerofoil 18 in the stream of the N.P.L. No. 3, 7ft. Tunnel, for which 


mr: 5 ae 
o/ a G) = ():0222, is 0-491, whereas the value obtained for the position given by the 
Uy 


Squire-Young method is about 0-590. No reliance can therefore be placed on the values of @ 


caw iL 
obtained for Aerofoils 10 and 11 in the very turbulent stream of the V.D.T. Vii () = 0-0425, 
0 
and they are not plotted in Fig. 3. The values of & for the Joukowski Aerofoils 14, 15 and 
18 for which the transition points were measured, are 0-512, 0-462 and 0-491 respectively, and 
the mean of these values is 0-49. The line & = 0-49 is drawn in Fig. 3, and it is seen that 
it represents the trend of the values of & obtained for the Aerofoils 13, 14, 15, 16, 17 and the 
flat plate 7 in the C.A.T., the flat plate 6 in a Bureau of Standards Tunnel, and the Aerofoils 


— 1 
19 and 20 in the N.P.L. 13 ft. x 9 ft. Tunnel. The values of vi ei for these tunnels are 
0 


0:0173, 0-0132 and 0-0092, respectively. The values of & for Aerofoils 8 and 9 in the 
D.V.L. 5 <* 7 metre Tunnel and Aerofoil 12 in the Full-Scale Wind Tunnel are however low. 
(See below). 


= iz 
10. As mentioned in section 6, the values of ao G) * taken in the calculation of &, for 
0 


the models in the wind-tunnel streams, Group II, were obtained from the relation 


um (Ds 
log io 100 Ni bie (2) = 0-991 — 2-193 log 4) [Rai X 10-°]. This relation was obtained 


0 
from tests on spheres in the turbulent streams behind square-mesh grids made for the range 


Rat, = 1:1 x 10° to 2-85 x 105, the corresponding range of NA — )° for a 6 in. sphere bein 
p § § U ie Pp § 
0 


0-092 to 0-011. The curve given by this relation represents the trend of the observations over 
this range (Fig. 4). It is known, however, that a sphere ceases to be a sensitive turbulence 
indicator as the limiting value R,., = 3-85 x 105 is approached, and there is no evidence to 
show how closely the above relation then holds. The values of R,;, for the N.A.C.A. Full-Scale 
Tunnel and the D.V.L.5 x 7 metre Tunnel, 3-3 x 105 and 3-51 x 105 respectively, are well 


——— aE 
5 
beyond the values covered in the sphere tests (see Fig. 4) and the values of J ad G) 


0 
calculated for these tunnels from the above relation may be too low and if so this would 


account for the lowness of the values of & for Aerofoils 8,9 and 12. It is apparent that a 


—- z 
sphere test in a wind-tunnel stream gives the value of e Ly (+) ; for a turbulent stream 
0 
behind a rectangular grid which has an equivalent effect on the boundary layer of a 
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sphere : but it does not follow that the two streams would have the same effect on other bodies. 
Further, casual changes in the direction and speed of a wind-tunnel stream may have a greater 


ay: 1 
disturbing effect on some models than on others at a low value of ve er 
5 NS 


11. The conditions in a boundary layer at transition depend on the local turbulence in the 
free stream and on the cumulative effect of turbulence on the flow in the boundary layer up to 
the transition point. This cumulative effect depends on the shape and size of a body, on the- 
mean velocity of the stream and on the turbulence in the free stream, and causes a departure 
from the conditions taken for a purely laminar layer, and, in particular, an increase in thickness. 
The measurements made on the models do not allow the actual thickness at the point of 
transition to be taken in the calculation of €,-but if this had been possible the consistency in the 
values of & might have been greater. 


12. Apart from the N.P.L. No. 3-7 ft. Tunnel, the tunnels in which the experiments on the 
models of Group II were made are of the return circuit type, with a honeycomb, a setting chamber 
and a contracting entrance nozzle, through which the air flows on its way to the working chamber. 
The turbulence spectrum for each of these streams differs from that for a turbulent stream 
behind a grid, and the conditions taken in the derivation of & are not likely to hold as closely for 
the models of Group II as for those of Group I. It is therefore not surprising that the repre- 
sentative value of & for the models in the wind-tunnel streams differs from that for the models 
in the turbulent streams behind grids: and in view of the diversity of the conditions in the 
wind-tunnel streams the agreement between the values of & obtained for the models in them is 
as close as would be expected. 


13. Finally, the analysis leads to the conclusion that when transition is caused solely by 


turbulence in the free stream, the value of 3-35 ( vn/ By ee ey denoted by &, at 
Sort il U4v 


the transition point in the boundary layer of a model in a turbulent stream behind a square- 
mesh grid is 0-985. A representative value of & for a model in a wind-tunnel stream is 0-49. 


The corresponding values of (s/s) en: e My are 0: 294 and 0-146. 
v if, Uv 
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SYSTEM OF AXES 


‘Symbol 


| Pee ee : 

Axes | eee | to fertdinal lateral normal 

: forward starboard downward 
direction | 

RORCE* Say abol ae ¥ Z 
Moment pa neool ye ee N 

Designation rolling pitching yawing 

Angle of 

Rotation shat Y y 
Velocit Linear u v w 
y Angular p g | r 
Moment of Inertia | cay: B 3a 


Components of linear velocity and force are positive in the positive diection 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 

The symbols for the control angles are :— 

€ aileron angle 

n elevator angle 
ny tail setting angle 
¢ rudder angle 
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AERODYNAMIC SYMBOLS 


1. GENERAL 


WO OS SS 


ZeHX OES oH 


Mass 

Time 

Resultant linear velocity 

Resultant angular velocity 

Density, o relative density 

Kinematic coefficient of viscosity 

Reynolds number, R = /V/» (where / is a suitable linear dimension) 


Normal temperature and pressure for aeronautical work are 15° C 


and 760 mm. 

For air under these ( e = 0:002378 slug/cu. it. 
conditions | » = 1:56 x 10-* ft.?/sec. 

The slug is taken to be 32-2 lb.-mass. 

Angle of incidence 

Angle of downwash 

Area : 

Span . / 

Chord 

Aspect ratio, A == 6?/S 

Lift, with coefficient C, = L/ZeV*S 

Drag, with coefficient Cy) = D/#eV?S 

Gliding angle, tan y = D/L 

Rolling moment, with coefficient C, = L/}eV?0S 

Pitching moment, with coefficient C,, = M/4eV?cS 

Yawing moment, with coefficient C, = N/3 VS 


2. AIRSCREWS 


S710 Ss 


Revolutions per second 

Diameter 

V/nD 

Power 

Thrust, with coefficient ky = T/on?D* 
Torque, with coefficient kg = QO/on*D° 
Efficiency, 7 == (PV /P = Jer/2ako 


ho 
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Summary.—Reasons for enquiry.—Tests were required to check the validity of the theoretical 
corrections for wind tunnel interference in the 5-ft. open tunnel when the wing is of large span 
or spans the jet. 

Range of investigation.—Two wings of 8-in. and I-ft. chord were tested on a horizontal diametric 
axis :— 

(a) Spanning the 5-ft. open jet wind tunnel. 
(6) Cut down successively to 0-9, 0-8 and 0-6 of the tunnel diameter. 
Lift, drag and pitching moment were measured in all cases. 


The following symbols are used in the text :— 


Wing area S 
Wing chord .. c 
Tunnel diameter .. Me a D 
Tunnel cross sectional area aH AP in, aD 


Conclusions.(a) Wings inside the jet diameter.—For values of S/C < 0-25 and c/D < 0-2 the 
standard corrections may be applied. Beyond this limit further systematic tests are required 
to derive satisfactory empirical corrections. Pending such tests it is better to calculate the slope 
of the lift curve in free air for the appropriate aspect ratio and use this to deduce the incidence 
correction. 


(b) Wings spanning the jet.—The corrections given by Glauert apply with a modified value of 
N = 1-6 forc/D < 0-2. Further tests are required for larger chord/diameter ratios. 


1. Introduction.—Theoretical estimates of the boundary interference on the characteristics of 
wings in open jet tunnels are based on simplifying assumptions which may not be justified in 
tests on comparatively large models. Tests have, therefore, been made in the 5-ft. open jet 
tunnel at the R.A.E. on two wings to check the validity of the theoretical corrections in the 
following cases :— ; 

(1) Wing chord/jet diameter + 0-2, span less than tunnel diameter. 


(2) Wing chord/jet diameter + 0-2, wing spanning the jet. 


' The wind tunnel corrections on these large models are complicated by “ shadow effect ’’! and the 


difficulty of measuring wind speed. In the tests in this report the effect of model on wind speed 
has been investigated (Appendix). 


*R.A.E. Report, February, 1941. 
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Fig. 2. Profile drag coefficients. 
R.A.F.34 (8 in. chord) and R.A.F.48 (1 ft. chord). 
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3.22 Profile drag (Fig. 2).—The profile drag curves for the 4:5 ft. and 4 ft. spans are identical, 
while the profile drag for the 3 ft. span is uniformly greater by about 6 per cent. for both aerofoils. 
This is probably an end effect!* 14 since it appears to be a uniform increment for both wings. 
For the case of the wing spanning the jet, the induced drag is overestimated. Estimated and 
empirical values are compared below. 


divs es oe 


Estimated N | Empirical N 


Condition | (Equation 9) 


Wing spanning jet. R.A.F.34 . a ie b 
R.A.F.48 4 1: 


3.23. Pitching moments (Fig. 3).—-The pitching moment curves shew good agreement for 
wings inside the jet boundary. For the wing spanning the jet there is a small increment of 
ACy = + 0-004 at C, = 0-6 for R.A.F.48. This difference may be more important on a 
thicker wing section. 


0-02 =o pease an eke oon 
[— 
o 3SPAN 
x o' SPAN 
Gy, + 45' SPAN 
aC & WING SPANNING JET. 


Fig. 3. Pitching moment coefficients. 
R.A.F.34 (8 in. chord) and R.A.F.48 (1 ft. chord). 


4, Discussion of results —The present tests shew that for wings of spans inside the jet diameter 
the correction given by equations (4) and (5) gives consistent lift and drag results for wing of 
varying span and chord for which S/C + 0-25 and the span is not greater than 0-9 x tunnel 
diameter. These conclusions were arrived at by Prandtl'®, Knight and Harris!*, Tani and 
Taimi‘* for span/diameter ratios up to 0-8. For profile drag measurements the effect of square 
ends on wings of small aspect ratio may be misleading, but further tests would be necessary 
before any recommendation, such as rounding the tips, could be made. 

For wings spanning the jet the value of M from equation (8) is correct, but the value of N is 
overestimated. 

There is evidence!’ that for wings of chords larger than 0-2D both spanning the jet and with 
spans inside the jet an additional correction must be applied. This correction is large and is very 
probably due to the effect of the finite dimensions of the working section. No satisfactory 
theoretical solution has been found and in the absence of any but scattered data it is hoped to 
make tests similar to these in this report on wings of 2 ft. and 1-5 ft. chord when an opportunity 
occurs. 

(59862) A* 
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For correcting the lift slope in 5 ft. tunnel tests on large chord wings it has been customary 
to estimate the value of 0C,/é« for the wing from other data and correct the measured values 
of « by an increment of the form — hC, to give the free air slope. Pending the systematic tests 
proposed above no better method can be suggested. 


Tunnel corrections on a complete model of span greater than the jet diameter are given in a 
further report. 


APPENDIX 


Experiments to determine the best position of the “‘ hole in the side’”’ of the tunnel with a large wing 
in the tunnel 


It has been customary in the 5 ft. tunnel to keep the speed steady by an automatic control 
working on the pressure difference between hole A, near the top of the tunnel side and atmosphere 
(see Fig. 4), whereas the calibration of the tunnel was made against three inter-connected holes 
in the roof. A wing at high lift was found to alter the pressure difference between these holes 
and atmosphere ; vertical traverses were, therefore, made to measure static pressure in the 
maximum section in the plane of the holes, and velocity before the model. This latter was measured 
well up the nozzle so as to avoid the velocity field of the wing itself. In Fig. 5 the velocity 
traverses are shown ; and also the static pressure in the maximum section divided by $oV? in 
the centre of the jet for two lift conditions of the wing and the tunnel without wing. It is seen 
that with the original hole in the side for the automatic (i.e. A,), the pressure is increased by 
1-3 per cent. by the presence of a wing at no lift, and by a further 1-7 per cent. when the wing 
has a lift of 62 lb. at 120 ft./sec. This corresponds to a reduction of speed of 1-5 per cent. in the 
jet when the pressure is kept constant by this hole. If, however, a lower hole A, is chosen (Fig. 4) 
a much more nearly constant result is obtained. A second hole was, therefore, fitted near A,, 
and these two holes are now used in place of the original roof holes and A. 
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TABLE 4 
R.A.F.34, (8 in. chord), } chord point 15 in. from the nozzle 


Measured values Corrected to A.R.% 
Condition Wind tunnel 
correction 

a? Wy ess Cm, a9 be | Can 
Wing spanning | —4 | —0-180 0:0149 | —0-0131 Corrected by —3-10 | —0-180 0-0124 
the jet —3 —0:118 0:0130 —0:0128 Glauert and —2 41 —0:118 0:0119 
—2Z | —0-063 0-0118. — Squire formulae | —1-69 | —0-063 0-0118 
—1 0-010 0-0106 —0-0084 equations (6)—(9) | —1-05 0-010 0-0106 
0 0-067 0-0100 —0-0058 —0-33 0-067 0-0100 
1 0-125 0-0101 —0-0021 0:38 0-125 0-0089 
2 0-192 0:-0118 — 1-04 0-192 0: 0086 
34 0-257 0-0138 0-0012 yp) 0-257 0-0088 
5 0-413 0-0209 —0-0012 3°00 0-413 0-0080 
7 0-573 0-0327 —0:0056 4-22 0-573 0:0079 

8 0-625 — — 4-97 0-625 — 
9 0-706 0:-0455 —0:0050 5:58 0-706 0-0080 
10 0-766 0-0524 —0-0031 6-29 0-766 0-0083 
Wing of 4-5 ft. —4 —0-190 0-0145 —0:0101 Equation —3-°14 —0:190 0-0127 
span. —3 —0-126 0-0125 —0-0094 (4) and (5) —2:44 —0:126 0:-0117 
ACR = 675 —] 0-010 0-0104 —0:0072 —1-04 0-010 0-0104 
0 0-069 0-0095 —0:0058 ==032 0-069 0-0093 
1 0-122 0-0096 —0:0021 0:45 0-122 0-0089 
2 0-197 0:-0108 —0:0001 1:12 0-197 0-0089 
3 0-271 0-0127 --0-0007 1:79 0-271 0-0090 
5 0-426 0-0182 —0-0029 3-09 0-426 0-0091 
f] 0-584 0 -0267 —0-0078 4-39 0-584 0-0097 
Od Oe Zs 0-0368 —(-0089 5:76 0-723 0-0107 
10 0-781 0:0424 —(0-0175 6-50 0-781 0:-0119 
Wing of 4 ft. =o —0-192 0-0150 —0-0100 Equation —3:10 —0:192 0-0129 
span =3 —0-122 0-0126 —0 0092, (4) and (5) —2 +43 —(122 0-0118 
OR 6 | —0-003 0-0108 —()-0069 —0:99 —0-:003 0-0108 
0 0-061 0-0098 —0:0050 —0-29 0-061 0-0096 
1 0-122 0-0096 —0-0017 0-43 0-122 0-0087 
2 0-192 0-0110 0-0008 1-10 0-192 0-0084 
3 0-258 0-0129 0-0010 1-80 0-258 0-0091 
os 0-413 0-0189 —0-0016 3:09 0-413 0-0094 
7 0-569 00-0282 —0:0078 4-35 0-569 0-0101 
9 0-718 0-0397 —0:0088 . 5:66 0-718 0-0107 
10 0-783 0-0457 —0-0083 6-36 0-783 0-0113 
Wing of 3 ft. —4 —0:179 0:-0150 —0:0090 | Equation —3-:05 —0:179 0-0127 
span —3 (0-121 0-0132 —0:0088 (4) and (5) | —2-35 —0=121 0-0121 
A-Rue 425 =i 0-006 0-0112 —0:0066 1-03 0-006 0-0112 
0 0-059 0-0104 —0-0061 —0:31 0-059 0-0101 
i 0-127 0-0107 —0-0006 0-33 0-127 0-0095 
2 0-182 0-0119 +0:0017 1:02 0-182 0-0095 
3 0-246 0:-0142 0-0020 1-69 0-246 0-0098 
5 0-397 0-0216 —0:0004 2°87 0-397 0-0101 
7 0-544 0-0323 —0-0073 4-09 0-544 00-0107 
9 0-684 0-0454 —0:0125 5°35 0-684 0:-0113 
10 0-748 0-0531 —0-0083 6-00 0-748 0-0122 


R.A.F.A8, (12 1n. chord), 


Measured values 


TABLES 


+ chord point 18 in. from the nozzle 


Wind tunnel 


Corrected to A.R. 0 


Condition : 
correction 

ae cy Cs CM, 00° C, Cao 
Wing spanning | —4 | —0-073 0:0122 | —0-0425 | Corrected by —3:°49 | —0-073 0-0116 
the jet —3 | —0:012 0-0112 | —0-0408 | Glauert and —2-92 | —0-012 0-O111 
ais, =='5 —1 +0-094 0:0113 | —0-0377 | Squire formulae. —1-66 0-094 0-0103 
0 0-147 0:0125 | —0-0359 | Equation —1-03 0-147 0-0101 
1 0-216 0:0146 | —0-0340 (6)-(9) | —0-51 0-216 0-0094 
3 0-327 0-0202 | —0-0301 Os71 0-327 0-0084 
5 0-429 0-0279 | —0-0262 1399 0-429 0-0075 
7 0-538 0:0389 | —0-0235 3°22 0-538 0- 0068 
9 0-665 0:0545 | —0-0241 4-33 0-665 0-0055 
11 0-773 0-0712 | —0-0215 5°56 0-773 0-0051 
Wing span 4-5 ft.| —4 | —0-074 0-0116 | —0-0403 | Equation —3*53 | —0-074 0-0112 
Pek. = 4:5 —3 | —0-024 0-0108 | —0-0392 (4) and (5) | —2-84 | —0-024 0-0109 
—1 0-093 0-0109 | —0-0371 —1-60 | +0-093 0-0103 
0 0-151 0-0118 | —0-0363 —0-98 0-151 0-0101 
1 0-209 0:0130 | —0-0347 —0-37 0-209 0-0098 
3 0-325 0-0170 | —0-0322 0-88 0-325 0-0100 
5 0-434 0-0230 | —0-0298 2°17 0-434 0-0092 
a 0-558 0-0324 | —0-0279 3°34 0-558 0-0097 
9 0-679 0:0436 | —0-0301 4-57 0-679 0-0099 
11 0-789 0:0566 | —0-0294 5°85 0-789 0-0111 
Wing span 4 ft. | —4 | —0-069 0-0116 | —0-0401 | Equation —3:54 | —0-069 0:0112 
AR: = 4 —3 | —0-016 0:0109 | —0-0393 (4) and (5) | —2-89 | —0-016 0-0109 
—1 0-096 0:0110 | —0-0368 —1-64 0-096 0-0102 
0 0-156 0:0109 | —0-0358 —1-05 0-156 0-0089 
1 0-202 0-0134 | —0-0344 —0-35 0-202 0-0101 
3 0-324 0-0183 | —0-0322 0-83 0-324 0-0097 
5 0-437 0:0252 | —0-0296 2-06 0-437 0-0096 
Z 0-556 0:0349 | —0-0284 3:26 0-556 0-0096 
9 0-684 0:0484 | —0-:0316 4-41 0-684 0-0100 
11 0-801 0:0638 | —0-0324 5-62 0-801 0:0113 
Wing span 3 ft. | —4 | —0-072 0-0123 | —0-0393 | Equation —3:44 | —0-072 0-0117 
ACR ta=3 —3 | —0-018 0-0117 | —0-0387 (4) and (5) | —2-86 | —0-018 0-0117 
—1 0-091 0:0118 | —0-0363 —1-72 0-091 0-0109 
0 0-146 0-0130 | —0-0355 —1-14 0-146 0-0107 
1 0-202 0-0149 | —0-0338 —0-58 0-202 0-0105 
3 0-296 0:0196" | —0-0311 0-67 0-296 0-0101 
5 0-402 0-0277 | —0-0285 1-84 0-402 0-0102 
7 0-519 0-0397 | —0-0270 2-92 0-519 0-0107 
9 0-639 0:0552 | —0-0304 3:97 0-639 0-0110 
11 0-756 0-0727 | —0-0335 5:06 0-756 0-0107 
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Results of traverses in the nozzle and the maximum section 


65-7 \? 


ak 2 aA wae ee 
OV ag pee (sep 


(S — A)/doV2 = 1-013 (S — A)/(P — S’) 


(P) S$ 00-087 (Pees) 


Traverse at J in nozzle Traverse at M in maximum section 
Tunnel 
condition | 

In. In. | LN 
above H-A P-S’ t0V? above H-A S-A ive 

C line C line 20 
No wing in tunnel 48 3-059 3-094 1-140 
24 3-057 3-365 3-321 45 3:058 3-093 1-139 
15 2-959 2-920 30 3-058 3-076 1-133 
6 2799 2-763 15 3-058 3-048 1-123 
0 B+ 7ol 27 Lo, 0 3-059 3-038 1-110 
— 6 2-774 2-738 —15 3-059 3-047 1-122 
—15 2995 2-956 —30 3-059 3-075 1-133 
—24 3-422 3-378 —45 3-059 3-107 1-144 
1 ft. chord wing 24 3-056 3°354 3-310 45 3-056 3-096 1-152 
spanning the jet 15 2-943 2-905, 30 3-057 3-079 1-145 
(at no lift) 6 2102 2-726 15 3-058 3-053 1-135 
0 2-725 2-690 0 3-058 3-043 1-132 
— 6 2 Tok 2-701 —15 3-058 3-050 1-134 
—15 2-955 2°OL7 —30 3-055 3-075 1-145 
—24 3-390 3-346 —45 3-055 3-104 1-154 
1 ft. chord wing 48 3-010 3-081 1-169 
spanning the jet 24 3011 3° 208 3-255 45 3-011 3-080 1-168 
ey ay 8 bt 15 3-011 2-888 2-850 30 3-014 3-063 1-162 
(L = 62 Ib.) 6 3-011 Of A 2-682 15 3-014 3-026 1-148 
0 3-066 2-671 2-636 0 3-011 2979 1-130 
— 6 3-010 2-641 2-607 —15 3-012 2-990 1-134 
—15 3-004 2-808 2-7/4 —30 3-012 3-006 1-140 
—24 3-009 3-309 3-266 —45 3-011 3-017 1-145 
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SYSTEM OF AXES 


Symbol ae y | i 
Axes cach eciatm longitudinal lateral normal 
BOsuLT forward starboard downward 
direction : 
Force Symbol x 4 Z 
; Symbol — We M any AN. 
: co Designation rolling» pitching yawing 
Angle of 
Rotation ae daowtia : | : 4 
Velocit Linear ui | v : w 
yi Angular p | q | r 
ee Sale A 
~ Moment of Inertia fa | B C 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. | 


Components of angular velocity and moment are positive in the cyclic order 
y to z about the axis of x, z to x about the axis of y, and x to y about the axis of z. 


The angular movement of a control surface (elevator or rudder) 1s governed 
by the same convention, the elevator angle being positive downwards and the rudder 
angle positive to port. The aileron angle is positive when the starboard aileron is 
down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 


The symbols for the control angles are :— 
€ aileron angle : 
n elevator angle 
ny tail setting angle 
¢ rudder angle 
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“Tests of Aerofoil N.A.C.A. 23012 in the Compressed 
Air Tunnel 


By 


Delis. WittiAMs, BocwA:P UR Ac SiwAt, Hi Betz AND E. Smytu, B.Sc., A.F.R.Ae.S. 
of the Aerodynamics Department, N.P.L. 


Reports and Memoranda No. 1898 
2nd .“fanuary, 1937 


Summary —N.A.C.A. 23012 is 12 per cent. thick and has a maximum camber of 2 per cent. at 15 per cent. aft of 
the leading edge. It is designed to have a very small pitching moment and a high maximum lift. The profile is 
superior to R.A.F.34 in that the minimum profile drag is lower, occurring at C, = 0-2 for N.A.C.A. 23012, whereas 
at.occurs. at) G,.—= 0 tor K.A.F 34. 


The section has also been tested in the Variable Density Tunnel and in the Full Scale Tunnel of the N.A.C.A. 
Results from the three tunnels are compared. Agreement is good at high Reynolds numbers, taking the “‘ effective ”’ 
Reynolds number for the Variable Density Tunnel. The very high maximum lift coefficient (1-61) observed in the 
V.D.T. was not obtained in the C.A.T., where the highest value was 1-49. The maximum lift values obtained in 
the C.A.T. and F.S.T. are in good agreement. 


Later measurements on the aerofoil with a 15 per cent. split flap at 90° have been included. 


As an extension of the study of a family of aerofoils in the Variable Density Tunnel, a series 
was developed having the maximum camber forward of the usual position. The most 
promising section, N.A.C.A. 23012, was tested! over a range of Reynolds numbers in the 
Variable Density Tunnel (V.D.T.) and in the Full Scale Tunnel (F.S.T.) of the National 
Advisory Committee for Aeronautics. It is claimed that this section is superior to sections in 

_ common use, in that it has a high maximum lift, a low minimum drag and a very small 
4; movement of the centre of pressure. In order to test its properties independently, the aerofoil 
€ has been tested in the Compressed Air Tunnel (C.A.T.) over. the complete range of Reynolds 
,, number available. The results are of particular interest, as they afford an opportunity of a 
ve. direct comparison of results from the Compressed Air Tunnel with those from the two American 
*~ tunnels at high Reynolds numbers. 

Pee 


The section N.A.C.A. 23012 has a maximum camber of approximately 2 per cent. (actual 
1-85 per cent.) Situated 15 per cent. aft of the leading edge with the usual N.A.C.A. 12 per cent. 
thick fairing, the maximum thickness occurring at 0-3c. The centre line is straight from 
a point just behind the position of maximum camber to the tail. The profile and centre line 
are shown in Fig. | and the ordinates of the section are given in Table 1. 


17 M4s 
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The model tested in the Compressed Air Tunnel was of hiduminium, with hand-finished 
surface. It was rectangular, of aspect ratio 6; the chord was 8 in. and the span 48 in. 
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Lift, drag, pitching moment and centre of pressure coefficients are given in Table 2, corrected 
for tunnel interference. The profile drag coefficient has been included, calculated from the 


formula C,= Co, + 0-0555 C,?. 


The most important characteristics have been plotted against Reynolds number in the figures, 
together with the corresponding results from the American wind tunnels. 


In general, the results substantiate the claims made for the aerofoil. 


Values of lift, drag and pitching moments on the aerofoil with a 15 per cent. split flap at 90°, 
fixed 0-15c forward of the trailing edge, have been given in Table 3. These measurements 
were made later than those on the plain wing. The results are taken from Table 7 of A.R.C. 
report 4191.* 


Comparison with R.A.’.34.—In N.A.C.A. 23012, a small pitching moment was obtained 
by placing the maximum camber far forward. In R.A.F.34, the same effect was achieved by 
giving the aerofoil a reflexed trailing edge. The series of aerofoils tested? in the Variable 
Density Tunnel included some with reflexed trailing edge and having the maximum camber 
far forward, but it was found that the series without reflexed trailing edge had the smaller 
minimum profile drag, hence the choice of this particular section for further test. 


The following table gives a comparison of certain characteristics of the two aerofoils, 
obtained’ in the Compressed Air Tunnel— 


N.A.C.A. 23012 R.A.F.34 
Reynolds number .. ‘a 43 $i 7:05 x 108 TAT 10S 
max: an ie a 4) ew. 1-49 1-41 
Cp, min. bid a oe by ee 0-0074 0-0080 
L/D max. a3 He :5 aif at 24-4 23°4 
C, at L/D max. % 5 te 5 0-35 0-35 
Coat G,=0 .: ye ie os 43 0-0078 0- 0080 
AC, |dag % i _ fe sip 0-100 0-103 
CG, max./Cp, min... 3, ¥ £ 191 176 
ing “tee ; € px ts —0-0093 —0-0041 
C. P, forward position A ig A 0-250 0-250 
C.P. at 1/4 Cmax. .. 7 5; 0-275 0-251 


The advantages of N.A.C.A. 23012 over R.A.F.34 are that it has a higher maximum lift, 
a lower minimum profile drag, a higher maximum L/D and a higher “ speed-range index ”’ 
(C;max./Cp,min.). The minimum profile drag occurs at zero lift for R.A.F.34, for the 
American section it occurs at about C; = 0:2, R.A.F.34 is superior in that it has a smaller 
value of Cm, though this coefficient is small for both aerofoils. R.A.F.34 also has a smaller 
C.P. movement at low values of C,. 


Comparison with American Tunnels.—At the outset it is necessary to explain the procedure 
in recent American reports!” for the presentation of results. Owing to the turbulence in the 
tunnel, the maximum lift values measured in the Variable Density Tunnel have been much 
higher than those obtained in less turbulent tunnels and on full scale, so a “ turbulence factor ”’ 
has been introduced. Multiplying the actual Reynolds number of the test by the turbulence 
factor, an “ effective Reynolds number ”’ is obtained, and it is this value of R that is used in 
the diagrams. The tentative turbulence factor adopted for the Variable Density Tunnel 


* Tests in the Compressed Air Tunnel on N.A.C.A.0015 with and without a flap and on other aerofoils with a flap.— 
R. Jones, A. H. Bell and A. F. Brown. August, 1939. (Unpublished.) 
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is 2-64, that for the Full Scale Tunnel is 1-1. The lift values have been applied directly at the 
effective Reynolds number but a correction has been applied to the drag results to allow for 
the change of Reynolds number. It is assumed that the measured drag coefficients are too 
high and may be diminished by an amount equal to the decrease in the turbulent skin friction 
coefficient due to moving from one Reynolds number (test) to a higher Reynolds number 
(effective). At the highest Reynolds number for the Variable Density Tunnel results, the 
correction to profile drag coefficient is 0-0011. 


The results for the Full Scale Tunnel are given in Tables in the American report, but the 
Variable Density Tunnel results are given only in small diagrams, except for isolated data for 
comparison with other results. In all the American diagrams, points are plotted at the 
effective Reynolds number. 


Maximum lift has been plotted in Fig. 2. The Compressed Air Tunnel values have been 
plotted at the actual Reynolds number, the Full Scale Tunnel results at the effective Reynolds 
number—the points should be moved slightly to the left for the actual Reynolds number—and 
the Variable Density Tunnel values have been plotted at both test and effective Reynolds 
number. There is good agreement between the results from the Compressed Air Tunnel and 
the Full Scale Tunnel, the results from the Variable Density Tunnel agree with these when 
plotted at the effective Reynolds number, except for high values of R, when the V.D.T. results 
are considerably higher than those of the C.A.F. The highest values obtained in the three 
tunnels are: 1-61 at R=8-05 x 108 (effective) in the V.D.T., 1-49 at R= 7-05 x 108 in 
the C.A.T., and 1-46 at R=4:0 x 108 and R=4-5 x 10° in the F.S.T. There is no 
similarity between the results on maximum lift from the Compressed Air Tunnel and the 
Variable Density Tunnel when the results from the latter are plotted at the test Reynolds number. 


The drag coefficient at zero lift is shown in Fig. 3. The values for the V.D.T. at actual 
Reynolds number have been obtained by estimating the corrections applied and adding them 
to the values given at the effective Reynolds numbers. The information given is sufficient 
for determining the position of only one of the four points accurately, but any errors-in the 
positions of the other three points must be small. For the highest Reynolds numbers, the 
results agree for all three tunnels, taking the corrected values for the V.D.T. at the effective 
Reynolds number. The C.A.T. results are surprisingly low at low values of Reynolds number. 


Three other quantities are given in Fig. 4, for as many points from the three tunnels as are 
available. They are minimum profile drag, moment coefficient at zero lift and the slope of 
the lift curve at zero lift (dC,/do, at C, = 0). 


Finally, the following table gives a comparison of data from the three tunnels. 


NALA. 20012 
Caer Bisuik Vi. Dele 
Test R.N. 
4°47 X 10° | 3°1 K 10% | 4-50 x 10° | 3:09 x 10% | 3:09 x 10°f | 1-29 x 10% 
GC Miaxa ear a as ee 1-44 1-40 1-50§ 1-40 161 1-43 
a, degrees a EM es —1:-0 —1-0 —1-2 —1-2 —1:2 —1-2 
AC, /dxg aye a 0-100 0-098 0-101 0-099 0-104 0-102 
Cp,min. .. a i: a 0-0072 0-0074 0-0069 0-0072 0-0074 0-0084 
Cing <3 a 3 a —0-0100 —0-0100 —0-0087 —0-0053 —0-0075 —0-0065 
C,max./Cp,min. .. ie oe te ou 189 208]| 195 217 170 
Gk riiine < Fie Pe aie 3 0-0077 0-0077 0-0078 0-0081 0-0077 0-0086 
L/D max... He ts - 24-1 22:8 25 - O|| 24-1 23-8 22-5 
eats hy tix, Bae a 0-35 0-35 0-3 0-30 0-36 0-40 
C.P. forward pos. fs ai 0-249 0-248 0-250 0-250 0-256 —- 
GPratcy= 0:3 a “4 0-264 0-262 0-260 0-250 — _— 
* Results interpolated from those at R = 2:36 x 108 and R = 4-47 x 10°. + Effective R.N. of 3-40 x 108. 
t Effective R.N. of 8-16 x 108. § Highest value given in Tables is 1-46. |p As given. 
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The values given for the Variable Density results at the test Reynolds number of 3-09 x 108 
are taken from the American report! giving results from the V.D.T. and the F.S.T. An 
American Technical Note* giving results on a series based on N.A.C.A. 23012 with various 
thickness ratios gives a higher value of maximum lift coefficient (1-68) and a lower minimum 
drag coefficient (0:0072) at an effective Reynolds number of 8-37 x 106. These values have, 
however, had corrections applied for turbulence and for tip effect, the latter correction increases 
maximum lift by 3 per cent. 
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TABLE 
N.A.C.A. 23012 


Dimensions 1n Terms of Chord 


Distance from Distance above Distance below 
Leading Edge Chord Chord 
0 0 0 

0:0125 0-0267 0:0123 
0-025 0-0361 0:0171 
0-050 0:-0491 0-0226 
0-075 00-0580 0-0261 
0:10 0:0643 0-0292 
0-15 0:-0719 0-0350 
0-20 0-0750 0:0397 
0-25 0-0760 0-0428 
0-30 0:0755 0:-0446 
0-40 0:0714 0:0448 
0-50 0:0641 0-0417 
0-60 0:-0547 0-0367 
0-70 0-0436 0-0300 
0-80 0-0308 0-0216 
0-90 0-0168 0-0123 
0-95 0-0092 0-0070 
1-00 0 0 


L.E. Rad. : 0-0158 


Aerofoil N.A.C.A. 23012. Aspect Ratio 6 


TABLE 2 


R='0-315 x 10% P = 1 atmos, 40V? = 6:62,, V = 74:9 f./s. 


a? CG, Cs Can C,,(c/4) G.P: 
— 4-5; —0-233 0-0176 0-0146 —0-0211 0-160 
— 2-4, —0-084, 0-0124 0-0120 —0-0156 +0-066 
— 1-4 —0-022 0-0106 0-0106 —0-0108 —0- 236 
— 0:3; +0-042, 0-0101 0-0100 —0-0069 +0-412 
+ 0-6, 0-113 0-0104 0- 0097 —0-0030 0-276 
LF 0-188, 0-0126 0-0106 —0-0008 0-254 
3°75 0-366 0-0197 0-0123 —0- 0066 0-268 
6:8, 0-617 0-0351 0-0140 —0-0123 0-270 
9-9, 0-812 0-0532 0-0166 —0-0059 Q- 257 
13-0 1-000 0-0770 0:-0215 +0-0021 0-248 
15-0; 1-086 0-0951 0: 0297 0- 0080 0-243 
16-1 1-104 0-106 0-038 +0-0079 0-243 
16-3* 0-893 ~- —- —0-0408 — 
17-1, 1-080 0-131 0-066 —0-0075 0-257 
17-3,* 0-793 0-227 0-192 —0-0711 0-336 
18*2, 1-058 0-154 0-092 —0-0208 0-270 
18-4* 0-744 0-249 0-218 —0-0769 0-348 
20-5; 0-657 0-280 0-256 —0-0805 0-362 
23-7, 0-642 0-326 0-323 —0-0850 0-368 
* Decreasing incidence. 
Stalling angle 16-0° 
Re=,0-802-. 10% P29 atmos.) FoVa— los) Vj == 66: 6,1. /s: 

— 4:6 —0: 241 0-0131 0- 0099 \ —0-0189 0-172 
— 2-4, —0-095 0- 0096 0-0091 —0-0139 +0: 103 
— 1-4 —0:021 0: 0088 0-0088 —0-0115 —0-293 
— 0:3, +0-044 0-0077 0-0076 —0-0077 +0-425 
+ 0:6; 0-113 0-0078 0-0071 —0-0045 0-290 
1:7 0-186, 0-0094 0-0075 —0-0017 0-259 
3°8 0-334 0-0141 0-0079 +0-0015 0-246 
6-8, 0-576 00277 0-0093 —0-0008 0-251 
10-0 0-795 0-0477 0-0126 -+0-0006 0-249 
12-0 0-932 0- 0620 0-0138 0: 0033 0-246 
14-1; 1-065 0- 0800 0-0170 0-0058 0-244 
15-1, I-iil 0-0908 0-0224 0- 0070 0-243 
16:1, 1-162 0-101 0-026 +0-0091 0-242 
17-2, 1-115 0-134 0-065 —0-0148 0-263 
19-3, 1-059 0-176 0-114 —0-0344 0-282 
21-4, 1-005 0-225 0-169 —0-0548 0-304 


TABLE 2—continued 
R= 1-26 210°, P= 4228 atmos:i40 V7—"2o" 4) Vee COE. e: 


* Decreasing incidence. 
Stalling angle 19-4° 


* 


ot C,, Cp Cho C,,(c/4) CrP: 
— 4:6 —0Q-249 0-0123 0-0089 —0-0172 0-181 
— 2:55 —0-100 0-0091 0-0085 —0-0130 +0-120 
— 1:4 —0-022 0-0082 0-0082 —0-0112 —0-255 
— 0:3; +0-049, 0-0080 0-0079 —0-0089 +0-430 
+ 0-6, 0-121 0-0078 0-0070 —0-0057 0-297 
1-7, 0-191 0-0089 0- 0069 —0-0034 0-268 
3°8 0-344 0-0137 0-0071 —0-0007 0-252 
6-9 0-576 0-0261 0-0077 +0-0018 0-247 
10-0; 05799 0-:0461 0-0107 0-0036 0-245 
12-1 0-943 0-0618 0-0125 0-0044 0-245 
14-1; 1-073 0-0804 0-0165 0-0050 0-245 
16-2, 1-193 0-101 0-022 0-0063 0-245 
17-2, 1-244 0-113 0-027 +0-0065 0-245 
18-3; 1-130 0-158 0-087 —0-0223 0-270 
20-4; 1-045 0-202 0-141 —0-0465 0-294 
22°65 0-968 0-251 0-199 —0-0639 0-314 
Stalling angle 18-0° 

R= 2;36 810°, P= 7°84 atmos, 40V*— 0451.) Vie 7/8 Gis 

— 1:4, —0-032 0- 0080 0-0079 —0-0108 —0-086 
— 0:3, +0-041; 0-0076 0-0075 —0-0090 +0-467 
+ 0:7 0-118 0-0084 0-0076 —0-0071 0-310 
1°75 0-192 0-0099 0-0079 —0-0052 0-277 
3°9 0-345 0-0157 0-0091 —0-0024 0-257 
7-0, 0-583 0-0280 0-0091 +0-0001 0-250 
10-2, 0-801 0-0478 0-0122 0-0019 0-248 
12-3 0-959 0-0650 0-0140 0-0024 0-248 
14-4; 1+ 102 0-0844 0-0171 0-0019 0-248 
16-5 1-242 0-106 0-020 0-0014 0-249 
17-5; 1-295 0-118 0-025 0-0013 0-249 
18-6 1-347 0-130 0-029 +0-0010 0-249 
18-6* 1-148 0-165 0-092 —0-0308 0-277 
19°65 1-083 0-195 0-130 —0-0509 0-297 
20-6; 1-043 0-215 0-155 —0-0566 0-304 
22°75 0-931 0-258 0-210 —0-0736 0-327 


23-0 


C, 


—0-034 
+0:046 

0-128 

0-212 
. 0-374 

0-616 

0-852 
-007 
-156 
*zol 
352 
-410 
-140 
-052 
-005 
0-898 


ee ee ee 


* Decreasing incidence. 
Stalling angle 20-1° 
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Ke= 57000 10% P = 13: Latmos,, 4aV*— 125: 


TABLE 2—continued 
R’='4"47 x 10°) P =='14-83 atmos., $6V" = 102°5, V = 78:7 f./s. 


C,,(¢/4) C.P. 
—0-0108 —0-066 
—0-0090 +0:446 
—0-0073 0-307 
—0-0057 0-277 
—0:0029 0-258 
+0-0002 0-250 

0-0009 0-249 
+0-0004 0-250 
—0-0003 0-250 
—0:0017 0-251 
—0- 0023 0-252 
—0-0031 0-252 
—0-0380 0-283 
—0-0558 0-302 
—0-0626 0-311 
—0-0803 0-336 


2., V = 78-2 f./s. 


—0-035; 
+0- 046, 
0-131 
0-214, 
0-382 


0-890 


* Decreasing incidence. 
Stalling angle 20-1° 


0: 0077 
0: 0078 
0: 0086 
0-0103 
0-0158 
0: 0301 
0-0516 
0- 0690 
‘0913 
-113 
-128 
-141 
-187 
-211 
232 
-275 


SSS Sr ar=> 


0: 0076 
0-0077 
0: 0076 
0:0077 
0-0077 
0- 0088 
0-0101 
0-0107 
0-0161 
0-018 
0-024 
0-028 
0-122 
0-153 
0-177 
0-231 


—0-0105 
—0- 0086 
—0-0071 
—0-0055 
—0-0030 
+0-0001 
+0-0002 
—0- 0006 
—0-0021 
—0-0039 
—0-0041 
—0-0045 
—0-0420 
—0-0557 
—0-0620 
—0- 0806 
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TABLE 2—continued 
R = 7-05 x 108, P = 22-6 atmos., }9V? = 160-5, V = 78-7 f./s. 


on? CG, Cp Coo C,,(c/4) 2.P, 
— 1:5 —0- 036, 0- 0080 0-0079 —0-0101 —0-022 
— 0:3; +0: 047, 0-0077 0-0076 —0-0083 +0-426 
+ 0:8 0-133, 0- 0084 0-0074 —0-0067 0-300 
1-9; 0-220 0-0101 0-0074 —0-0051 0-273 
4-2. 0-389 0-0161 0- 0083 —0-0023 0-256 
7°55 0-638 0-0309 0- 0083 —0-0002 0-250 
11-0 0-884 0-0545 0-0112 —0-0005 0-251 
13-2 1-048 0-0737 0-0127 —0-0020 0-252 
15-4, 1-198 0- 0969 0-0172 —0:0036 | 0-253 
17-6 1-343 O2a9 0-019 —0-0054 0-254 
18-7 1-398 0-134 0-025 —0-0060 0-255 
19-7, 1-460 0-147 0-029 —0-0065 0-255 
20:2, 1-485 -— ~- — mas 
20-5; 1-024 0-219 0-161 —0-0599 0-308 
21-5, 0-980 0-236 0-183 —0-0647 0-315 
23:4, 0-875 0-285 0-242 —0-0786 0-337 
Stalling angle 20-2,° 
AAD 
N.A.C.A. 23012 with 0:15c split flap at 90°, fixed at 0:15c from trailing edge 
Asbect Ratio 6 
R= 0-315 49108, P = 1 atmosipipV* =6 71/2 R= 0°81. x 10°, P = 3-0 atmos mey tom 
V = 76°12.)/s. V = 69°61./s 
ra | 6 Cy ‘ a? CE | C8 | C,, 
—6:9, 0-449 0-168 —0-166 —7-0 0-447 0-164 —0-163 
—5:1 0-616 0-182 —0-180 —5-1 0-606 7S —0-174 
—1-1,; 0-897 0-212 —0-183 —1-1, 0-923 0-212 —0-186 
+2:°9, 1-149 0-244 —0-180 +3:0 1-202 0-245 —0-185 
Weds 1-402 0-283 —0-174 7-1, 1-471 0-283 —0-181 
10°1, 1-602 0-319 —0-174 10-2 1-656 0-320 —0-179 
12 1-713 0-344 —0-176 \ places 1-769 0-346 —0-178 
12°4.% 1-262 0-447 —0-235 13°+.3, 1-839 0-361 —0-:179 
13-3 1-772 0-359 —0-177 14-3; 1-898 0-376 —0Q-181 
MO Th 1-136 0-475 —0- 234 15-4 1-957 0-388 —0-180 
14-3 1-838 0-368 —0-180 135° 7* 1-366 0-507 —0-232 
14-7, 1-086 0-501 —0- 237 16-4 2-018 0-405 —0-183 
15-8, 1-058 0-523 —0- 237 16-7,;* 1-248 0-670 —0-229 
17-9; 1-069 0-587 —0-249 17-4 2:072 0-417 —0-186 
17-8, 1-162 0-592 —0- 247 
* Decreasing incidence. 19-0 1-146 0-633 —(0: 256 
Stalling angle 14-7° 20-0 1-156 0-666 —0-262 
22:1; 1-152 0-722 —0-269 


* Decreasing incidence 
Stalling angle 17-8° 


: ; 
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TABLE 3—continued 


R= 1-27 X 10°, P = 4:4 atmos., toV* = 25°3 
V = 70-1 f:/s 
a° te, Cs | Ce 
| 
—7:0 0-445 0-168 —0-166 
—5:1 0-602 0-181 —0-173 
—1-1, 0-929 0-214 —0- 187 
+3:-0,; 1-206 0-248 —0-186 
7°2 1-468 0-290 —0-184 
10-2, 1-665 0-325 —0-182 
12:3 1-785 0-351 —0-181 
14-4 ‘1-910 0-378 —0-181 
15-5 1-961 0-395 —0- 182 
15-6" 1-628 0-423 —0- 204 
16-5 2-037 0-406 —0-183 
16°75 1-494 0-454 —0Q-211 
17-5 2-098 0-421 —0-184 
17° 3* 1-360 0-580 —()- 249 
18:5, 2-153 0-439 —0- 187 
_18-9,;* 1-115 0-604 —0-245 
19-9, 1 +242 0-645 —0-251 
22:1 1-235 0-705 —0- 258 
* Decreasing incidence 
Stalling angle 18-9° 


Rota 47a, 10?) Pi 
V= 


14-6 atmos., 4¢V2 = 101-7 
78-8 f./s. 


n° oF Cp ee 
| 

—7-1, 0-450 0-169 —0-169 
—5-2 0-600 0-181 —0-174 
—1-0, 0-909 0-209 —0-177 
+3:3 1-211 0-245 —0-178 

7:6 1-489 0-289 —0-177 
10-7; 1-688 0-327 —0-177 
12-9, 1-814 0-355 —0-177 
15-1 1-933 0-383 —0-177 
7+ 2, 2: 067 0-414 —0-180 
07:08 1-536 0-471 —0-215 
18-3 2-134 0-430 —0-181 
18-0* 1-469 0-520 —(- 230 
19-4 2-192 0-446 —0-183 
19-0,* 1-322 0-535 —0- 235 
20-4, 2-244 0-465 —0-184 
20-1* 1-262 0-556 —0- 234 
21 - T* 1-218 0-585 —0-242 
23°1, 1-212 0-712 —0- 268 

* Decreasing incidence 


Stalling angle 21-4° 


I= 2-47. x 10° P= 8-0 atmos., ¢eV72 = 53°9 
= 76°71 


= 76-7. 64s. 

a? Ce Cp Cy 
—7:0; 0-438 0-167 —0-166 
—5:1; 0-591 0-179 —0-170 
—1-] 0-909 0-212 —0-179 
+3-1 1-203 0-248 —0-181 

7°35 1-479 0-291 —0-180 
10-4, 1-675 0-327 —0-180 
12-5, 19797 0-352 —0-178 
14-7 1-923 0-384 —0-179 
15:7 1-984 0-395 —0-179 
16:7, 2-045 0-411 —0-180 
16-8* 1-552 0-458 —0-214 
17:8 2-094 0-425 —0-179 
17-8* 1-432 0-483 —0-223 
18-9 2°175 0-440 —0-182 
18:94 1-299 0-515 —(-224 
19:9 2:215 0-456 —0-184 
20° 0," 1-258 0-545 —0- 235 
yA Nea ke P27 0-577 —0- 238 
221 1-226 0-637 —0- 256 
23°1, 1-246 0-697 —0-272 

* Decreasing incidence 
Stalling angle 20: 4° 


R = 5-43 x 108, P = 18-1 atmos., }oV2 = 124 


V = 78-4 f./s. 
a? ‘of Gs Gy 

zy. 9 0-453 0-171 —0-171 
—5-2 0-599 0-180 —0-173 
—1:0 0-917 0-210 —0-178 
43-4 1-221 0-245 —0-178 
7-7; 1-501 0-290 —0-178 
10-9, 1-705 0-329 —0-177 
13-0; 1-834 0-358 —0-177 
15:3 1-967 0-387 —0-178 
17°45 2-076 0-419 —0-180 
17-05% 1-522 = a 
18-5 2-137 0-436 —0-182 
18-0;* 1-477 0-511 —0-229 
19-6, 2-210 0-453 —0-185 
19-1* 1-329 0-545 —0-236 
20-7 2-263 0-467 —0-185 
20+ 1* 1-258 0-563 —0-234 
21-1,* 1-204 0-600 —0+245 
22-15 1-216 0-656 —0+259 
24-05 1-232 0-750 —0:284 


* 


Decreasing incidence 
Stalling angle 21-7° 
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® CAT. Test RN. 0 V.D.T. Test R.N, 
x FST. EFP. RIN A V.0.T. EFF. RN. 


L.-T. fa 
Test RN. 
Ve 


Ty 


V.D.T. 
EFFective RN. 


cAr. 
Test R.N. 


NACA. 23012. 


oT t 


os 


RxlO&. 
MAXIMUM LIFT AGAINST REYNOLDS NUMBER 


lables 24 
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®@ CAT. Tesl R.N. Oo V.O-7T. Test R.N. 
Xe fal. Creer, ; A V.BT. EFF. R.N. 


O-4 os 0-7 ie 2 3 q 5 7 | Ta} 
DRAG COEFFICIENT AT ZERO}LIFT » 


Hicso 


2 
COMPARISON OF ¢ omin.s Cre AND (dc, /da,) IN CAT. AND 


AMERICAN WINO ON Ee 


o's os OT U Rxlo~> a 3 4 5 T 


(60168) Wt.6/R8 375 8/42 Hw. G.377/1 


SYSTEM OF AXES 


ss Symbol | 


| 

| : : x : OY z 

Axes | es oe longitudinal | lateral . normal 
dean forward starboard | downward 

Force Symbol ees w: Z 

Moment Designation rolling pitching yawing 
Angle of | | 

Rotation ial > ‘ if 

TBE t Linear us v w 

coe y Angular ? | q 4 

By C 


Moment of Inertia A 


Components of linear velocity and force are positive in the positive direction 
of the corresponding axis. 


Components of angular velocity and moment are positive in the cyclic order 
¥ to z about the axis of x, z to x about the axis of y, and x to y about the axis of 2. 


The angular movement of a control surface (elevator or rudder) is governed 
by the same convention, the elevator angle being positive downwards and the rudder 


i angle positive to port. The aileron angle is positive when the starboard aileron is 


down and the port aileron is up. A positive control angle normally gives rise to a 
negative moment about the corresponding axis. 
The symbols for the control angles are :— 
€ aileron angle 
n elevator angle 
nr tail setting angle 
¢ rudder angle 


Publications of the 
Aeronautical Research Committee | ¢€ 


TECHNICAL’ REPORTS »-OF °° THE, “AERONAUTICAL 
RESEARCH COMMITTEE— 


1933-34 Vol. I. Aerodynamics. {1 $s. 
Vol. II. Structures, Engines, Instruments, etc. {1 10s. 
1934-35 Vol. I. Aerodynamics. £2. 
Vol. Il. Seaplanes, Structures, Engines, Materials, etc. 
£2. 
1935-36 Vol. I. Aerodynamics. £1 Ios. 
Vol. II. Structures, Flutter, Engines, Seaplanes, etc. 
PUTO. 
1936 Vol. I. Aerodynamics General, Performance, Air- 
- screws, Flutter and Spinning. {2. 
Vol. II. Stability and Control, Structures, Seaplanes, 
Engines, etc. £2 10s. 
1937 Vol. I. Aerodynamics General, Performance, Air- 
screws, Flutter and Spinning. 2. 
Vol. Il. Stability and Control, Structures, Seaplanes, 
Engines, etc. £3. 


ANNUAL REPORTS OF THE AERONAUTICAL RESEARCH 


COMMITTEE— 
1933-34 Is. 6d. 
1934-35 1s. 6d. 
April1j¢993'5 to: December. 21) /4930. 7-4 
1937 Of 
1938 Is. 6d. 


INDEX TO). THE “4 BCHNICAL (REPORTS Or, Wik 
ADVISORY COMMITTEE ON AERONAUTICS— 


1909-1919 Reports and Memoranda No. 1600. 8s. 


Prices are net and postage extra. 


His Majesty’s Stationery Oihice 
London W.C.2: York House, Kingsway 


Edinburgh 2: 120 George Street Manchester 2: 39-41 King Street 
Cardiff: 1 St. Andrew’s Crescent Belfast: 80 Chichester Street 


or through any bookseller. 


S.0O. Code No. 23-1898 
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